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ABSTRACT 

In  this  article  we  present  a  self-contained  account  of  the  recent 
theory  of  micropolar  elasticity.  Micropolar  elastic  materials,  roughly 
speaking,  are  the  classical  elastic  materials  with  extra  independent  degrees 
of  freedom  for  the  local  rotations.  These  materials  respond  to  spin  inertia 
and  body  and  surface  couples,  and  as  a  consequence  they  exhibit  certain  new 
static  and  dynamic  effects,  e.g.,  new  types  of  waves  and  couple  stresses. 

The  theory  is  fully  deterministic  as  against  the  background  of  the  recently 
popular  indeterministic  couple  stress  theory  (cf.  Art  23).  The  mechanics 
of  certain  classes  of  materials  with  fiberous  and  elongated  grains  (e.gt, 
dumbbell  types  '1  grains)  represents  a  potential  field  of  application  of  the 
theory . 

The  geometry  of  deformation  and  its  measures  are  introduced  on  a  more 

effects 

general  background  of  materials  exhibiting  granular  and  microstructure^(micro- 
morphic  materials).  Various  types  of  microstrains  and  microrotations  are 
discussed.  Compatibility  conditions  for  the  micropolar  strains  are  derived 
The  kinematics  of  strains,  microstrains  and  rotations  are  presented  Basic 
laws  of  motion,  conservation  of  mass,  conservation  of  microinertia,  balance 
of  momentum,  balance  of  moment  of  momentum,  conservation  of  energy  are  postu¬ 
lated  and  their  local  forms  are  obtained.  The  thermodynamics  of  micropola- 
solids  is  formulated  and  the  consequences  of  the  entropy  inequality  are 
discussed.  Constitutive  equations  are  found  for  the  linear  theory  of  micro- 
polar  elasticity.  The  basic  field  equations  and  initial  and  boundary  conditions 
are  given. 

The  indeterminate  couple  stress  theory  is  shown  to  result  as  a  special 
case  of  the  theory  when  the  motion  is  constrained  Several  static  and  dynamic 
problems  are  solved  to  reveal  seme  of  the  new  physical  phenomena  exhibited  by 
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the  theory.  These  include  propagation  of  waves  in  infinite  micropolar 
elastic  solids,  reflection  of  various  types  of  micropolar  waves  in  a 
half  space,  surface  waves,  the  stress  concentration  around  a  circular 
hole  in  a  tension  field,  and  force  and  moment  singularities  in  an  infinite 
solid.  The  Papkovitch  and  Galerkin  representations  are  presented. 

The  article  is  based  mostly  on  the  works  of  Eringen  and  his  co-workers 
published  during  the  past  several  years.  Many  parts,  however,  contain  new 
compositions  and  several  other  results  are  presented  for  the  first  time. 
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1.  INTRODUCTION 

Classical  continuum  mechanics  is  based  upon  the  fundamental  idea 
that  all  material  bodies  possess  continuous  mass  densities,  and  that  the 
laws  of  motion  and  the  axioms  of  constitution  are  valid  for  every  part 
of  the  body  no  matter  how  small  they  may  be.  Accordingly,  a  small  volume 
AV  enclosed  within  a  surface  AS  possesses  a  mass  density  p  defined  by 

(1.1)  _  lim  Am 

P  '  AV  -►0  AV 

where  Am  is  the  total  mass  contained  in  AV.  Here  p  is  independent  of  the 

size  of  AV  and  depends  only  on  the  position  vector  x  of  a  point  in  AV  and 

time  t.  Consider  the  following  experiment  for  the  measurement  of  p:  The 

mass  density  of  a  homogeneous  material  may  be  calculated  approximately  by 

weighing  a  large  number  of  pieces  having  different  volumes  and  calculating 

the  ratio  Am/AV  for  each  piece.  If  the  resulting  numbers  p  are  plotted 

against  AV,  one  finds  that  this  ratio  is  nearly  constant  when  AV  is  greater 

* 

than  certain  critical  volume  AV  and  begins  to  show  dependence  on  AV  when 
*  * 

AV  <  AV  .  The  size  of  AV  depends  on  the  constitution  of  the  material.  As 
AV  approaches  zero,  this  dependence  becomes  violent,  Fig.  1.1.  This  situa¬ 
tion  is  well-understood  when  we  remember  the  granular  and  molecular  nature 
of  materials.  The  classical  continuum  theory  may  therefore  not  represent 

a  good  mathematical  model  for  the  approximation  of  a  physical  theory  in  the 
* 

range  AV  <  AV  . 

The  loss  of  accuracy  in  classical  continuum  mechanics  may  stem  from 
another  important  reason.  If  the  response  of  the  body  is  sought  to  an  external 
physical  effect  in  which  the  length  scale  is  comparable  to  the  average  grain 
or  molecular  size  contained  in  the  body,  the  granular  or  molecular  constituents 
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of  the  body  are  exited  individually.  In  this  case,  the  intrinsic  motions 
of  the  constituents  (microelements)  must  be  taken  into  account.  This  point 
becomes  clear  especially  in  connection  with  the  propagation  of  waves  having 
high  frequencies  or  short  wave  lengths.  When  the  wave  length  is  of  the  same 
order  of  magnitude  as  the  average  dimension  of  the  microelements,  the  intrinsic 
motions  of  the  microelements  of  AV  with  respect  to  the  center  of  mass  of  AV 
can  affect  the  outcome  appreciably.  This  situation  previals  in  practical 
applications  when  the  material  under  consideration  is  a  composite  material 
containing  macromolecules,  fibers,  and  grains.  For  such  materials,  the  critical 
volume  AV*  is  of  the  order  of  magnitude  of  the  cube  of  a  fraction  of  an  inch 
to  several  inches.  Solid  propellant  grains,  polymeric  materials,  and  fiberglas 
are  but  a  few  examples  for  such  materials. 

Another  example  is  the  anomalous  behavior  of  blood  when  flowing  through 
capillaries.  Blood  consists  of  a  fluid  (plasma)  in  which  are  dispersed  elements 
of  microscopic  dimensions  (corpuscles).  When  blood  flows  through  capillaries 
whose  diameter  is  comparable  to  that  of  the  corpuscles,  the  flow  characteristics 
differ  from  those  for  large  vessels.  As  a  further  example,  experiments  have 
shown  that  the  resistance  of  a  solid  to  surrounding  fluid  can  be  reduced  by  as 
much  as  1/3  when  a  minute  amount  of  additives  are  cast  into  fluid. 

It  may  be  conjectured  that  a  rational  treatment  of  surface  tension,  micro¬ 
cracks,  microfracture,  and  the  mechanics  of  granular  media  and  composite  materials 
ultimately  will  have  to  be  based  on  the  theory  of  microcontinua.  The  nature  of 
experimental  work  on  the  properties  of  such  materials,  no  doubt,  will  be  affected 
with  these  developments. 

Presently  there  exist  several  approaches  to  the  formulation  of  micro¬ 
mechanics.  Some  of  these  theories  are  very  general  in  nature  but  incomplete 
and  not  closed.  Others  are  concerned  either  with  special  types  of  material 
structure  and/or  deformations.  Fundamental  ideas  contained  in  some  of  these 
theories  can  be  traced  all  the  way  to  Bernoulli  and  Euler  in  connection  with 
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their  work  on  beam  theories.  In  the  elementary  beam  theory,  with  each  section 
of  the  bar  there  is  associated  two  sets  of  kinematical  quantities,  namely, 
a  deformation  vector  and  a  rotation  vector  and  two  types  of  internal  loads, 
namely,  the  tractions  and  couples.  In  plate  theory,  we  have  a  similar 
situation.  Bar  and  plate  theories  involving  these  independent  quantities 
were  recorded  by  Kelvin  and  Tait[1879],  The  existence  of  the  stress  couples 
independent  of  tractions  is  essential  to  these  theories.  For  three-dimensional 
bodies,  this  concept  is  found  in  the  work  of  MacCullagh  [1839]  in  connection 
with  his  work  on  optics.  Lord  Kelvin  went  as  far  as  building  a  model  of  what 
he  called  "quasi-rigid"  ether  which  is  supposed  to  provide  a  mechanical  model 
for  Maxwell's  theory  of  electromagnetism.  The  existence  and  basis  of  couple 
stress  in  elasticity  was  also  noted  by  Voigt  [1887]  in  connection  with  polar 
molecules . 

In  a  remarkable  monograph,  E.  and  F.  Cosserat  [1909]  gave  a  unified 
theory  for  the  deformable  bars,  surfaces,  and  bodies.  A  Cosserat  continuum 

is  defined  as  a  three-dimensional  continuum,  each  point  of  which  is  supplied 

1  2 
with  a  triad  .  By  the  use  of  a  principle  which  they  call  "euclidean  action" 

and  by  calculating  the  variation  of  the  internal  energy  density,  they  gave 

the  equations  of  local  balance  of  momenta  for  stress  and  couple  stress  and 

the  expressions  of  surface  tractions  and  couples.  In  the  work  of  the  Cosserat 

brothers,  we  find  that  the  effect  of  couple  stress  on  the  motion  of  deformable 

bodies  is  fully  taken  into  account. 

Some  fifty  years  elapsed  after  the  work  of  the  Cosserats  with  very  little 
activity  in  this  field.  The  idea  of  a  Cosserat  continuum  was  revived  in  various 

In  the  terminology  of  Truesdell  and  Noll  [1960,  Arc  256]  "directors",  the 
same  terminology  and  similar  ideas  were  used  by  Toupin  [1964],  Green,  Rivlin 
and  Naghdi  [1965]. 

2 

Equivalent  to  the  principle  of  objectivity,  ci .  Eringen  [1962,  Art.  27]. 
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special  forms  by  GUnther  [1958],  Grioll  [1960],  Aero  and  Kuvshinskii  [1961], 
and  Schafer  [1962]  of  whom  GUnther  also  remarked  on  the  connection  to  the 
theory  of  dislocations.  The  question  of  couple  stress  was  reopened  with 
an  incomplete  theory  of  Cosserat  bars  and  surfaces  included  in  Truesdell 
and  Toupin  [I960].  Mindlin  and  Tiers  ten  [1962],  Toupin  [1962],  Eringen 
[1962,  Arts.  32,  40]  recapitulated  a  special  Cosserat  continuum  now  known 
as  the  (indeterminate)  couple  stress  theory.  In  these  theories,  the  rotation 
vector  is  not  an  independent  vector,  consequently  the  antisymmetric  part  of 
the  stress  and  symmetric  part  of  the  couple  stress  remains  undetermined 
(cf.  Art.  23  below).  Eringen  and  Suhubi  [1964a, b]  and  Eringen  [1964  ]  intro¬ 
duced  a  general  theory  of  a  nonlinear  microelastic  continuum  in  which  the 
balance  laws  of  continuum  mechanics  are  supplemented  with  additional  ones, 
and  the  intrinsic  motions  of  the  microelement  contained  in  macrovolume  AV 
are  taken  into  account.  This  theory,  in  special  cases,  contains  the  Cosserat 
continuum  and  the  indeterminate  couple  stress  theory.  Independently,  a  micro¬ 
structure  theory  of  elasticity  was  published  by  Mindlin  [1964]  and  a  multipolar 
continuum  theory  by  Green  and  Rivlin  [1964],  Both  of  these  theories  appear 
to  have  contacts  with  these  of  Eringen  and  Suhubi,  in  special  situations. 
Following  these  works,  an  intense  activity  began  and  literature  now  contains 
several  hundred  papers  in  this  and  in  related  fields  A  proper  assessment 
of  these  works  with  appropriate  references  is  beyond  the  scope  of  this  article. 

This  article  is  concerned,  basicially,  with  special  types  of  continua 
called  micropolar  continua .  The  theory  was  initiated  by  Eringen  and  Suhubi 
[1964b,  Art.  6]  as  a  special  case  of  their  work  rp  the  microelastic  solid 
and  was  named  couple  stress  theory.  Later,  Eringen  [1966a, b]  recapitulated 
and  renamed  it  micropolar  theory  and  proved  several  uniqueness  theorems.  A 
similar  theory  appears  to  be  given,  independently,  by  Palmov  [1964]  for  the 
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linear  elastic  solid.  While  the  theory  is  fresh  and  no  experimental  work 
has  been  published  as  yet,  we  believe  that  the  results  obtained  so  far  are 
sufficient  to  strengthen  the  future  of  the  theory. 

The  theory  of  micropolar  elasticity  is  concerned  with  material  media 
whose  constituents  are  dumbbell  molecules.  These  elements  are  allowed  to 
rotate  independently  without  stretch  The  theory  is  expected  to  find  appli¬ 
cations  in  the  treatment  of  mechanics  of  granular  materials  with  elongated 
rigid  grains  and  composite  fiberous  materials. 

The  first  eight  sections  (Arts.  2  -  9)  of  this  article  give  a  treat¬ 
ment  of  the  geometry  of  deformation  and  microdeformation,  strain  and  rotation 
measures,  compatibility  conditions,  and  some  illustrative  examples  of  deformation. 
Sections  10  -  13  are  devoted  to  kinematics  and  rate  measures.  External  and 
internal  loads  and  the  balance  laws  are  discussed  in  Sections  14  to  17,  and 
energy  and  entropy  in  Sections  18  and  19.  The  constitutive  equations  of  the 
theory  of  micropolar  elasticity  and  restrictions  on  the  coefficients  are 
derived  in  Sections  20  and  21.  The  field  equations,  boundary  and  initial 
conditions  are  prescribed  and  discussed  in  Section  22.  In  Section  23,  we 
show  how  the  indeterminate  couple  stress  theory  arises  as  a  special  case  of 
the  micropolar  theory.  Sections  24  through  29  are  devoted  to  the  solutions 
of  various  problems. 

Micropolar  continuum  mechanics  is  In  the  stage  of  its  infancy.  The 
linear  theory  is  reasonably  simple  and  it  lends  itself  to  the  solution  of 
some  important  boundary  and  initial  value  problems.  A  large  class  of  unsolved 
problems  and  experimental  work  offer  a  challenge  to  future  workers. 
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2.  DEFORMATION  AND  MICRODEFORMATION 


A  material  point  P  of  a  body  B  having  volume  V  and  surface  S  in  its 
undeformed  and  unstressed  state  may  be  located  by  its  rectangular  coordinates 
X^,  X2,and  X^  (or  simply  X^,,  K  ■  1,2,3),  Fig.  2.1.  If  the  body  is  allowed 
to  move  and  deform  under  some  external  loads,  it  will  occupy  a  region 
having  volume  \J  and  having  surface  S.  Referred  to  the  same  rectangular  frame 
of  reference,  the  new  position  of  the  point  P  will  be  x^,  x2,  (or  simply 
x^,  k  -  1,2,3).  Under  the  assumption  of  indestructibility  and  impenetrability 
of  matter,  each  material  point  in  the  undeformed  body  B  will  occupy  a  unique 
position  in  the  deformed  body  8. 

Conversely,  each  point  in  t'  can  be  traced  back  to  a  unique  point  in  B. 
Thus,  the  deformation  of  the  body  at  time  t  may  be  prescribed  by  a  one-to-one 
mapping 

(2.1)  Xj^  -  xk(X1,X2,X3,t),  k  -  1,2,3 


or  its  inverse  motion 

(2.2)  -  XK(x1,x2,x3,t),  K  -  1,2,3 

We  assume  that  (2.2)  is  a  unique  inverse  of  (2.1)  for  all  points  contained  in 
the  body  except  possibly  some  singular  surfaces,  lines, and  points.  For  this 
to  be  valid,  the  three  functions  x^(X^ ,X2,X3 , t)  must  possess  continuous  partial 
derivatives  with  respect  to  X^,  X2,and  X3  for  all  times,  and  the  Jacobian 


J  :  det 


3Xk 


3xl 

'  3X1 

3xi 

3X1 

9X2 

3X3 

3x2 

3x2 

3x2 

axj 

3X3 

3x3 

3X1 

ax2 

(2.3) 
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must  not  vanish.  Henceforth,  we  assume  this  to  be  the  case. 
The  partial  derivatives 


(2.4) 


„ 

*k,K  "  3XR 


are  called  deformation  gradlentst  and  they  are  basic  in  the  study  of  continuum 
mechanics . 

i 

We  now  consider  a  volume  element  AV  enclosed  within  its  surface  AS  in 
the  undeformed  body.  Let  the  center  of  mass  of  AV  have  the  position  vector  $. 

All  materials  possess  certain  granular  and  fibrous  structure  with  different 
sizes  and  shapes.  If  the  physical  phenomena  under  study  has  a  certain  charac¬ 
teristic  length  (such  as  wave  length)  that  compares  with  the  size  of  grains 
in  the  body,  then  the  microstructure  of  the  material  must  be  taken  into  consid¬ 
eration.  In  such  situations,  classical  continuum  mechanics  must  be  modified 
by  considering  the  effect  of  the  granular  character  of  tne  medium.  Suppose 
that  the  eJement  AV  +  AS  contains  N  discrete  micromaterial  elements  AV^  +  AS^a\ 
(a  ■  1,2,...N),  each  with  a  mass  density  The  position  vector  of  a 

material  point  in  a  microelement  may  be  expressed  as 


(2.5) 


X(“>  -  X  +  E(a) 


(a) 

where  z  is  the  position  of  a  point  in  the  microelement  relative  to  the 
center  of  mass  of  AV  +  AS,  Fig.  2.2.  Upon  the  deformation  of  the  body, 

AV  +  AS  goes  into  Av  +  As  with  the  microelement  displaced  with  respect  to 
the  mass  center.  Because  of  the  rearrangements  and  the  relative  deformations 
of  the  microelements,  the  center  of  mass  P  may  now  move  to  a  position  p  and 
the  material  point  Q  to  a  new  position  q  in  the  deformed  body.  The  final 
position  of  the  particle  will  therefore  be 


x 


(a) 


x  + 


(a) 


(2.6) 
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•  (a) 

where  x  is  the  new  position  vector  of  the  center  of  mass  of  Av  and  §  is 

the  new  relative  position  vector  of  the  point  originally  located  at  . 

The  motion  of  the  center  of  mass  P  of  AV  is  expressed  as  usual  by  (2.1)  or 

simply 


(2.7)  x  -  x(X, t) 

The  relative  position  vector  1  however  depends  not  only  on  X  and  t  but 
also  on  E^  ,  i.e.  , 


(2.8) 


(a) 


,t) 


A  microstructure  theory  must  lean  heavily  on  the  assumption  charac¬ 
terizing  the  nature  of  the  functional  character  of  ^a\  Eringen  and  Suhubi 
[1964a, b]  and  Eringen  [1964]  have  constructed  a  general  tneory  in  which  (2.8) 
is  linear  in  E  .  The  theory  so  constructed  was  later  called  by  Eringen 
the  theory  of  micromorphic  materials .  The  basic  assumption  underlying  this 
theory  is  the 

Axiom  of  Affine  Motion.  The  material  points  in  AV  +  aS  undergo 
homogeneous  deformation  about  the  center  of  mass .  thus  > 


(a) 


xx(x,t) 


(a) 
'  1 


+  x2(x,t)  =  (2a) 


X3(X,t)E(3a) 


a  *  1,2, .  .  .N 


or  simply 

(2.9)  £(a)  =  XK(X,t)5iJa) 


where  summation  over  repeated  indices  is  understood  This  assumption  is  justi¬ 
fiable  on  the  physical  grounds  that  when  AV  +  AS  is  small  enough,  its  motion 
consis ts  of  a^  trans lation .  a_  rotation  about  i ts  center  of  mass .  and  a^  homogeneous 
deformation .  Note  that  in  all  classical  theories  of  continuous  media,  this  last 


assumption  is  missing  (cf.  Eringen  [1962,  Art.  10]). 
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Clearly,  the  theory  can  be  extended  to  include  quadratic  and  higher  order 
terms  in  5^  .  However,  it  quickly  takes  a  very  complex  form,  losing  its 
simplicity  and  usefulness.  In  this  article,  in  fact,  we  shall  be  concerned 
mostly  with  a  much  simpler  case,  namely,  the  theory  of  micropolar  elasticity^. 

This  latter  theory  admits  only  rigid  microrotations  of  the  micro¬ 
volume  elements  about  the  center  of  mass  of  the  volume  element.  In  other 
words,  the  basic  assumption  (2-9)  is  further  simplified  by  placing  further 

restrictions  on  the  three  vector  functions  x.,*  As  we  shall  see  later,  this 

■*  K 

will  amount  to  reducing  the  number  of  the  microdeformation  functions  xK  from 
three  to  one.  In  fact,  we  shall  also  be  dealing  mainly  with  the  linear  theory 
In  classical  continuum  mechanics,  the  problem  is  centered  around  the  determina 
tion  of  the  spatial  position  x  of  all  material  points  of  the  body  at  a  given 
instant.  This  means  that  when  we  are  through  with  all  calculations,  we  will 
have  the  three  functions  x^(X,t)  determined.  In  the  theory  of  micromorphic 
materials,  in  addition  the  three  vector  functions  xv(x»t)  (equivalently 
nine  scalar  functions)  must  also  be  determined.  The  complicated  nature  of 
the  problem  and  the  necessity  of  additional  physical  concepts  and  laws  are 
now  apparent. 

In  coordinate  form  for  the  spatial  position  of  a  material  point  in  a 
microelement,  we  have 


(2.10) 


r 


xk(X,t)  +  xkK(X,t)E<Q) 


k,K  =  1,2,3 


It  is  now  clear  that  we  have  3  +  9  *  12  functions,  x^(X,t)  and  xkK(X»t),  to 
determine  the  spatial  position  of  the  material  point 

Just  as  in  (2.2),  we  introduce  the  inverse  micromotions  such  that 


(2.11) 


XkKXKJ! 


■'k£* 


"kK^Lk 


KL 


1 


The  theory  was  developed  in  the  original  paper  of  Eringen  and  Suliubi  [1964b, 
Sec.  6]  as  a  special  case  of  the  general  theory  and  was  called  the  couple 
stress  theory.  Eringen  [1965]  later  named  it  the  micropolar  theory. 
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Here  and  throughout  this  article,  summation  is  understood  to  be  over  the 
repeated  indices,  e.g. , 

xkn\l  XklXli  +  Xk2X2£  +  Xk3X3£ 

The  symbols  5  and  6  are  the  Kronecker  deltas  which  are  1  when 

K  X>  rvL 

the  indices  take  the  same  numerical  value  and  zero  otherwise. 

In  component  form  (2.9)  reads 

(2.12)  Ck“}  “  XkK(^t)HKa) 


Upon  multiplication  of  both  sides  by  and  using  (2.11),  we  also  get 
(2.13)  E<“>  -  \k(x.t)E<0) 


In  vector  form  this  reads 


(2.14) 


-(a) 


4<x.tKk 


<t>) 


The  rectangular  components  of  are  denoted  by  X^  and  those  of  xk 
by  XkR,  i.e.  , 

(2. 1.5)  XR  -  XkK(X,t)ik  ,  xk  -  xRk(x,t)IK 

where  I  and  i,  are,  respectively,  the  unit  base  vectors  for  the  material 
coordinates  X^,  and  the  spatial  coordinates  x^ 

The  motion  and  the  inverse  motion  of  a  material  point  in  a  micro¬ 
element  are  therefore  expressed  by 

(2.16)  x<°°  -  x^X.t)  +  xkK(X,t)^a) 

<a)  ■  V*'0  +  V!'dl“ 


(2.17) 
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In  vectorial  form  these  read 

(2.18)  x(a)  -  x(X,t)  +  xK(X,t)=^a) 

(2.19)  X(o)  -  X(x,t)  +  Xk(x,t)^“) 

It  is  clear  that  we  may  employ  either  representations  (2.16)  or 
(2.17),  and  that  whenever  either  set  of  functions  (xjc»X]cj^)  or 
are  determined,  the  problem  is  completed  since  the  other  set  is  soluble 
from  the  one  that  is  found. 
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3.  STRAIN  AND  MICROSTRAIN  TENSORS 


The  differential  line  element  in  the  deformed  body  is  calculated 
through  (2.18) 

O.i)  dx  -  (x,R  +  XL>K=L  +  xRd-R 


where  an  index  followed  by  a  comma  denotes  partial  differentiation.  This 
convention  will  be  used  throughout  this  article.  Thus,  for  example, 

3x  _  dXj 

(3,2)  -  *K  =  »  XL,K  =  3X^ 


(3.3) 


3xk 


are  used  for  brevity.  Note  that  (3.2)^  and  (3.3)^  are  the  classical  defor¬ 
mation  gradients,  and  (3.2)2  and  (3.3)2  are  the  microdeformation  gradients 
of  the  present  theory. 

The  square  of  the  arc  length  is  now  calculated  by  forming 


dx 


(a) 


dx 


(a) 


Upon  using  (3.1)  and  forming  the  scalar  product,  we  find 


^  <dS<“,)2  '  (CKL  +  2  rKKL=M  +  *kH,K*kN,LSM;»>  dV\ 


+  2(yKL  +  XkLXkM,K“M)  dXKd"L 


+  XkKXkLd"Kd"L 


where  we  also  dropped  the  superscript  a  on  E  and  dE  for  brevity,  since 

K  K, 

this  is  understood  whenever  we  use  the  letter  E(and  also  O.  In  (3.4) 
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we  introduce  the  notations 


(3.5) 

cKL(X,t)  = 

(3.6) 

YKL(X,t)  =  \(KXkL 

(3.7) 

Of  these,  is  the  classical  Green  deformation  tensor  and  4^ 

and  f .  are  new  microdeformation  tensors  of  the  present  theory. 

KLM 

/  \ 

We  now  introduce  the  displacement  vector  u  as  the  vector  that 
extends  from  to  x^a\  Fig.  (3.1).  Thus  we  write 

(3.8)  u(a)  -x-X+C-==u+§-E 

where 


(3.9) 


u  =  x  -  X 


is  the  classical  displacement  vector, 
are,  respectively. 


tne  components  or  wnicn  in 


"K 


Tc 


(3.10) 


K 


=  u 


'K 


’SAk  "  \ 


(3.11) 


\  5  U  ’  *k  "  *k  '  VKk 


where 


(3.12) 


Since  the  spatial  and  material  frames  are  taken  to  be  the  same  rectangular 
frame  of  reference,  <5  is  none  other  than  the  Kronecker  'elta  which  has 
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the  value  one  when  the  two  indices  take  the  same  numerical  value  and  zero 

i  x^  and 

the  convention  of  majuscule  indices  for  the  material  frame  and  miniscule 
indices  for  the  spatial  frame  of  reference.  This  convention  is  espec  xv 
useful  in  the  nonlinear  theory. 

From  (3.10)  and  (3.11),  by  partial  differentiation,  we  obtain 


(3.13) 

\,K  ‘  <6LK  +  UL,K)SkL 

(3.14) 

^K.k  ■  <Stk  '  ut,k*6Kl 

otherwise. 


It  is  possible  to  write  x^6 


kK 


VlCk  E  V  But  We  keep 


Similarly,  we  introduce  the  microdisplacement  tensors  4>  (X,t)  and  .  (x,t) 

by 


(3.15) 


'kK 


(6LK  + 


\K)6kL 


(3.16) 


\k  ‘  (6Jlk  "  *£k*6K£ 


By  use  of  (3.9)  and  (3.15)  and  (3.16),  we  see  that  (3.8)  may  also  be  expressed 
as 


(3.17) 


(a) 


'  +  5  *  =  (UK  +  ^KITL^'K 


(3.18) 


(a) 


u  +  =  (uR  +  ^^)ik 


Upon  substituting  (3.13)  and  (3.15)  into  (3.5)  to  (3.7),  we  find  that 


(3.19) 

CKL  “ 

^KL  +  UK,L 

(3.20) 

\l  “ 

6KL  +  *KL  ' 

(3.21) 

rKLM  " 

^KL.M  +  UN 
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So  far  all  these  expressions  are  exact.  For  a  linear  theory,  one  assumes 
that  the  product  terms  are  negligible  so  that 


(3.22) 


CKL  6KL  +  UK,L  +  UL,K 


(3.23) 


V  -6  +  4>  +  U 

KL  KL  KL  L,K 


(3.24) 


r  4- 

KLM  KL,M 


In  this  case,  the  difference  between  the  spatial  and  material  representations 
disappears  so  that  one  may  use  in  place  of  U^,  and  4^  in  place  of  4^,  etc., 
a  fact  which  is  well  known  in  the  classical  continuum  theory  (cf.  Eritigen 
[1962,  Art.  14]).  For  the  microdeformation,  this  may  be  seen  as  follows:  If 
we  use  (3.15)  and  (3.16)  in  (2.11),  we  obtain 


(3.25) 


4>  x  (&  +  4>  )  4>  6  6 

KL  v  KM  KM/vm£MmzL 


Neglecting  the  product  terms,  we  see  that 


(3.26) 


4  -  0  <5  i 

KL  m*  Kir.  IL 


which  is  the  proof  of  our  statement. 

Since  we  will  be  dealing  with  the  linear  theory,  we  shall  not  distin¬ 
guish  between  material  and  spatial  representations  except  when  it  becomes 
necessary  for  clarity  in  the  development 

The  material  (or  lagranRian)  strain  tensor  E„.  and  the  material  nicro- 
strain  tensors  E  and  f  are  defined  in  the  linear  theory  by 


(3.27) 


EKL  2(CKL  "  "kL)  =  2(UK,L  +  UL,K) 


(3.28) 


E  if  -  5  =4  +U 

KL  KL  KL  KL  L,K 


(3.29) 


r  =4' 

KLM  KL ,M 
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In  the  light  of  what  has  been  shown  above,  we  may  also  introduce 
the  spatial  (or  eulerlan)  strain  tensor  e^  and  spatial  mlcrostraln 
tensors .  anc*  in  a  similar  fashion: 


(3.30) 


ekl  2  Cuk,t  +  “l.k' 


(3.31) 

(3.32) 


Set  =  Set  +  U£,k 


Y  *  —  (h 

ktm  k£,m 


Clearly,  when  these  tensors  are  known,  changes  in  arc  length  and  angles 
during  any  deformation  of  the  body  can  be  calculated. 

For  the  linear  theory,  the  difference  between  the  squares  of  arc 
length  in  the  deformed  and  undeformed  body  follows  from  (3.4)  and  using 
(3.15)  and  (3.27)  to  (3.29) 

2  2 

<dS<a>)  -  Ms(a))  »  2(ekl  + 


(3.33) 


+  2(EKL  +  rLMK"M)dXKd"L 


+  (tKL  +  ELK  '  2EKL)d::Kd;:L 


In  classical  continuum  mechanics,  only  the  first  term  on  the  right  involv 

ing  E  survives. 

K.L 

From  (3.33)  it  is  clear  that  when  E  ,  E  and  r  vanish,  there 

K.L  KL  KLM 

will  be  no  change  in  the  arc  length  after  a  deformation.  In  such  a  situa 
tion,  the  body  is  said  to  undergo  a  rigid  motion, 
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4.  MICR0P0LAR  STRAINS  AND  ROTATIONS 


We  now  consider  a  special  class  of  materials  in  which  the  state 
of  the  microdeformation  can  be  described  by  a  local  rigid  motion  of  the 
microelements.  A  large  class  of  materials  exist  in  which  the  micro¬ 
material  elements  are  dumbbell  type  molecules.  Materials  consisting  of 
rigid  fibers  or  elongated  grains  fall  into  this  category  For  example, 
wood,  certain  rocks,  and  minerals  contain  elongated  molecular  elements. 
Among  fluids,  blood  possesses  dumbbell-like  molecules.  For  such  media, 
the  micromorphic  material  theory  becomes  much  simpler.  Mathematically, 
this  specialization  is  obtained  by  setting 


(4.1) 


4> 


LK 


or  in  the  spatial  notation  <J> 
skew-symmetric  second-order 
<t>K  defined  by 

(4.2)  4- 


o  —  d) 

kH. 

tensor  $ 


ik' 

KL 


In  three-dimensional  space 
can  be  expressed  by  an  axial 


G 

KLM  ML 


,  every 
vector 


where  e  is  the  alternating  tensor  defined  as 
K.LM 

G  m  €  m  £  “  —  r  =  —  r  =*  —  r 

123  231  312  213  132  321 

c,„ *  0  otherwise 
KLM 

Expression  (4.2)  is  a  compact  expression  of 


A  ^  s  $  d>  ««|) 

1  32’  2  13*  3  21 


The  solution  of  (4.2)  for  <t>  is  given  by 

K.L 


KL 


(4.3) 


r  'f' 

KLM  M 
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Subs  tituting  this  into  (3.15),  we  sea  that 


(4.4) 


y  “6  —  e  0 

xkK  kK  kKM  M 


In  the  classical  theory  we  have  the  rotation  tensor 


(4.5) 


.  1 


^CL  "  -RLK  =  2  (UK,L  ”  UL,K) 


The  axial  vector  corresponding  to  this  is  given  by 


(4.6) 


\  "  2  £'KLMRML  ”  2  EKLMUM,L’ 


EKLMRM 


Using  (3.27)  and  (4.6)2,  we  find  that 


(4.7) 


U”  ’  "  E'”  +  *SCL  “  EKL  “  EKLMRM 


K,L  KL 


When  this  and  (4.2)  are  substituted  into  (3.28)  and  (3.29)^,  we  get 


(4.8) 


ekl  •  ekl  +  cklm(rm  -  V 


(4.9) 


KLM 


_EKLM<fN,M 


When 


and  the 


*M  *  V  “ee  See  that  EKL  '  EKL  an<l  rKLM  ’  V,M 

microstrains  are  no  longer  independent  of  the  classical  strains  and 
rotations.  Thus  the  micropolar  theory  assumes  that  the  classical  rota¬ 
tion  Rj,  is_  different  from  the  microrotation .  In  the  micropolar  theory, 
we  have  therefore  six  functions  to  determine,  namely  U  (X,t)  and  4>  (X,t) 
Once  this  is  done,  length  and  angle  changes  can  be  fully  calculated. 

For  the  micropolar  theory,  the  spatial  position  of  the  a*’*1  point 
x(a)  is  obtained  through  (3.8),  (3.18),  and  (4.3),  i.e., 
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(4.10)  x(a)  -X+E+u-E*$ 

From  this  it  is  apparent  that  $  represents  an  angular  rotation  of  a 
microelement  about  the  center  of  mass  of  the  deformed  macrovolume  element 
(£  is  the  moment  arm  from  this  centroid),  Fig.  4.1.  Accordingly,  we  also 
have 


(4.11) 

which  shows  that  aside  from  a  rigid  body  translation  the  relative  position 
of  E  of  a  material  point  after  deformation  is  obtained  by  translating  5 
parallel  to  itoelf  to  the  center  of  mass  x  of  the  deformed  macrovolume 
element,  and  then  rotating  it  an  amount  E  *  $.  Since  we  are  dealing  with 
linear  theory,  we  also  have 


(4.12) 


S  +  £  *  0 


where  #  -  is  the  spatial  microrotation.  In  fact,  we  have  the  complete  dual 
to  equations  (4.1)  to  (4.10)  which  we  record  here  for  future  convenience. 


(4.13) 


^k  2  ''kim'^mi  ’  ^ki  Ckim^m 


(4.14) 


^Ki  Ki  'Kim^m 


(4.15) 


r  ■  —  r  r  .  r  *  —  r  r 
k  2  kirn  mi  ’  ki  kirn  m 


2  kimUm,i 


(4.16) 


\,i  eki  "  ekimr 


m 


(4.17) 


'  ek«  +  £kt»(rm  '  *m> 


(4.18) 


Y  *  C  <t) 
kim  kin  n,m 


-23- 


and 

(A. 19)  x(a)  =  X+  5  +  u-  5x$ 

Now  consider  the  deformation  of  an  infinitesimal  vector 

( 

dx'  ■  dX  +  dE  at  .  4 E.  Upon  deformation,  this  vector  becomes 
(A. 20)  dx(a)  -  dx  +  d£  -  dX  +  dE  +  u,  dX^  -  d:  *  t  -  :  *  *,  dX^ 
By  use  of  (A. 7)  and  (A. 6)  we  may  write 


'*KdXK  ”  UL,KdXK^L  "  ELKdXK^L  + 


=  Wl  “  d?  *  * 


Similarly  using  (A. 9) 


' ’NdXN  “  EKLM"L$M,NdXN^K  ”  rKLN"LdXN~K 


For  convenience  we  now  introduce  the  notation 


(A. 22) 


r 


KLM“L 


so  that 

(4,23)  =  *  *»NdXN  “  “r:KM)dXMIK  _  r[KM)dXMIK 

where  indices  in  parenthesis  (and  brackets)  indicate  the  symmetric  (and 
antisymmetric)  parts  .  Carrying  (A. 21)  and  (A. 23)  into  (A. 20)  we  rearrange  it  into 

(A.2A)  dx(a)  =  dX  +  dE  -  (dX  »  R  +  d:  M  +  dX  <  [) 

+  (EKL  +  r(KL))dXK'L 


where  we  also  defined  a  new  microrotation  vector  f  by 
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(4.25) 


k'KLMrML 


r[MK]  _CKLMrM 


We  name  this  vector  minirotation  for  distinction  from  the  microrotation 

<t>.  If  we  carry  (4.9)  into  (4.25),,  we  also  find 
-  1 


(4.26) 


2  K  L,L  K  K,L  1/ 


Equation  (4.24)  reveals  that  the  deformation  of  the  vector 

(a) 

dX  =  dX  +■  d—  may  be  achieved  by  the  following  three  operations: 

(a)  A  rigid  translation  of  dX  +  dS  from  the  material  centroid 
X  to  the  spatial  centroid  x. 

(b)  Rigid  rotations  of  dX  and  dE  by  the  amounts  dX  *  (R  +  T)  and 
dE  *  $,  respectively. 

(c)  Finally,  a  stretch  represented  by  the  strains  and 

in  (4.24). 

The  following  special  cases  help  to  visualize  these  deformations: 

(i)  When  :  =  0  we  have 

(4.27)  d$(a)  *  dx  =  dX  -  dX  *  R  +  E^dX^ 


This,  of  course,  is  a  well-known  theorem  in  classical  continuum  mechanics 
attributed  to  Helmholtz  (cf.  Eringen  [1962,  Art.  10]). 

(ii)  When  E  «  constant  vector 

(4.28)  dx(Q)  *  dX  -  dX  *  (R  +  [)  +  (E^  +  r(KL))dXK?L 

Here,  oi  course,  we  have  no  rotation  of  d:  0. 

(iii)  <t>  =  constant  microrotation.  In  this  case  we  have 

r(KL)  "  r[KL]  “  rK  =  °  and  We  °btain 
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(4.29) 


dx 


(a) 


dX  +  d:  -  dX  i  R  +  d:  *  $  + 


EKLdV-L 


In  this  case,  the  rotation  consists  only  of  a  macrorotation  of  dX  and  a 
microrotation  of  d£.  We  can  also  write  it  as 


(4.30) 


d£  ■  dx 


(a) 


dx  =  d:  -  dX  «  r  -  d5  *  R  +  r,„T  ,  dXL,IT 

-  \  KL )  K.  L 


.(a) 


In  this  form  we  see  the  difference  between  the  deformation  of  dX  and 
that  of  dX  of  which  the  latter  is  known  to  us  from  the  classical  theory. 
This  difference,  therefore,  is  the  result  of  the  composition  of  a  mini¬ 
rotation  of  dX,  a  macrorotation  of  dE  ,  and  the  ministraining  of  dX 

characterized  by  r.  ..  The  terminology  of  minirotation  is  being  used 
(.K.L; 

for  T  and  ministraining  for  r.^).  Of  course,  R  is  the  classical  rotation 
for  which  we  use  the  terminology  macrorotation. 

For  the  spatial  representation  the  dual  of  (4.10)  is 


(4.31) 


x+C-u+S* 


From  this,  in  the  same  way  as  in  the  case  of  (4.24),  we  obtain 
(4.32)  dX^1^  *  dx  +  d£  +  (dx  "  r)  +  d£  *  $  +  dx  *  y) 

-  (eki  '  v<ki))dxk1t 


where  r  is  the  spatial  macrorotation  vector  defined  by  (4.15)  and  e^  is 
the  spatial  macrostrain  tensor.  The  spatial  minirotation  vector  y  is  given 
by 


(4.33) 

and 


2  'kim^mJl 


2  ^1,1^ k  ^k.JlV 


(4.34) 


km 


^kim^it 


Lk£nvn,m  l 
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Finally,  below  we  record  the  component  form  of  various  strains  and 
rotations  in  spatial  rectangular  coordinates  x,  y,  z.  For  the  components 
of  the  displacement  vector  u  and  micropolar  rotation  vector  $  we  write, 
respectively,  u,  v,  w  and  $x»  $  ,  tp^ . 

Macrostrain  tensor: 


3u 

>  m  —  0 

'xx  3x  ’  xy 


A.  cAil.  +  3v  . 
2  {dy  3x 


(4.35) 


3v 

'yy  “  3y  »  6yz 


1  ,_3v  3w, 

2  C3z  3x; 


3w  _  1_  ._3w  _3u. 

ezz  3z  ’  ezx  2  ^3x  3z' 


Micropolar  strain  tensor: 


(4.36) 


3u  3v  3w 

Cxx  3x  *  Cyy  3y  *  Ezz  3z 


3v  3u  . 

c  =  - - <p  ,  c  «  $ 

xy  3x  z  yx  3y  Tz 


3w  3v  , 

c  =  —  ~  .  e  =  - —  +  it 

yz  3y  x  zy  3z  x 


3u  3w  , 

e  *  - - ,  e  =*  —  +  $ 

zx  3z  y  xz  3x  y 
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Mlcropolar  strain  tensor  of  third  order: 


3$ 


3 


34 


yzx 


zyx  3x  ’  'yzy  ' zyy  3y  ’  'yzz 


-Y 


zyz 


x 

3z 


(A. 37)  Y  -  -Y 


zxx 


34> 

r  -  — 1  Y 

xzx  3x  ’  'zxy 


3  $ 

V 

-Y  -  -r-  .  Y 

xzy  3y  zxz 


-Y 


xzz 


34> 

3z 


3  4 


34 


3  4 


y  m  — y  m  Y  b  — y  =  —  *  Y  «  —  Y  «  "  " “■ 

xyz  yxx  3x  *  xyy  yxy  3y  ’  xyz  yxz  3z 


all  other  y,  „  *  0. 

'k£m 


Macrorotation  vector: 


_1  ,  3w  3v. 

rx  2  ^3y  3z 


(4.38) 


JL  /JLy. 

2  Sz 


3w. 

3x; 


1  .  3v 

2  ^3x 


3u. 

3y; 


Microrotation  vector: 


(4.39) 


i  +  4  i  +  4  k 

V*  Z- 


Minirotation  vector: 


34  34  34  34 

v  -  i  [  ( — 3L  +  - *r - ] 

Yx  2  3y  3z^x  3y>  3z  zJ 


34  34  34  34 

v  -  -  r  ( — -  +  — £.)  r _ If _ If  ] 

y  2  3x  3z  y  3x^x  3z  z 


34  34  34  34 

1  ,  , _ x  ,  __y .  _  z.  _ z _  , 

YZ  ■  2  ((—  +  -S?5*  -  — '  ~l ’y] 


(4.40) 
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5.  GEOMETRICAL  MEANING  OF  MICROPOLAR  STRAINS  AND  ROTATIONS 

The  geometrical  significance  of  various  strain  and  rotation  measures 
can  be  understood  quickly  if  we  consider  the  change  dx^  in  x^  as  the 
vector  sum  of  three  changes,  namely, 


(5.1) 

dx^  ■  dx  +  dy  +  dx 

where 

(5.2) 

dx  ■  dX  -  dX  *  R  +  E^dX^I. 

(5.3) 

d*=-dX*E+  F (KL) dXK'L 

(5.4) 

dz  =  d;  -  d;  x  4 

Here  dx  is  that  known  to  us  from  classical  continuum  theory  (cf.  Eringen 
[1962,  Art.  6]).  Accordingly,  on  the  right-hand  side  of  (5.2)  the  first 
term  dX  represents  translation  of  dx  from  X  to  x;  the  second  is  a  rotation; 
and  the  last  term  represents  the  straining  of  the  body.  More  specifically, 
consider  the  vector  dX  at  the  point  X  of  the  undeformed  volume  element  dV. 
This  vector  after  deformation  becomes  dx. 

Writing  (5.2)  in  another  form 

(5.5)  dx  *  Cj,dX^, 


where 

(5.6) 


CK  DXj,  lK  +  UM,K-M 


we  see  that  a  parallelepiped  with  side  vectors  I^dX^,  I^dX,,,  and  I ;dX^ 
after  deformation  becomes  a  rectilinear  parallelepiped  with  side  vectors 
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C^dX^,  and  C^dX^,  Fig.  5.1.  The  stretch  and  extension  are 

defined  by 


(5.7) 


e(n)  a(n) 


Now  calculate  the  extension  of  one  of  the  sides  of  the  undeformed  parallele¬ 
piped  at  X,  e.g.,  I^dX^: 

|dx| 2  =  C1dX1  •  C1dX1  =  Cn(dX1)2 


Hence 


E 


(1) 


-  1  - 


1 


But  we  have 


1  +  2E 


11 


Thus 


(5.8) 


E(l)  "  \l) 


1  =  /I  +  2E  -  1 


From  this  it  follows  that 


(1  +  E(i)r 


1 


For  small  extensions  E 


(1) 


<<  1. 


and  this  approximates  to 


(5.9) 


-  E 


(1) 


which  provides  a  meaning  for  the  normal  components  E^,  ^22*  ant^ 
infinitesimal  strain  tensor.  For  the  shear  strains  E^*  ^23*  ^1  We 
a  geometrical  meaning  by  calculating  the  change  of  angle  between  two  side 
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vectors  such  as  I^X^  and  The  angle  0^12)  between  these  two  vectors 

after  deformation  is  calculated  by 


cos  0 , 


C1dX1  •  C2dX2 
'  C,  dX,  IcJdX, 


The  change  of  angle  "  J  “  ®(12)  between  tbe  original  and  the  final 

angle  follows  from  this 


(5.10) 


sin  r(12) 


For  small  strains  E^  <<  1,  E22  <<  1  and  sin  ^q2j  ~  ^(^2)  S0  tbat  we  bave 
the  approximation 


(5.11) 


r  =  ?e 
(12)  "l2 


which  provides  the  geometrical  meaning  for  the  shear  strain  E^.  Similar 
results  are,  of  course,  valid  for  E^  and  E^. 

For  the  rotation  vector  R  we  have  the  following  geometrical  inter¬ 
pretation:  Let  be  a  unit  vector  in  the  X^ ,  X2~plane  at  £.  After  defor¬ 
mation,  becomes  a  vector  n^  at  x.  Bring  n;  to  X  and  obtain  its  projection 
n^  on  X^ ,  X2~plane,  Fig.  5.2.  Let  the  angle  between  and  n^  be  denoted  by 

9y  We  can  show  that  (cf.  Eringen  [1962,  Art.  10))  the  average  <  tan  0^> 

over-all  angles  <J>  that  can  make  with  the  X^-axis  is  related  to  R^  by 


(5.12) 


<  tan  0^> 
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where  E^,  E.^,  and  E^  are  the  infinitesimal  strain  components  and  R  is 
the  rotation  component  in  the  X^-direction,  i.e., 

1  3U2  3U1 

D  m  _R  .  —  ( _ —  _  _ L\ 

3  12  2  9X2; 


For  small  strains  (5.12)  can  be  approximated  to  give 


(5.13) 


<  o3>  *  r3 


which  provides  a  geometrical  meaning  for  R3>  Similar  interpretations  are, 
of  course,  valid  for  and  R^ .  Accordingly,  for  small  deformations  and 
rotations,  R  represents  the  average  rotation  angle  of  the  diagonal  vector 
dX  of  the  undeformed  parallelepiped. 

The  above  discussion  concerning  the  deformed  parallelpiped  with 
diagonal  dx  at  x  may  be  used  to  provide  a  geometrical  meaning  for  the 
deformed  parallelepiped  with  diagonal  dy  at  x.  To  this  end,  from  (5.3) 
we  observe  that  T  takes  the  place  of  the  rotation  R  and  r  that  of  E^. 
Accordingly,  the  undeformed  parallelep'ped  with  diagonal  vector  at  X  +  : 
with  S  fixed,  undergoes  an  additional  rotation  T  and  length  and  angle 
changes  represented  by  the  strains  aC  the  spatial  point  x  +  £  that 

the  point  X  +  E  occupies  after  deformation.  This  deformation  takes  place 
with  E  being  displaced  parallel  to  itself.  This  deformation  emanates  pri¬ 
marily  from  the  parallelepiped  with  diagonal  dX  being  constructed  at  the 
point  X  +  E  with  E  fixed  imtead  of  being  at  X,  as  in  the  first  case. 

Finally,  a  parallelepiped  with  diagonal  dE  constructed  at  X  +  E 
with  X  fixed  is  rigidly  displaced  to  a  spatial  point  x  +  £  and  rotated  with 
the  microrotation  $  at  that  point.  Thus  the  final  rectilinear  parallelepiped 
with  diagonal  dx^  at  x  +  (  is  constructed  as  a  result  of  these  translations. 


rotations,  and  strains. 
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The  picture  of  the  deformation  in  these  latter  cases  is  clarified 
further  if  we  write  (cf.  (3.1)) 

(5.14)  dx(a)  -  (CR  +  CR) dXR  +  xKd=R 

where 


(5.15) 


*L,K”L 


and  Xv  related  to  $  by  (3.15). 

In  the  form  (5.14),  the  deformation  of  a  parallelepiped  with  diagonal 
dX  at  X  +  E  and  another  one  with  dE  are  shown  on  Figure  5.3.  According  to 
this  picture,  dX  at  X  +  5  with  fixed  E  becomes  (C^  +  C^JdX^,,  Fig  5.4,  and 
the  one  with  diagonal  dE  at  X  +  E  with  X  fixed  becomes  x^dX^,  Fig.  5.5.  The 
resulting  deformation  of  dX  +  dE  at  X  +  E  is  the  vector  sum  of  these  two 
deformations,  Fig.  5.3. 


-33- 


6.  INVARIANTS  OF  STRAIN  TENSORS 


The  state  of  local  deformation  at  a  point  X  +  E  of  a  micromorphic 
material  is  fully  determined  when  the  three  material  strain  tensors 


(6.1) 


'KL 


£  r 

KL  ’  KLM 


are  given.  With  knowledge  of  these  tensors,  we  can  calculate  length 
(cf.  (3.30))  and  angle  changes  (Art.  5)  and  construct  the  spatial  locations 
of  various  material  points  in  the  body.  Instead  of  the  list  (6.1),  we  may 
of  course  employ  the  spatial  tensors  e^,  and  .  An  important 

question  in  continuum  mechanics  is:  If  the  material  (spatial)  coordinates 
are  rigidly  rotated  at  X  (at  x) ,  are  there  some  functions  of  the  material 
(spatial)  strain  measures  which  remain  unchanged?  The  answer  to  this  ques¬ 
tion  is  provided  by  the  theory  of  invariants.  In  fact,  the  theory  of  invariants 
is  concerned  with  a  more  difficult  question,  namely:  To  determine  the  complete 
set  of  invariants  (called  integrity  basis)  of  a  given  set  of  vectors  and 
tensors  which  are  unchanged  under  an  arbitrary  group  of  transformations  of 
coordinates.  The  minimal  basis  is  a  subset  of  these  invariants  which  can 
be  employed  to  express  all  other  invariants  of  a  given  set.  The  answer  to 
this  question  is  important  on  two  accounts: 

(i)  Constitutive  equations  must  be  form-invariant  under  rigid  motions 
of  the  spatial  frame  of  reference.  This  is  known  as  the  principle  of  objec¬ 
tivity  .  Investigation  of  this  restriction  on  the  constitutive  equations  often 
requires  knowledge  of  the  invariants  of  constitutive  variables  such  as  strain 
measures . 

(ii)  The  Material  symmetry  places  restrictions  on  the  form  of  the  con¬ 
stitutive  equations  when  the  material  frame  of  reference  is  transformed 
according  to  some  group  of  transformations.  For  example,  when  the  material 
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has  a  plane  of  symmetry,  the  constitutive  equations  should  not  change  their 
forms  when  a  reflection  of  axes  is  performed  according  to  the  plane  of 
symmetry.  Similarly,  for  isotropic  materials,  the  constitutive  equations 
remain  form-invariant  under  the  full  group  of  orthogonal  transformations 
of  the  material  frame  of  reference.  For  anisotropic  materials  with  some 
symmetry  axes,  the  group  of  transformations  is  less  restrictive. 

The  invariants  of  a  symmetric  second-order  tensor  (such  as  E^)  in 
three  dimensions  are: 


(6.2) 


XE  E 


JKK 


II. 


5  I  <ekkell 


“  eklelk} 


III  =  det  E 


KL 


Instead  of  this  set,  one 


(6.3) 


may  also  employ 


tr  E  zE 


KL 


tr  ?2  5  eklelk 


tr  E  ^eklelmemk 


The  above  two  sets  are  related  to  each  other.  In  fact 


(6.4) 


tr  E  -  I£  .  tr  E2  -  I2  -  2I1£ 
tr  E3  ■  I2  -  3IeIIe  *  311Ie 


The  determination  of  the  minimal  integrity  basis  of  (6  1)  is  much  more 
complicated.  In  fact,  to  our  knowledge  the  basic  invariants  of  third  and 
higher-order  tensors  have  not  been  studied  to  date  Fortunately,  in  the 
present  theory  T  always  occur  in  the  form  =  H  r  so  that  we  may 
instead  search  for  the  invariants  of 


e  E  r 

KL  ’  KL  ’  KL 


(6.5) 
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For  the  micropolar  theory,  the  situation  is  simplified  further 


since  E^^  »  E^.  The  invariants  of  two  symmetric  and  two  antisymmetric 
second-order  tensors 


(6*6)  L(KL)  “  EKL  *  £[KL]  *  r(KL)  *  r[KL] 

would  be  sufficient  for  this  purpose.  The  integrity  basis  for  the  proper 
orthogonal  group  for  such  a  set  has  been  studied  by  various  authors.  Below 
we  give  a  table  for  the  construction  of  these  invariants  For  the  sake 
of  simplicity,  we  introduce  the  symbols  a,  b  for  the  symmetric  tensors  and 
u,  v  for  the  antisymmetric  tensors.  In  Table  6.1  we  give  the  invariants 
of  these  tensors  in  ascending  order  of  the  integrity  basis  of  various  sub¬ 
sets  of  a,  b,  u  and  v.  The  integrity  basis  of  the  quantities  listed  in  each 
entry  of  the  first  column  includes  all  entries  on  this  row  and  the  integrity 
basis  of  all  subsets  of  these  quantities.  Thus,  for  example,  the  integrity 

basis  for  a,  b  in  the  second  row  includes  that  of  a  and  b,  namely,  tr  a, 

2  3  2  3 

tr  a  ,  tr  a  ,  and  tr  b,  tr  b  ,  tr  b  .  Also  an  asterisk  (*)  placed  on  the 
products  indicates  that  we  include  in  this  list  all  other  products  obtained 
from  this  by  cyclic  permutations  of  the  symmetric  matrices.  A  dagger  (t) 
indicates  the  inclusion  of  all  quantities  obtained  by  cyclic  permutation  of 
the  skew-symmetric  matrices  Thus,  for  example,  ab*  means  the  inclusion  of 
the  set 

ab  ,  ba 

2 

Similarly,  u  va  means  the  inclusion  of  the  set 

2  2 
u  va  ,  vu  a+ 

For  other  details  and  more  extensive  studies  on  invariant  theory  see  Spencer 
[196b j.  Since  we  will  be  interested,  generally,  in  a  linear  theory,  many 
entries  in  this  list  will  not  be  needed  in  the  construction  of  the  constitutive 


equations . 
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7 .  VOLUME  CHANGES 


Here  we  calculate  the  change  of  volume  with  deformation.  A  volume 
element  dV^  i  dX^dX^dX^  at  X  +  5  with  fixed  5  will  be  called  a  material 
macrovolume  element  and  one  dU^  dE^^dE^  with  fixed  X  will  be  called  a 
material  mlnlvolume  element.  After  deformation,  dV^  becomes  dv  and  dUn 
becomes  dv  given 

(7.1)  dv  =  JdX1dX2dX3 


(7.2) 


dv  =  jd-1d=2dE3 


where  J  and  j  are  the  jacobians  of  deformation  with  E  and  X  fixed  respec¬ 
tively,  i.e.  , 


(7.3)  J  =  det  (x^K  +  XkM>K=M) 

(7.4)  j  det  (XfcK) 


To  obtain  the  ratios  of  deformed  volume  elements  to  those  of  undeformed 
ones,  we  need  to  calculate  the  jacobians  J  and  j.  Since  the  determinant  of 
the  product  of  matrices  is  equal  to  the  product  of  the  determinants  of  matrices, 
we  have 


J  -  {det  [(x^K  +  xkM,K':M)  (xk,L  +  XkN,L“N^} 


1/2 


det  (xk,KXk,L  Xk,LvkM,K  M  +  Xk,KXkN,L "N  +  XkM,KXkN,L::M~N)  ' 


1/2 


Upon  using  (3.5)  and  (3,6)  this  becomes 


_1  1/2 

J  =  det  (CKL  +  ' LMK “M  +  XNL”N  +  fPMK' QNL  CPQ"M=N)  1 


(7.5) 
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where 


(7.6) 


CPQ  =  VkVk 


is  the  inverse  tensor  to  Cp^. 
right-har.d  side,  we  employed 


In  the  calculation  of  the  last  term  on  the 


XkL,M  =  FKLMXK,k 


which  follows  from  (3.7). 

The  above  result  (7.5)  is  valid  for  a  general  micromorphic  material 
and  can  be  further  simplified.  We  are,  however,  interested  only  in  a 
linear  theory.  In  this  case  (7.5)  can  be  approximated  by 

j  »  (dec  [SKL+  2EKL  +  2r(KU))1/2 
Expansion  and  linearization  of  this  gives 
O.7)  J  '  1  +  EKK  +  rKK 


Thus  the  macrovolume  change  with  fixed  5,  to  a  linear  approximation,  is 
given  by 


(7.8) 


dv 


1 


tr  E  +  (V  *  t)  '  ■ 


where  we  used  (4.9).  Here  the  first  term  on  the  right-hand  side  is  the 
classical  expression  of  the  dilatation, and  the  second  term  is  the  addi¬ 
tional  volume  change  due  to  microdeformation. 

Similarly,  we  can  calculate  the  minivolume  element  dv  by  determining 
j  given  by  (7.4).  In  this  case 

J  "  |det  <XkK*kL)ll/' 


(7.9) 
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For  the  linear  theory  from  (3.15)  we  have 

XkK  *6MK  +  ^MK^KM 
Hence  (7.9)  to  a  linear  approximation  is 
(7.10)  j  =  [det  (6rl  +  $RL  +  \K)]1/2  -  1  +  +  °(^2) 

For  the  micropolar  theory  4>  =  0  so  that  we  have 

K.K. 

- 1  -  °<»2> 

Hence ,  in  the  linear  micropolar  theory  there  will  be  no  minivolume  change. 
A  minivolume  (for  X  fixed)  rigidly  rotates  without  altering  its  value. 
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8.  COMPATIBILITY  CONDITIONS 


The  strain  tensor  and  micropolar  strain  tensors  and 
are  expressed  in  terms  of  the  displacement  field  and  microrotation 
field  $  by 


(8.1) 


ekt  •  2  <uk,l+  “l.k* 


(8.2) 


ek£  U£,k  +  F£km^m 


(8.3) 


Y  =  r  (1) 
k£m  kin  n,m 


Of  these, (8.1)  is  the  eulerian  linear  strain  tensor  known  to  us  from  the 
classical  theory  and  (8.2)  follows  from  (4.17)  by  use  of  (8.1)  and  (4.15) 
and  (8.3)  is  identical  to  (4.18).  When  the  six  quantities  u^  and  are 
prescribed,  these  strain  fields  are  determined  uniquely  through  (8.1)  to 
(8.3)  by  mere  substitution.  If  instead  the  six  strains  e^,  nine  micro¬ 
strains  ,  and  nine  non-vanishing  components  of  are  prescribed, 

then  the  determination  of  the  displacement  and  microrotation  fields  requires 
the  solution  of  twenty-four  partial  differential  equations  (8.1)  to  (8.3) 
for  the  six  unknowns  u^  and  Such  a  system  is  over-determined  and 

restrictions  must  be  imposed  on  e^,  ,  and  Yj^  •  These  conditions  are 

known  as  the  compatibility  conditions.  For  the  classical  strain  tensor  e^, 
the  compatibility  conditions  are  (cf.  Sokolnikoff  [1956]): 


(8.4/ 


G  +  e  —  Q  —  0  a  o 

k£,mn  mn,k£  km,£n  in, km 


If  we  notice  that 


(8.5) 


ek  =  e(k£)  2  (ck£  +  eik) 
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then  (8.4)  can  also  be  written  as 


(8.6) 


C (k£) ,mn  +  C (mn) ,k£ 


-  £ 


-  C  , 


(km) ,  £n  (in) ,km 


An  alternative  derivation  of  (8.6)  is  instructive:  To  this  end  chrough 
(8.2)  we  calculate  the  displacement  field  u^  from 

(8'7)  "k  '  “k  +  ^Cctk  +  <Wn,)dx« 

0  On 

where  is  the  value  of  u^  at  one  end  of  an  open  curve  C  in  the  body. 

Integrating  by  parts  the  second  term  under  the  integral,  we  also  write 


(8.8) 


\  +  e 


£km^n/X£ 


X£)  + 


/  l^,  +  £„,._(* 


'ik 


’  £km 


-  X?  )( 
£  £ 


m,lldXl 


For  the  displacement  field  to  be  independent  of  the  path  C  followed 
between  the  points  x^  and  x^,  the  integrand  must  be  expressible  as  a 
total  differential,  i.e., 


ik 


C£km(x£ 


m,  i 


where  F^(x,t)  are  single-valued  and  possess  continuous  partial  derivatives 
with  respect  to  x^,  through  second  order.  From  this  it  follows  that 

Fk,ij  =  Fk,ji 


Consequently 


c ik. j  ' 


Wl’.j  •  €Jk,i  +  [rekm<x£  -  VVj’.l  *  ° 


Expanding  the  second  term  and  using  (8.3)  we  obtain 


£ik,j  C j  k , i  +  Yikj  jki  ° 


(8.9) 
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These  are  the  necessary  and  sufficient  conditions  for  the  displacement  field 
u  to  be  single-valued  and  continuous  in  a  simply-connected  region.  Elimina¬ 
ting  y  from  (8.9)  by  differentiating  and  using  (8.2)  and  (8.3),  we  obtain 
(8.6)  again. 


A  similar  method  can  be  applied  to  Y,  «  as  follows.  First  we  solve  for 


4*  by  multiplying  both  sides  of  (8.3)  by  E  .  Hence 
IT  & in  kx,  r 


(8.10) 


>  «  —  £  Y 

r,m  2  k£r  kim 


Integration  of  (8.10)  along  a  smooth  open  curve  C  gives 

<f>  ■  /  e.  .  Yi  »  dx 

r  r  2  kir  kirn  m 

The  condition  of  single-valuedness  for  now  reads 


Thus  we  have  proved  the  following. 

Theorem:  The  necessary  and  sufficient  conditions  for  the  integrabilit 
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9.  SOME  SPECIAL  DEFORMATIONS 

In  this  section  we  present  a  lew  special  deformations  for  illustrative 
purposes  . 

(i)  Rigid  Deformation .  The  deformation  of  a  body  is  called  rigid  if 
the  distance  between  every  pair  of  points  X^  and  ^  in  the  body  remains 
unchanged.  From  (3.33)  it  is  clear  that  the  necessary  and  sufficient  condi¬ 
tion  for  the  rigid  deformation  of  a  linear  nicromorphic  body  is 

(9.1)  E  *  E  =  0  ,  [  =  0 

Alternatively,  in  terms  of  the  spatial  measures  of  strains 

(9.2)  e  =  e  =  0  ,  y  =  0 

For  a  micropolar  body  E  =  0  implies 

KL 

(9.3)  U(,  -  +  Br 

where  R^,  is  an  arbitrary  skew-symmetric  tensor  and 
both  of  which  are  independent  of  X.  The  condition  T 

K.L 

independent  of  X.  Finally  E  =  0  gives 

K.L 

(9.4)  \  ^ 

where 

(9.5)  ^  “  2  CKLMR>L 

is  a  rotation  vector  independent  of  X. 

(ii)  Isochorlc  Deformations.  The  deformation  will  be  called  macroisochoric 
if  the  material  maorovolume  remains  unchanged.  It  will  be  called  miniisochoric 


is  an  arbitrary  vector, 

=  0  implies  that  C>  be 
K 


if  the  minivolume  is  unchanged.  The  necessary  and  sufficient  condition  for 
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macroisochoric  deformations  according  to  (7.7)  or  (7.8)  is 

(9.6)  ekk  +  -  0 

The  condition  for  miniisochoric  deformation  follows  from  (7.10).  Since, 
for  the  linear  micropolar  bodies  $  =  0,  we  see  that  the  linear  micro- 

polar  bodies  undergo  only  miniisochoric  deformation.  We  note  that,  in 
general,  this  is  not  true  for  general  micromorphic  materials. 

For  the  condition  (9.6)  to  be  valid  for  all  i,  it  is  necessary  and 
sufficient  that 

(9.7)  ERK  "  0  ’  rKK-° 


(iii)  Homogeneous  Strain,  rhe  state  of  strain  in  a  body  will  be  called 
homogeneous  when  the  deformation  is  linear  and  homogeneous  in  the  position 
vectors  X  and  5  of  the  material  points,  i.e.  , 


(9.8) 


?<a)  ■  ?K*K  +  ?K;K 


where  D  and  P  are  constant  vectors.  In  terms  of  the  components  of  these 
-  K 

vectors,  (9.8)  is  equivalent  to  the  system 


(9.9) 


*k  =  DkRXK 


(9.10) 


Ck  =  ^kK~K 


where  and  P^  are  constants  for  static  deformations  and  functions  of  time 
only  for  the  dynamical  motions  and 


D,,  =  D,  i,  ,  P  =  Vi. 
~k  kk  k  k  kk  k 


(9.11) 
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Equation  (9.9)  is  the  expression  of  homogeneous  strain  in  classical  continuum 
mechanics.  The  deformation  described  by  this  set  carries  straight  lines  into 
straight  lines,  ellipses  into  ellipses,  ellipsoids  into  ellipsoids.  The 
microhomogeneous  deformation  (9.10)  is  new,  and  it  possesses  the  same  kind  of 
properties  with  respect  to  E.  For  the  micropolar  body,  D  can  be  replaced 
by  a  single  vector  V  given  through  the  solution  of 


(9.12) 


kL 


kL 


—  £  V 

kLM  M 


namely 


(9.13) 


V  =  --  £  o 

M  2  kLM  kL 


Equation  (9.12)  results  from  equations  (3.15)  and  (4.3)  with  $  =  V.  Thus  the 
microhomogeneous  deformation  can  also  be  expressed  by 


(9.14) 


E  -  E  *  D 


where  the  vector  V  is  independent  of  X  and  of  course  S. 

The  material  deformation  tensors  follow  from  (3.5)  to  (3.7)  and  (9.9) 
and  (9.10)  with  D  =  Thus 


(9.15) 


CKL  “  \KDkL 

y  =  n  ^  f 

KL  kK  kL  ’  KLM 


It  is  now  clear  that  the  strain  measures  are  homogeneous.  For  a  micropolar 
solid  we  use  (9.12)  so  that 


(9.16) 


KL 


°LK  “  ELMNDNK^M 


where  we  put 
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ULK  =  6kLDkK 


Below  we  give  several  special  cases: 

(iiia)  Uniform  Macrodilatation .  In  this  case,  D  is  a  diagonal  matrix 
having  the  same  entries,  i.e., 


(9.17) 


D  0  0 
0  D  0 
0  0  D 


0  <  D  <  00 


The  deformation  tensors  in  this  case  take  the  forms 


(9.18) 


CKL  °  5KL  *  ?KL  '  Dl?kL6kK  ’  ‘  KLM  =  ° 


The  deformation  carries  a  parallelepiped  having  edge  vectors  I^dX^  at  X  + 
to  one  with  edge  vectors  DdX^I^,  Fig.  9.1.  The  ratio  of  the  length  of  an 
edge  to  its  original  value  (the  macrostretch)  is  therefore  given  by 


l(k)  ■  D 


The  angle  between  any  two  edge  vectors  of  the  deformed  macroelement  is  90°. 

The  deformation  carries  a  macrocube  of  unit  volume  to  a  macrocube  of  volume 
3 

D  .  The  microelement  changes  according  to  the  values  of  V  .  From  (9.16), 

K 

for  this  case,  we  have 


(9.19) 


Y  .  o(5  -  r  V  ) 

KL  V  KL  KLM  M' 


From  this  it  is  clear  that 


<9-20)  ’ll  •  ’22  ■  ’ll  ■  “  •  ’KL  -  -CKLMU0M  -  (K  *  L) 


Thus,  the  microstretch  ^(^)»  ratio  of  the  edge  vector  of  the  deformed 
microelement  to  that  of  the  deformed  element,  is  given  by 


UP 
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(9.21) 


Hence,  the  edges  of  the  microelement  are  stretched  the  same  amount  as  those 
of  the  macroelement.  For  D  ■  1  we  have  no  macrodeformation  or  microdeforma¬ 
tion.  However,  the  microelement  undergoes  a  rotation  described  by  V.  For 
the  general  micromorphic  materials,  the  situation  is  much  more  complicated 
and  it  is  possible  to  have  microdeformations  even  when  the  macrovolume  remains 
unchanged.  This  situation  is,  of  course,  a  familiar  one  in  molecular  theories 
of  crystal  lattices. 

(iiib)  Uniaxial  Strain.  Consider  the  homogeneous  strain  characterized 
by 


(9.22) 


0  <  D  <  ® 


Macrostretch  L; 


and  microstretch  l, 


(9.24) 


L(l)  =  °  ’  L(2)  =  L ( 3)  "  1 
*(1)  =  °  ’  *(2)  =  *(3)  “  1 


For  0^  =  =  0 ,  we  have  the  classical  uniaxial  strain  according  to 

which  a  bar  of  length  dX^  after  deformation  becomes  a  bar  of  length  DdX^, 
Fig.  9.2.  For  non-vanishing  V  we  see  that  the  microelement  is  stretched 
ly  the  amount  D  in  the  X^  direction  in  addition  to  having  a  rotation.  The 
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geometry  becomes  particularly  simple  for  the  two-dimensional  case  for  which 
we  have 


(9.25) 


D  0 

i 

o 

M 

O 

_ 1 

,  T  = 

d  -£>  d] 

>  c  * 

3 

0  1 

L°  ij 

I 

o 

u> 

o 

h- 

1 - 

A  sketch  of  the  deformation  is  shown  on  Fig.  9.2.  The  macroelement  OACB 
after  deformation  becomes  OA'C'B  elongated  an  amount  (D-l)dX^  along  the 
X^-axis  with  no  change  in  the  lateral  directions.  The  microelement  Cacb 
is  stretched  by  the  amount  DdE^  in  X^  direction  becoming  Da'c'b'  After¬ 
wards,  a  rigid  microrotation  occurs  about  the  X^-axis  at  C.  The  final 
shape  of  Cacb  is  marked  by  Ca"c"b".  For  the  general  micromorphic  materials, 
it  is  possible  also  to  have  microstretches  independent  of  the  macrodeforma¬ 
tions  . 

Generally,  when  a  bar  is  stretched  in  one  direction  without  any 
constraints  on  its  sides  it  will  also  change  its  lateral  dimensions.  This 
situation  is  characterized  more  realistically  by 

D:  0  0 

0  d2  0 

0  o  d3 

In  this  case  the  state  of  strain  is  called  simple  extension .  For  simple 
extension  we  have 


(9.26) 


'of 

0 

0 

°1 

-°1P3 

°1P2 

(9.27) 

c  ■ 

0 

0 

,  a  = 

“lP3 

°2 

-°2P1 

0 

0 

4 

;“ip2 

°2P1 

V 
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From  (9.27)^  it  is  clear  that  a  macroelement  in  the  shape  of  a  rectangular 
parallelepiped  after  deformation  becomes  another  rectangular  parallelepiped 
with  its  sides  elongated  proportionally  to  D^,  D2»  and  D^.  A  microelement 
in  the  shape  of  a  rectangular  parallelepiped  changes  its  sides  in  the  same 
proportions,  however,  it  also  rotates. 

(iiic)  Simple  Shear.  In  classical  continuum  mechanics  a  homogeneous 
strain  characterized  by 

(9.28)  D  = 


ISO 
0  10 
0  0  1 


—  00  <  S  <00 


is  called  a  simple  shear.  Here  S  is  independent  of  X^.  The  spatial  position 
x(“}  of  any  material  point  X  +  5  after  deformation  is  given  by 

xja)  =  Xx  +  SX2  +  =!  +  p253  “  P3H2 


(9.29) 


(a) 


x2  +  -2  +  P3=1  -  P-^3 


x(a)  =x  +  =  +  V  -  -  V  - 

3  3  ”3  1~2  2'1 


In  the  case  of  V  =  0,  simple  shear  rotates  X.  =  const  planes  rigidly  about 
K.  1 

their  lines  of  intersection  with  X^  =  0-plane,  by  an  amount  equal  to  the 
angle  of  shear  y  given  by 


(9.30) 


arc  tan  S 


The  X2  =  const,  and  **  const,  planes  are  unchanged  (cf  Eringen  [1962,  Art.  15]) 
The  deformation  tensors  C,  Y,  and  1  are  given  by 
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(9.31) 


S  1+S  0 


0 


i 

-S3 

V2 

!  “ 

(s+i)P3 

1-SP3 

St>2-Vl 

V , 

1 

2 

1 

r  -  o 

From  (9.31)  we  see  that  even  when  there  is  no  macroshear,  that  is  when  S  ■  0, 
we  shall  have  a  microrotation  prescribed  by  the  microdeformation  tensor 

(9.32)  4-  - 


The  picture  for  a  plane  microdeformation,  in  this  case,  is  similar  to  the 
one  described  by  Fig.  9.2. 

(iv)  Plane  Strain.  In  classical  continuum  mechanics  when  the  defor¬ 
mation  of  a  body  is  identical  in  a  family  of  parallel  planes  and  vanishes 
in  the  directions  of  their  common  normal,  we  say  that  a  state  of  plane  strain 
exists.  This  plane  strain  is  thus  characterized  by 

(9.33)  *  xk(X1,X2)  ,  (k  =  1,2)  ,  x^  *  X3 


For  the  micropolar  deformations,  we  define 
similarly.  Using  (4.11)  and  setting  :  9 


-  <KXrX2)E 


C2  “  =2  +  *<X1,X2)5 


C3  "  =3 


the  state  of  plane -micropolar  strain 
0,  '  $  we  have 

2 

3 


(9.34) 
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Thus,  in  a  plane  macrostrain  and  microstrain  the  deformation  is  fully 
prescribed  when  the  two  unknown  displacements 


(9.35) 


u1(x1,x2)  =  ^  -  x2 

U2(X1,X2)  =  x2  -  X2 


and  the  microdisplacement  4>(X^,X2)  is  determined.  The  field  equations 
should,  therefore,  consist  of  three  partial  differential  equations 
replacing  the  two  equations  of  the  classical  theory.  This,  as  we  shall 
see,  is  the  case  (cf.  Arts.  25  and  26). 

The  spatial  deformation  tensors  for  the  plane  strain  are 


(9.36) 


C11 

C12 

0 

V 

11 

T 

12 

0 

C21 

C22 

0 

,  ?  c 

21 

Y 

22 

0 

_0 

0 

1_ 

_0 

0 

1_ 

ss  —  r  s  —  <t> 
12M  21M  »M 


(M  -  1,2) 


Equations  (3.27)  to  (3.29)  and  (4.3)  provide  the  relations  between  defor¬ 
mations,  tensors,  and  strains  (or  displacement  vectors): 


CKL  "  6KL  +  2EKL  “  6KL  +  UK,L  +  UL,K 


(9.37) 


4*  =  £  +6  -  6  -  c  $+U 

KL  KL  KL  KL  KLM  M  L,K 


T  «  -e  $ 

KLM  KLN  N,M 

Thus,  for  the  strain  tensors  E  and  we  have 

KL*  K.L 


E11 

E12 

0 

’En 

CM 
* — < 
LU 

0 

E21 

E22 

0 

.  §  " 

E21 

522 

0 

_0 

0 

0 

-J 

_0 

0 

0_ 

(9.38) 
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when 


V  B 

9Ui 

_1  ^2 
E21  ”  2  (3X2  +  3X1) 

3U„ 

F  E  £ 

L11 

3X1 

9 

E12 

*  22  3X2 

3U 

3U2 

(9.39) 

E11  " 

1 

3X1 

> 

'■12 

=  -$  +  — -- 

3XX 

E21  " 

4>  + 

3U1 

3X2 

9 

£  .22l 

22  3X2 
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10.  MOTION,  MICROMOTION,  MATERIAL  DERIVATIVE  OF  TENSORS 

A  material  point  at  X  +  E  at  time  t  ■  0  is  carried  to  a  spatial 
point  x  4-  5^°^  at  time  t.  The  motion  of  this  material  point  in  a  body 
is  described  by  one  parameter  family  of  transformations. 

(10.1)  x(a)  -  x(X, t)  +  aX.E.t) 

where  x(X,t)  is  the  place  occupied  by  the  center  of  mass  X  of  a  macro¬ 
volume  element  dV  +  dS  at  time  t  and  £(X,E,t)  is  the  relative  position 
of  the  point  X  +  E  at  time  t  with  respect  to  the  center  of  mass.  For 
micromorphic  bodies  we  have  (3.9)  or 

(10.2)  x(X,t)  ■  X  +  u(X,t) 

for  x(X,t)  and  (2.9)  for  £(X,E,t).  For  micropolar  bodies, (10.2)  remains 
valid  but  (2.9)  is  replaced  by  (4.11),  i.e., 

(10.3)  §(X,E,t)  =  E  -  E  *  $(X,t) 

In  these  expressions  U(X,t)  and  i>(X,t)  are  respectively  the  macrodisplace¬ 
ment  vector  and  microdisplacement  vector.  The  parameter  t  is  real  repre¬ 
senting  time. 

According  to  the  axiom  of  continuity  and  indestructibility  of  matter, 
the  inverse  motions  X(x,t)  and  E(x,£,t)  are  assumed  to  exist.  Thus  we 
may  also  write 

(10.4)  X(a)  -  X(x, t)  +  E(x,  ,t) 

where  -  is  given  by  (2.13)  for  micromorphic  bodies.  We  also  have 

=  x  -  u(x,t) 


(10.5) 


X(x, t) 
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However,  for  micropolar  bodies  (2.13)  is  replaced  by  (4.12),  t.e., 

(10.6)  5  =  §  +  §  x  $(x,t) 

For  displacement  vectors,  we  employ  the  same  symbol  u  and  $  both  in  material 
and  spatial  descriptions.  However,  in  the  former  case  u  and  <i>  are  assumed 
to  be  functions  of  X  and  t,  and  in  the  latter  case  functions  of  x  and  t 
since  we  may  substitute  X  =  X(x,t)  for  X  in  u(X,t)  and  $(X,t)  to  pass  from 
material  description  to  spatial  descriptions.  Single-valued  inverses  X(x,t) 
to  x(X,t)  and  E(X,£,t)  to  (  =  £(X,E,t)  are  assumed  to  exist  at  a  neighborhood 
of  X  at  all  times  except  possibly  some  singular  points,  lines,  and  surfaces 
in  the  body.  A  sufficient  condition  for  this  is  the  continuity  of  partial 
derivatives  of  these  functions  with  respect  to  X  and 

(10.7)  det  (j^  R)  i  0  det  (*R  R)  +  0 


in  some  neighborhood  of  X^  at  all  times.  We  assume  that  such  is  the  case. 

In  the  kinematics  of  continuous  media,  the  time  rates  of  vectors  and 
tensors  associated  with  material  points  play  an  important  role. 

Def .  1_.  The  material  derivative  of  any  tensor  is  defined  as  the 
partial  derivative  of  that  tensor  wi th  respect  to  time  with  the  material 
coordinates  X^  and  held  constant .  The  material  derivative  is  indicated 
either  by  placing  a  dot  on  the  letters  or  by  D/Dt .  Thus,  for  example, 

DF  3  F  ( X , t ) 

-~V-  s 

A 


(10.8) 


Df  3f  (x,t) 


D‘k 


Dt 


at 


X, 


k(x,t)  ~k 


where  subscripts  attached  to  a  bar  indicate  that  those  variables  are  held 
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constant  during  the  differentiation.  Since  through  the  centroidal  motion 
we  have 

(10.9)  xk  -  x^X.t)  or  \  “  Xj^x.t) 

we  have 


so  that 


fk(x,t)  -  fk(x(X,t)  ,t) 


5c 

3xt(X,t) 

at 

x  3*( 

t  3t 

X 


In  short,  without  ambiguity  we  write 


(10.10) 


Dfk  3fk 

fk  fk,^ 


The  first  term  on  the  extreme  right  of  this  equation  represents  the  time 

x 

rate  of  change  that  occurs  at  a  place^at  time  t.  The  second  group  of 
terms  is  known  as  the  convective  change.  These  arise  from  the  motion  of 
the  material  point  X  through  the  place  x. 

In  the  case  when  the  tensors  involve  the  macromotion  in  the  material 
differentiation,  we  also  consider  the  relative  location  vector  5  held 

Iv 

constant,  cf.  (10. 8) Another  example  is  provided  by  differentiation  of 


(4.11) . 

(10.11) 


§ 


ayx.t) 

*k  =  “ekLM"L 


X  $ 


or 


at 
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11.  VELOCITY,  ACCELERATION,  MICROROTATION,  SPIN 

Def .  _1 .  Veloci tv  Is  the  time  rate  of  ch ange  of  the  position  vector 
of  £  material  point.  Thus 

•  • 

(11.1)  v  =  x(X,t)  or  vk  “  \ 

where  X  is  held  constant.  Since  x  ■  x(X,t)  we  have 

3x(X,t) 

v  -  — - - -  x(X, t) 

Note  that  in  a  body  with  microstructure,  X  is  the  position  of  the  center 
of  mass  of  a  macrovolume  element  and  it  may  or  may  not  be  actually  occupied 
by  a  material  point.  Nevertheless,  in  defining  the  time  rates  of  vectors 
and  tensors  associated  with  the  body,  we  refer  to  X  as  the  material  point. 
Upon  replacing  X  by  (10.9)2  we  also  write 
•  • 

(11.2)  v  -  x(X(x,t) ,  t)  -  xk(x,t)ik  -  vk(x,t)ik 


where  ik  are  the  spatial  rectangular  unit  base  vectors.  This  equation 
defines  the  velocity  field  vk  at  a  spatial  point  x  at  time  t.  This  is  the 
eulerian  concept  of  the  velocity  field  which  is  prominent  in  hydrodynamics. 
For  the  lagrangian  viewpoint,  we  express  the  velocity  vector  in  the  material 
frame  of  reference  X^.  Thus 


(11.3) 


v  -  yx,t)iK  .  vK 


jyx.t) 

3t 


kK 


Def .  2^.  Acceleration  is  the  time  rate  of  change  of  the  velocity 
vector  of  a_  material  point.  Thus 


Dv 


(11.4) 


a  =  v 


or 


a.  -  v, 
k  k 


_k 

Dt 
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For  the  lagrangian  viewpoint  we  have 

3V  (X, t) 

(ii.5)  a  -  AR(X,t)IK  ,  AR  H  ~ - 


and  for  the  eulerian  viewpoint 

Dv, 


(11.6) 


ak  = 


_k 

Dt 


3vk  + 

■  +  v  V 

at  k,n 


Here  we  notice  the  appearance  of  the  convective  terms  v^  ^v^ 

The  velocity  v  and  acceleration  a  defined  above  are  the  kinematical 
quantities  describing  the  motion  of  the  center  of  mass  X  in  a  macromaterial 
element  V  +  S.  We  now  proceed  to  obtain  the  relative  velocity  and  accelera¬ 
tion  of  a  material  point  X  +  -  with  respect’  to  the  center  of  mass  X.  For 
these,  we  take  the  time  rates  of  the  relative  motion  for  a  micropolar  body 
given  by  (2.9),  namely, 


(11.7) 


§  -  xK(x,t)=R 


Thus 


(11.8) 


•  • 

§  ■  XR(X,t)EK 
•  •  •  • 

§  -  XK(X,t)EK 


Alternative  expressions  are  obtained  by  replacing  5  by  its  expression 
(11,9>  "K  "  XKk(',t)Ck 


Thus ,  for  example  • 


£  as  \)  £ 

T,  £k*k 


(11.10) 


<  - 


or 
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where 

(11,11)  -k(*,C^  =  XK(?,t)XKk(x,t)  ,  vak  =  X£K\k 


t 

It  is  understood  that  X  appearing  in  the  argument  of  is  also  replaced 
by  (10.9)2. 

Def .  _3.  The  three  vectors  defined  by  (11.11)  are  called  the 
gyration  vectors .  and  their  components  form  the  gyration  tensor . 

When  the  gyration  tensor  is  given,  we  can  calculate  the  eulerian 
microvelocity  §  by  (11.10).  For  the  microacceleration  in  a  similar  fashion, 
we  obtain 


Vk^k 


Vk 


Vk  + 


VkVk l 


where  we  used  (ll.lO)^.  Thus 


(11.12) 


?  -  ?k<*.tKk  or  Ek  -  »k,«* 


where 


(11.13)  a,  (x,t)  =  v.  +  v  v  ,  or  a.,  =  v„,  +  v.  v 

k  -m  mk  u*  ^  ■ 


ilk  Ik  Im  mk 


Def .  4^  The  three  vectors  a^fx.t)  defined  by  (11.13)  are  called 
the  spin  tensor. 

The  total  velocity  v^  and  acceleration  a^  of  a  material  point 
X  +  5  can  now  be  calculated  by 


(11.14) 


» 


X  +  §  =  v  +  vk£k 


(11.15) 


» 


a  +  a  =  v  + 


“kCk 
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For  a  micropolar  body,  these  expressions  are  modified  by  use  of  (A. 4)  and 
(4.14).  Thus,  for  example, 


(11.16) 


VkH  "  “ekJLM4>M  +  CkKHCK£mVm 


For  the  linear  theory  this  is  simplified  to 


(11.17) 


v.  .  -  -e,  . 
kt,  k£m  m 


Upon  introducing  an  axial  vector  v^,  called  microgyration  vector,  by 


(11.18) 


1 

v.  m  t  e,  „  v  .  ,  v,  .  =  -e,  .  v 

k  2  kirn  mS,  kl  k£m  m 


we  see  that 

(11.19)  vk  -  ^ 

and  (11.10)^  now  reads 

(11.20)  §  -  -§  x  y 

•  • 

Similarly,  one  can  calculate  §.  An  alternative  approach  that  may  be  instruc¬ 
tive  is  through  taking  the  time  rate  of  (11.20) 

•  •  ■  •  • 

(11.21)  §=-£xv-£xv--£xv+(£*v)  x  v 


If  we  recall  a  vector  identity 

(11.22)  (a  x  b)  x  c  ■  (a  •  c)b  -  (b  •  c)a 

the  above  expression  can  be  written  as 

•  •  • 

§  =  -C  x  v  +  (£  •  v) v  -  (v  •  v)£ 


(11.23) 
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whose  component  form  is 


(11.24) 


where 


(11.25) 


ki 


=  -e,  „  v  + 
kirn  m 


v,  v.  -  v  v  6,  . 
k  1  m  m  kl 


It  can  be  seen  that  this  is  identical  to  (11.13)2* 

The  total  velocity  and  acceleration  vectors  of  a  material  point 
X  +  S  in  a  micropolar  body  can  now  be  expressed,  respectively,  by 

(11.26)  v(a)  -  x(a)  -x+{-v-{M 

(11.27)  a^°^  ■  v(a)  ■a+£'«v-5xv+(£xv)  *v 
• 

Here  v  and  v  refer  to  the  centroidal  point  of  the  macrovolume  element, 
and  the  remaining  terms  on  the  extreme  rights  of  (11.26)  and  (11.27)  are 
the  relative  velocity  and  acceleration  about  the  centroid.  Equation 

(11.27)  can  be  linearized  further  by  dropping  the  triple  vector  product 


term  on  the  extreme  right. 
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12.  MATERIAL  DERIVATIVE  OF  ARC  LENGTH 

In  continuum  mechanics,  the  time  rates  of  arc  length,  elements  of 
surface,  and  volume  in  the  deformed  configuration  are  often  required. 

Here  we  prepare  the  groundwork  for  this,  while  at  the  same  time  introducing 
certain  new  concepts  essential  to  the  study  of  motion. 

Fundamental  Lemma  1_.  The  material  derivative  of  the  displacement 
gradient  is  given  by 


(12.1) 


Dt  “  vk,£Xl,K 


The  proof  of  this  is  immediate  since  D/Dt  and  9/3X^  can  be  exchanged,  i.e. 
D  Dxk 

Dt  lxk,K<*>t)1  ‘  <_Dt)tK'  '  vk,tXl,K 

where  we  used  x^  -  v^(x,t)  and  the  chain  rule  of  differentiation.  Another 
useful  expression  that  follows  from  (12.1)  by  multiplying  it  by  dX  is 


(12.2) 


dXk  ’  Vl4*! 


A  corollary  to  Fundamental  Lemma  1  is 


(12.3) 


Dt  (XK,k)  "  _XK,£VJl,k 


which  is  proved  by  differentiating  x^  ^  *  6^.  Thus 


*k,K  \,l  +  \,K  V*  a  0 


Now  multiply  this  by  X^  This  gives  (12.3). 
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Theorem  _1 .  The  material  derivative  of  the  square  of  the  arc  length 


ds  i.i  given  by 


(12.4) 


ds  ■  2dudxkdx* 


where 


(12.5) 


dk£  =  V(k,£)  =  2  (vk,£  +  V£,k) 


is  called  the  deformation  rate  tensor. 


To  prove  (12.5)  we  take  the  time  rate  of  ds  : 


<u2.f 

Dt 


(dW  • 2  drk  d\  ■  2vk,tdxkd,1t 


<vk,t  +  vt.k,dVxt 


Hence  the  proof. 

In  the  material  description,  (12.4)  can  be  written  as 

d“2  '  2dklVt.LdW. 

If  we  use 

ds2  -  '  (<KL  +  2EKL)dXKd\ 


then 


(12.6) 


d,!2  ’  'W\  ■  2EKLdXKd\ 


By  comparing  this  with  the  foregoing  expression,  and  since  d,  and  C  and 

K  It  KL 

are  symmetric  tensors,  we  find  that 
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(1Z,7)  CKL  “  2EKL  "  2dkiXk,KXH,L 

This  is  the  material  derivative  of  the  lagrangian  strain  measures  C  and  E. 

2 

When  d  -  0  we  have  D(ds  )/Dt  ■  0.  Conversely,  when  for  arbitrary 
2 

dx,  D(ds  )/Dt  -  0  we  must  have  d  ■  0.  Hence  we  have 

Theorem  2.  (Killing)  .  The  necessary  and  sufficient  condition  for 
the  macromotion  x(X,t)  be.  rigid  is  d  ■  0. 

Note  that  macrorigid  motion  does  not  imply  microrigid  motion.  As 
we  shall  see  below,  the  microelements  may  undergo  non-rigid  motions  even 
though  macroelements  may  be  moving  rigidly. 

Fundamental  Lemma  2.  The  material  derivative  of  the  micromotion  is 
given  by 

<12-8)  XkK(X,t)  -  vfcJlxiK 

This  result  follows  from  (11. 11)^  by  multiplying  it  by  x^L  and  using  (Z.ll)^. 
A  corollary  to  (12.8)  is 

U2.9) 


which  is  obtained  by  taking  the  material  time  rate  of  (2.11)  and  multiplying 


the  result  by  X. 


Ik’ 

Theorem  3..  The  material  derivative  of  the  microdisplacement  gradient 


Xj^L  is  Riven  b£ 


(12-10)  Dt  (xkK,L)  =  XkK,L  "  VkHX£K,L  +  Vki,mX)lKXm,L 


To  prove  this,  we  take  the  partial  derivative  of  (12.8)  with  respect  to 
and  exchange  D/Dt  and  3/3X^  since  this  is  permissible. 


*L 
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A  corollary  Co  this  theorem  is 


(12.11) 


<!’  -  v.  d£  +  v.  .  £„dx 

km  m  k£,m  £  m 


which  is  proved  by  taking  the  material  derivative  of 


(12.12) 


d^k  “  XkK,L“KdXL  +  XkKd=K 


and  using  (12.8),  (12.10),  and  CR  » 

Theorem  4_.  The  material  derivative  of  the  square  of  the 
(ds^)  is  given  by  ^ 

^((ds(<>))2]  -  [vkjt  +■  vt>k  +  (Vjl  +''lr,kHrldxkd,‘)l 


(12.13) 


+  2(vi,k  +  \£  +V,kCr)dxkd^  +  (V  +  v*k)d5kd5£ 


To  prove  thio,  we  take  the  material  derivative  of 


(a)  2 

(12.14)  (dsv  ')  -  dxRdxR  +  2dxkdCk  +  d^d^ 


and  use  (12.2)  and  (12.11). 

If  we  now  introduce  the  mi erode  format ion  rates 


(12.15) 

(12.16) 


bk£  =  Vk£  +  V£,k 


.  (a)  f(a) 

ak£m  Vk£,m  *  ak£  ’’kr.i.  r 


arc  length 


1 


Lringen  [1964c] 
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equation  (12.13)  can  be  expressed  in  the  form 


2 

to  l(d5<  '  >  1  ■  2|dk»  +  a<k£)  ]dXkdX4  +  2(bkt  +  °Kk)dlViSt 


(12.17) 


+  2Ib(kt)  -  dktldEkd{l 


In  a  region  of  a  micromorphic  body  when 


(12.18) 


d  *  0  ,  b  ■  0 


we  have  the  general  solution 


(12.19) 


\  ■  w +  bk  ’  \t m  “ki 


where  id^  is  an  angular  velocity  and  b^  is  a  velocity,  both  of  which  are 
independent  of  x,  and 


(12.20) 


“k!  +  “ik  '  ° 


Upon  substituting  (12.19)  into  (12.16),  we  see  that  a^^  *  0.  Conversely, 

(a)  ^ 

we  can  show  that  the  vanishing  d  and  b  is  also  necessary  for  D[(ds  )  ]/Dt»0. 
Hence  we  have  proved 


Theorem  5.  The  necessary  and  sufficient  conditions  for  a  micromor- 


(12.21) 


b  -  d  -  0 


This  theorem  replaces  the  well-known  Killing's  theorem  (Theorem  2  above)  for 
micromorphic  bodies. 

For  linear  micropolar  bodies,  considerable  simplification  is  achieved 
in  the  foregoing  results.  To  this  end  we  recall  (4.4)  and  (11.18).,  namely» 


(12.22) 


XkK  kK  "  kKMvM  *  kl  "£k)Un  n 
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Substitution  of  these  into  (12.8)  and  (12.10)  and  linearization  gives, 
respectively , 


(12.23)  x,  ,,  -  -e.  „  u  ,  Xi  w  ,  ~  ~c,  v  6  , 

kK  kKm  m  kK,L  kkm  m,r  rL 


The  expression  of  (12.11)  for  a  linear  micropolar  body  is 


(12.24)  d£,  =  -e,  „  d£.v  -  c,  .  v  £.dx 

k£m  *  tn  k£m  m,r  £  r 


or  in  vector  notation 


(12.25) 


d£  -  -d£  x  v  -  £  *  v,  dx 
-  -  -  -  r  r 


The  material  derivative  of  the  square  of  the  arc  length  for  this  case  is 


(12.26)  i  (<ds(a))2]  -  2[dkI  +  a(la)]dxkd*(  +  20,^  +  a^ldyl^ 


where 


(12.27) 


f ,  .  v 
kirn  m 


t,  v  «£ 
krm  m,£  r 
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13.  RATES  OF  STRAIN  MEASURES 


Def.  1.  The  i:ime  rates  of  various  strain  measures  are  the  same 


as  their  material  derivatives.  Thus,  for  example, 


DE 


KL 


De 


kS. 


KL  Dt 


(13.1) 


.  dekl 


'k  l  Dt 

De 


ta 


'KL  Dt 

Dr 


KLM 


'kl  '  Dt 
Dy, 


k£m 


KLM  ‘  Dt  *  'kim  '  Dt 


We  now  proceed  to  give  explicit  expressions  for  these  quantities. 
Theorem  _1 .  The  lagrangian  s train  rates  are  given  by 


(13.2) 

EKL  “  dk£Xk,KX^,L 

(13.3) 

• 

EKL  “  bk£Xk,KX£L 

(13.4) 

• 

fKLM  “  bk«-Xk,KX£L,M 

The  proof  of  (13.2)  has  already  been  given  in  Art.  12.  To  prove  (13.3), 
we  calculate  the  maferial  derivative  of 


(13.5) 


Hence 


EKL  \l  KL  *k,KXkL  "  6 


KL 


EKL  “  \l  "  Xk,K  XkL  +  *k,K  X 


kL 


Upon  using  (12.1)  and  (12.8)  we  obtain  (13.3).  The  proof  of  (13.4)  is 
constructed  similarly  by  taking  the  time  rate  of 
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(13.6) 

and  using  (12.1) 
Theorem  2 . 


rKLM  ' 

and  (12.10). 

The  eulerian  strain  rates  are  given  by 


(13.7) 

(13.8) 

(13.9) 


e,  „  =  d,  .  -  (e.  v  .  +  e  .v  ,  ) 
k£  ki  km  m ,£  mi  m,k 

e,  ,■  b.  ,  -  (t.  v  .  +  e  ,v  ,  ) 
ki  ki  km  mi  mi  m,k 

y.  *  -a,  „  +  £.  a  „  -  (y,  .  v  +  y,  v  .  +  y  .  v  ,  ) 

kirn  Heim  kr  rim  kir  r,m  krm  r£  rim  r,k 


The  proofs  of  these  are  somewhat  lengthy  and  will  not  be  given  here.  They 
are  obtained  by  differentiating  the  expressions  for  the  strains  and  using 
various  results  obtained  in  Art.  12.  For  the  proof  of  (13.7)  see  Eringen 
(1962,  Art.  22],  and  for  (13.8)  and  (13.9)  see  Eringen  [1967]. 

Equations  (13.2)  to  (13.4)  and  (13.7)  to  (13.9)  indicate  that  the 
eulerian  strain  rates  are  not  the  same  as  the  deformation  rates.  If  at 
time  t  the  medium  is  unstrained  and  the  motion  is  just  beginning,  we  can 
set  x  *  X  and  t  =  ■  so  that 


EKL(^0)  "  dk£4kK6£L 


(13.10) 


EKL(X’0)  =  bki  kK°£L 


KI.M 


(X,0) 


ak£m*kK6£L  m.M 


and 

eu (x-0)  ■  du 

(13-11)  ;kn(*-0)  ■  bk. 

Wx’0)  ‘ 
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For  the  infinitesimal  deformation  theory,  the  terms  enclosed  in  parentheses 
on  the  right-hand  side  of  (13.7)  to  (13.9)  can  be  neglected.  Therefore, 
in  this  case  (13.11)  should  be  valid  approximately  for  all  times,  i.e., 


(13.12) 


*  du 

‘ki(*-0  *  bkl 


Ykim 


(x,t) 


akim 
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14.  EXTERNAL  AND  INTERNAL  LOADS 

A  material  body  subject  to  the  external  and  internal  forces  undergoes 
a  deformation.  These  forces  may  be  of  mechanical,  electrical,  chemical, 
and  other  origin.  Here  we  are  only  concerned  with  the  mechanical  forces. 

In  the  particle  mechanics  of  Newton,  the  force  F  acting  on  a  particle  is 
considered  to  be  a  function  of  the  position  of  the  particle  x,  its  velocity 
v ,  and  time  t ,  i ,e.  , 

(14.1)  F  =  F(x,v, t) 

When  we  have  a  collection  of  particles,  then  for  each  particle  we  may  write 

(14.2)  F  =  F  (x  ,v  ,t)  ,  (a  =  1,2,...) 

"  ci  -  ci  ~  a  **  a 

In  a  volume  element  AV  of  a  continuum,  we  have  a  large  number  of  particles 
interconnected  with  such  forces.  If  the  particles  of  a  continuum  are  not 
free  to  move  independently,  then  the  interparticle  forces  are  balanced  among 
themselves  in  pairs.  This  then  places  restrictions  on  (14.2)  so  that  the 
number  of  independent  forces  is  much  smaller  than  the  free  collection  of 
particles.  In  conformity  with  these  restrictions  and  with  the  basic  postu¬ 
lates  of  continuum  mechanics,  the  forces  acting  on  a  body  are  resolved  into 
a  resultant  force  F  and  a  resultant  couple  M  given  by 

(14.3)  F  -  I  F  ,  M  =  J  x  xF 

L  -  Q  L  'Cl  -  Cl 

a  a 

The  first  of  these  equations  gives  the  vector  sum  of  all  forces  acting  on 
each  material  point  with  position  x^,  and  the  second  gives  the  vector  sum 
of  the  moments  of  these  forces  about  a  point  which  constitutes  the  origin 
of  x^.  In  a  continuum,  the  force  field  is  usually  considered  to  be  continuous 
and  (14.3)  may  be  replaced  by 
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(14.4)  F  *  J  dF  ,  M  *  J  (x  x  dF  +  dM) 

1/  1/  ' 

where  dF  is  the  force  density  at  a  point  x  and  dM  is  a  couple  density. 

This  latter  terms  arises  in  classical  mechanics  from  coupling  various 
particles  in  the  form  of  doublets  or  rigid  blocks  so  that  some  of  the 
applied  forces  on  particles  produce  a  couple  also.  This  physical  pic¬ 
ture  can  be  used  in  the  construction  of  the  theories  of  micromechanics 

From  a  continuum  viewpoint,  whatever  the  origin  may  be,  the  forces 
and  couples  may  be  divided  into  three  categories. 

a.  Extrinsic  Body  Loads .  These  are  the  forces  and  couples  that 
arise  from  the  external  effects.  They  act  on  the  mass  points  of  the 
body.  They  appear  in  the  form  of  body  forces  and  body  couples  per  unit 
mass  of  the  body.  The  force  of  gravity  is  an  example  of  a  body  force, 
and  an  electromagnetic  moment  in  a  polarized  medium  is  an  example  of  a 
body  couple.  A  body  couple  can  also  arise  from  the  uneven  distribution 
of  the  mass  among  microvolume  elements,  Fig.  14.1  and  Fig.  14.2. 

b.  Extrinsic  Surface  Loads  (Contact  Loads) .  These  loads  arise  from 
the  action  of  one  body  on  another  through  the  contacting  surfaces.  At  a 
small  macrosurface  they  are  equipollent  to  a  force  and  a  couple.  Thus, 
for  example,  forces  acting  on  a  macrosurface  of  Fig.  14.3  are  equipollent 
to  a  force  and  a  couple  placed  at  the  centroid  of  the  ma  rosurface  element 
Aa,  Fig.  14.4. 

When  a  macrovolume  Av  is  allowed  to  tend  to  dv,  in  general  the  body 
couple  L  vanishes  since  the  moment  arms  cf  forces  tend  to  zero  while  the 
forces  are  assumed  to  remain  bounded.  Similarly,  when  Aa  is  allowed  to 
tend  to  da,  the  surface  couple  M  will  approach  to  zero.  This  is  the  classical 
picture  in  continuum  mechanics 
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Because  of  the  granular  nature  of  the  bodies,  the  mathematical  limits 
dv  and  da  for  the  surface  and  volume  elements  are  approximations  which  may 
not  be  admissible  for  various  physical  phenomena  in  which  the  applied  loads 
produce  effects  with  some  typical  lengths  (e.g.,  wave  length)  that  are 
comparable  to  distances  and  sizes  of  the  microelements.  In  such  situations, 
Av  and  Aa  are  not  infinitesimals,  and  the  granular  nature  of  the  bodies  must 
be  taken  into  account  in  some  form.  This  then  requires  that  we  consider  the 
existence  of  both  forces  and  couplesfor  macrovolume  and  surface  elements. 

c.  Internal  Loads.  Internal  loads  arise  from  the  mutual  action  of 
pairs  of  particles  that  are  located  inside  the  body.  According  to  Newton's 
third  law,  the  interparticle  forces  cancel  each  other  so  that  the  resultant 
force  is  zero. 


In  continuum  mechanics,  the  internal  effects  are  found  by  isolating  a 
small  macroelement  from  the  body  and  considering  the  effect  of  the  rest  of 
the  body  as  forces  and  couples  on  the  surface  of  the  macroelement  as  illus¬ 
trated  in  Figs.  14.1  to  14.4.  Internal  forces  give  rise  to  the  stress  and 
couple  stress  hypotheses,  as  we  shall  see  below. 

Let  the  surface  force  per  unit  area  at  a  point  x  on  the  surface  of  a 


body  having  exterior  normal  n  be  denoted  by  t^,  and  the  surface  couple 
per  unit  area  by  m^  .  Let  the  body  force  and  body  couple  per  unit  mass  at 
an  interior  point  of  the  body  be  respectively  represented  by  f  and  l.  The 
total  force  F  and  the  total  couple  M  about  a  point  0  acting  on  the  body 


can  be  calculated  by,  Fig.  14.5 


(14.5)  F  ■*  Set/  sda  +  /  ofdv 

(n) 

Concentrated  loads  are  imagined  as  resulting  from  a  limiting  process  in  which 
the  surface  loads  or  body  loads  are  distributed  over  a  very  small  region. 


]da  +  j  p  (S,  +  x  *  f)dv 
\! 


(14.6)  M  ■  +  x  *  t 
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15.  MECHANICAL  BALANCE  LAWS 


The  mechanical  balance  laws  -  conservation  of  mass,  balance  of 
momentum,  and  balance  of  moment  of  momentum  -  are  obtained  by  a  process 
of  averaging  applied  to  a  macrovolume  element  containing  N  microelements 
for  which  the  classical  balance  laws  are  postulated  to  be  valid.  Each 
of  the  microelements  is  considered  to  have  a  uniform  mass  density. 

I .  Principle  of  Conservation  of  Mass .  The  total  mass  of  each  micro¬ 
element  remains  constant  during  any  deformation.  Thus,  x  and  P ^ 

respectively  denote  the  mass  densities  of  a  microelement  a  before  and  after 
deformation,  and  AV^  and  Av^  their  volume,  Fig.  15.1,  then 

(15.1)  p^AV^  “  P^Av^  ,  (a  not  summed) 


The  total  mass  of  a  macrovolume  before  and  after  deformations  is  respectively 
given  by 


(15.2) 


P0AV0 


N 

V 

L 

a-1 


avJ“) 


(15.3) 


04  5  l  0  <a>4V <a) 

a*l 


These  equations  in  effect  define  the  mass  densities  P0(X)  and  p(x,t)  of  the 
undeformed  and  deformed  macrovolume  elements.  In  view  of  (15.1)  we  see  that 


(15.4) 


P^AVq  =  oAv 


If  we  let  AV(  and  AV  approach  their  limiting  values  dV  and  dV,  then 


0ndV 
0  0 


pdv 


(15.5) 


or 

dv 

T  '  dv"  J  ‘  det  (xk,K> 
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which  are  the  equivalent  expressions  of  the  principle  of  conservation  of 
mass  for  the  macrovolume  element. 

In  Section  2  we  said  that  X  is  the  position  vector  of  the  center  of 
mass  of  a  macroelement.  Accordingly, 

[C<“)H<“Va)  -  0 

a 


Upon  using  (15.1)  and  (2.14)  this  gives 

■  0 
a 


Since  4  0,  this  shows  that  the  position  vector  x  is  the  center  of  mass 
of  the  deformed  macrovolume.  Consequently, 

Theorem  1_.  The  motion  carries  the  center  of  mass  of  the  undeformed 
macrovolume  to  the  center  of  mass  of  the  deformed  macrovolume* . 

Nect,  we  calculate  the  second  moments 


(15.6) 


(a) _(a) _(a)  (a) 

0  "K  "L  0 

* 


upon  substituting  (15.1)  and  (2.13),  this  may  be  written  as  (cf.  Eringen 


[1964c]) 

(15.7) 


*KL  “  ikHXKkXL 


where 

/ 1  c  o\  .  -  v  (a) .  (a) .  (a)  (a) 

(15.8)  pi^Av  =  i  t,  Av 

a 

Quantities  1^  and  i^  are  respectively  called  the  material  and  spatial 
microinertia  tensors.  Equations  (15.7)  may  be  stated  as 


1 


Eringen  [1964c] 
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Theorem  2_.  The  microlnertla  is  conserved,  i.e., 


(15.9) 


Dt  "  0 


Using  (12.9),  this  may  also  be  expressed  as  (Eringen  [1964c]) 


— —  +  i.  v  -  i,  v.  -  i  .v.  **0 

3t  k£,m  m  km  £m  m£  km 


In  a  micropolar  continuum,  a  combination  of  i,  and  I  occur  more 

K  K.i_i 


frequently.  These  are 


(15.10) 


JKL  =  IHM6KL  "  :KL 


j,  i  i  6.  -  i.  . 

Jk£  mm  k£  k£ 


These  tensors  are  identical  to  the  inertia  tensor  encountered  in  rigid- 
body  dynamics. 

Upon  linearization  and  using  (4.14),  (15.7),  and  (15.11),  we  get 


(15.11) 


JKL  ^  Jk£6Kk6L£ 


Global  balance  equations  for  mass  and  microinertia  are  obtained 
by  integrating  (15.4)  and  (15.7)  over  the  volume  of  the  body.  Thus 


(15.12) 


/ p  dV  ■  / pdv 

y  0  0 


(15.13) 


/°0IKL<iV0  ’  /,0iktXKkXUdV 


where  V  is  the  undeformed  volume  and  U  is  the  deformed  material  volume. 


II.  Principle  of  Balance  of  Momentum.  The  time  rate  of  change  of  momentum 


is  equal  to  the  sum  of  all  forces  acting  on  a  body. 
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/fl\ 

The  mechanical  momentum  of  a  microelement  Av'  is  the  product  of 
its  mass  with  the  velocity,  namely,  p  ^Av^  .  The  total  momentum  of 
a  macroelement  is  the  vector  sum  of  the  micromomenta  of  its  microelements. 
For  a  micropolar  body  we  have 

.  p  (a)  (a)  (a)  r  (a),  ,  L.  .  <<*> 

Ag  ■  Zp  Y  Av  ”  Zp  (y  +  C)Av 

a  a 

■  y2p(‘‘,Av(‘‘,  +  v  x  £p<“Wo) 
a  a 

The  last  term  vanishes  and  in  the  limit  we  write 

dp  ■  pvdv 

The  total  momentum  of  the  body  is  therefore  given  by 

(15.14)  p  -  / pvdv 

'  1 1  ' 

The  principle  of  balance  of  momentum  is  expressed  by 

(15.15)  /pvdv  -  kt,  .  •  da  +  /pfdv 

dt  '  §'(n)  '  ' 

Here  is  the  surface  traction  per  unit  area  acting  on  the  surface  of 

the  body  S  with  an  outward  directed  normal  n  so  that  the  surface  integral 
is  the  vector  sum  of  all  forces  acting  on  S.  The  vector  sum  of  the  body 
forces  is  given  by  the  volume  integral  on  the  right.  Equation  (15.15)  is 
none  other  than  that  given  in  classical  continuum  mechanics. 

Ill .  Balance  of  Moment  of  Momentum.  The  time  rate  of  change  of  moment 
of  momentum  about  _a  point  is  equal  to  the  sum  of  all  couples  and  the  moment 
of  all  forces  about  that  point . 
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The  mechanical  moment  of  momentum  of  a  microelement  is  defined  as 
the  moment  of  its  momentum,  namely, 

(a)  (a)  (a)  (a) 

x  x  p  v  Av 

The  total  moment  of  momentum  of  a  macroelement  is  calculated  by 

c  (°)  (<*)  ( a ) .  (a) 

dH  -  2?  x  p  Y  “v 

a 

-  [(x  +  §)  x  P(a)(y  +  OAv(a) 
a 

Upon  carrying  out  the  multiplication,  we  get 

Am  ■  x  *  v£p(a)Av(a)  +  x 
a  a 

+  s  «  Ip(a)54v<a>  -  v  «  Io(“,{4v(“) 
a  a 

The  last  two  summations  vanish  since  £  is  measured  from  the  center  of  mass 
of  the  deformed  macroelement.  Upon  carrying  the  expression  £  from  (11.20), 
in  the  limit  we  may  write 

(15.16)  dm  -  px  x  vdv  +  podv 

where 

(15.17)  poAv  5  £p(a)S  x  (v  x  pAv^ 

a 

is  called  the  intrinsic  spin .  In  component  form  thi3  reads 

(15. 18)  V  Jk£\ 

where  we  used  (15.8)  and  (15.11)2  after  expanding  the  triple  product. 

The  total  moment  of  momentum  of  a  macroelement,  therefore,  is  the 
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vector  sum  of  its  angular  momentum  and  the  intrinsic  spin.  The  total  moment 
of  momentum  of  a  micropolar  body  is  now  calculated  by 


M  -  /(x  x  pv  +  po)dv 
1 1  ' 


(15.19) 

The  principle  of  moment  of  momentum  is  expressed  by 


^  /(?  x  pv  +  pq) dv  -  x  +  m^Jd 


l / 


h' 


(n)  "‘(n)' 


(15.20) 


+  fp(l  +  x  *  f)dv 
1/ 


The  right-hand  side  gives  the  sum  of  all  moments  about  the  origin  as  in  (14.6) . 
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16.  STRESS  AND  COUPLE  STRESS 

The  state  of  Internal  loads  and  their  connections  to  surface  loads 
may  be  found  by  applying  the  principles  of  global  balance  of  momenta  to 
small  regions  fully  and  partially  contained  in  the  body,  Fig.  16.1.  To  this 
end,  we  first  consider  a  small  macrovolume,  v  +  s,  fully  contained  in  the 
body.  At  a  point  x  of  s,  the  effect  of  the  remainder  of  the  body  is  equi¬ 
pollent  to  a  surface  force  per  unit  area,  t^,  called  the  stress  vector, 
and  a  couple  per  unit  area,  ,  called  the  couple  stress  vector.  These 

loads  depend  on  the  position  x,  time  t,  and  the  orientation  of  the  surface 
s  at  x  which  is  described  by  the  exterior  normal  n  to  s  at  x.  This  latter 
dependence  can  be  found  explicitly  by  applying  the  mechanical  principles  of 
momenta  to  a  region  v  +  s  adjacent  to  the  surface  of  the  body  This  approach 
does  in  fact  also  provide  the  connection  of  surface  loads  to  the  internal 
loads.  Consider  a  small  tetrahedron  with  three  faces  taken  as  the  coordinate 
surfaces  and  the  fourth  face  being  a  part  of  the  surface  of  the  body,  Fig. 
16.2.  We  denote  the  stress  vectors  on  any  coordinate  surface  x^  =  const,  by 
-t^  and  on  s  by  t^.  The  equation  of  balance  of  momentum  (15.15)  can  be 
applied  on  this  tetrahedron.  Using  the  mean-value  theorem  to  estimate  the 
volume  and  surface  integrals,  we  write 

d  *  *  *  *  * 

-  (pv  Av)  -  t(n)Aa  -  tkAak  +  pf  Av 

where  the  quantities  marked  with  asterisks  are  the  values  of  those  without 
asterisks  at  some  points  of  v  +  s  The  volume  is  denoted  by  Av  and  the  sur¬ 
face  areas  by  A a.  and  Aa.  The  mass  is  conserved  so  that 


(16.1) 


“7"  (pdv)  =  0 
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Upon  dividing  both  sides 
and  Av  approach  zero,  we 

(16.2) 

The  four  surfaces  of  the 
limit  of  the  sum  of  area 

(16.3) 


of  the  foregoing  equation 
see  that  Av/Aa  •+  0  and  we 

hn)da  ‘  Vak 

tetrahedron  form  a  closed 
vectors  da^  must  add  up 

da  ■*  nda  ■  da^i^ 


by  Av  and  letting  Aa 
obtain 

surface,  therefore,  the 
to  da.  Hence 


From  this  we  get 

(16.4)  da^  »  n^da 


Substituting  this  into  (16.2)  we  get 
(16'5>  -(n)  '  W 

where  t^  is  independent  of  n.  Thus  we  found  that  the  stress  vector  t^  is 
a  linear  function  of  n.  At  two  sides  of  a  surface,  n  changes  sign.  From 

(16.5)  we  therefore  see  that 


(16.6) 


E(-n)  ‘  ~'(n) 


which  proves  that  the  stress  vectors  on  opposite  sides  of  the  same  surface  at 
a.  given  point  are  equal  in  magni  tude  and  opposite  in  sign . 

The  application  of  the  above  method,  with  the  use  cf  the  equation  of 
balance  of  moment  of  momentum  (15.20)  leads  to 


(16.7) 

%) 

=  Vk 

(16.8) 

m(-n) 

E  'm(n) 
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The  concepts  of  stress  tensor  t^  and  couple  stress  tensor  m^  now  follow 
from  the  decompositions 


(16.9) 


t,  ■  t,  „i„ 
~k  k  I'i. 


(16.10) 


,  th 


Thus  t^  is  the  i  component  of  the  stress  vector  t^  which  acts  on  the 
surface  x^  ■  const,  and  m^  is  the  component  of  the  couple  stress  vector 
which  acts  on  the  same  surface.  The  positive  directions  of  t^  and  those 
of  m^  are  shown  on  Figs.  16.3  and  16.4  respectively.  We  use  double-headed 
arrows  for  m, 

From  (16.5),  (16,7),  (16.9),  and  (16.10)  it  follows  that 


(16.11) 


'(n) 


(16.12) 


”(n)  “  Vnk'f 


It  is  thus  clear  that  the  moment  vectors  for  the  couple  stress  have 
the  identical  sign  convention  to  those  of  the  stress  vectors.  The  plane  of 
each  couple  is  of  course  perpendicular  to  the  couple  vector,  and  the  direc¬ 
tion  is  ns  described  by  the  right-hand  screw  rule. 

The  expanded  form  of  the  components  of  t ^  and  m^  in  rectangular 
coordinates  are 


s 

t  n 

+ 

t  n 

+ 

t  n 

(n)x 

XX  X 

yx  y 

zx  z 

_ 

t  n 

+ 

t  n 

+ 

t  n 

(n)y 

xy  x 

yy  y 

zy  z 

s 

t  n 

+ 

t  n 

+ 

t  n 

(n)z 

XZ  X 

yz  y 

zz  z 

(16.13) 
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(16.14) 


m  \  “ 

m  n 

+ 

m  n 

+ 

m  n 

(n)x 

XX  X 

yx  y 

zx  z 

m  /  N  “ 

m  n 

+ 

m  n 

+ 

m  n 

(n)y 

xy  x 

yy  y 

zy  z 

m(n)  - 

m  n 

+ 

m  n 

+ 

m  n 

-  z 

XZ  X 

yz  y 

zz  z 
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17.  LOCAL  BALANCE  LAWS 

Local  balance  laws  are  obtained  by  postulating  that  the  global 
balance  laws  are  valid  for  every  part  of  the  body.  For  the  conservation 
of  mass,  we  convert  the  volume  integral  over  V  to  V.  This 

(17.1)  /(p0  -  pJ)dv  -  0 

V 


where 


(17.2) 


J  =  det 


(,k,K> 


is  the  Jacobian  of  the  transformation.  Postulating  that  (17.1)  is  valid 
for  every  part  of  the  body,  we  obtain  the  equation  of  local  mass  conser¬ 
vation. 


(17.3) 


co/o  -  J 


Another  form  often  used  in  hydrodynamics  is  obtained  from  this  by  taking 

the  material  derivative  of  (17.3).  Thus 

•  • 
pj  +  pJ  ■  0 

and  we  can  show  that  (Eringen  [1962,  Art.  19]) 


J  "  Jvk,k  ’  0 


Consequently 


0  +  0Vk_k  .  0 


'  8p  , 

P  -  +  P.kvk 


or  since 
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this  reads 


(17.4) 


ft  +  (pV*k  ‘  0 


This  is  the  celebrated  equation  of  continuity.  One  may,  of  course,  equally 
employ  (17.3)  in  place  of  (17.5).  Alternative  forms  which  follows  from  (15.4) 


are 


(17.5) 


p.dV.  -  pdv  ,  — ■ —  =  0 

0  ’  pdv 


An  expanded  form  of  (17.4)  in  rectangular  coordinates  is 


(17.6) 


i£.+ 

3t  3x 


3(pv  )  3(pv  )  3(pv  ) 


3y 


3z 


where  (v^,  v  ,  v^)  are  the  rectangular  components  of  the  velocity  field. 

Equations  of  local  balance  for  the  microinertia  are  already  given  by 

(15.7)  and  (15.9).  For  the  linear  theory  using  (15.11)  we  get 

DJ, 


(17.7) 


Dt 


0 


The  local  balance  of  momenta  follows  from  (15.15)  and  (15.20).  Upon 

carrying  out  the  indicated  differentiation  and  using  (17. 5)^,  and  writing 
•  * 

a  *  v,  v  *  x,we  see  that 


(17.8) 


/ padv  =  it,  .da  +  /pfdv 
V  '  s(t?)  V 


(17.9)  /(x  x  pa  +  po)dv  =  i (x  x  t.  .  +  m.  .)da  +  fp(l  +  x  *  f)dv 

y  g  y 


These  are  other  forms  of  the  global  balance  of  momenta.  We  now  take  V  +  S 

from 

to  be  a  small  internal  portion  v  +  s  of  the  body.  Subs  ti  tutinfy|  (16 . 5) 

(16.7)  we  write 


and 
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(17.10) 

Jpadv  ■  ^t^n^da  +  Jpfdv 

vs  V 

(17.11)  / p(x  x 

V 

g  +-  o)dv  -  $(x  x  t^  +  m^)n^da  +  J p(l  +  x 
s  V 

x  f)dv 

In  rectangular  coordinates,  the  Green-Gauss  theorem  is  expressed  as 

(17.12) 

Vkda  ■  /?t.kdv 

s  V 

If  we  now  apply  this  theorem  to  (17.10)  and  (17.11)  to  convert  the  surface 

integrals  to  volume  integrals,  we  obtain 

(17.13) 

/[tk  k  +  p(f  -  a)  ]dv  =  0 
v  ’ 

(17.14) 

/ln>k  k  +  4k  x  tk  +  p(l  -  q)  ]dv 
v  '  ’ 

+  /x  x  [tk  k  +  o(f  -  a)]dv  =  0 

v~  ’ 

For  these  equations  to  be  valid  for  any  arbitrary  volume  v 

the  necessary  and  sufficient  condition  is  the  vanishing  of 

Hence 

in  the  body, 

the  integrands. 

(17.15) 

tk,k  +  ( (f  -  v)  =0 

(17.16) 

^.k  +  ik  x  +  P(4  "  c)  =  0 

Note  that  the  second  integrand  in  (17.14)  vanishes  by  virtue  of  (17.15). 

These  equations  are  the  expressions  of  the  local  balance  of  momenta.  They 

are  identical  to  those  given  in  Eringen  [1962,  eqs .  (32.7)  and  (32.8)]  with 

the  exception  of  the  spin  inertia  term  pa.  This  term  arises  from  the  postulate 
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of  an  independent  microrotation.  In  fact,  without  such  an  internal  degree 
of  freedom,  the  existence  of  m^  and  £  is • questionable .  Upon  substituting 
(16.9)  and  (16.10)  into  (17.15)  and  (17.16),  we  obtain  the  component  form 
of  these  equations,  namely, 


(17.17) 


'£k,£ 


+  P(fk  - 


V 


0 


(17.18) 


m„.  „  +  e,  t  + 

£k,  kmn  mn 


P(£,  -  o,.) 


0 


These  are  the  first  and  second  laws  of  motion  of  Cauchy  which  express  the 
local  balance  of  momenta  for  micropolar  bodies^-.  When  the  body  is  nonpolar, 
that  is,  when  o  =  =  £  *  0,  then  (17.18)  gives  the  classical  result 


(17.19) 


tk£  =  C£k 


which  expresses  the  symmetry  of  the  stress  tensor.  For  micropolar  bodies, 
we  see  that  the  stress  is  in  general  nonsymmetrical  and  the  new  set  of  differ¬ 
ential  equations  (17.18)  must  be  employed  replacing  (17.19). 

In  rectangular  coordinates,  the  expanded  expressions  of  (17.17)  and 
(17.18)  are  recorded  below 


(17.20) 


3t  3t  3t 

xx  +  yx  zx 

3x  3y  3z 

3t  3t  3t 

,*y  +  -22.  +  _ ?2L 

3x  3y  3  z 

3t  3 1  at 

xz  vz  t  zz 

3x  3y  dz 


+  P(f 

X 


+  p(f 

y 


•*  o(f 


z 


V 


V  ) 

z 


0 


0 


0 


If  we  disregard  the  relation  of  to  given  by  (15.18),  these  balance 
laws  are  valid  for  the  nonlinear  theory,  and  (17.17)  and  (17.  .8)  are  exact 
expressions,  cf.  Eringen  and  Suhuoi  [1964a  &  b]. 


1 
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t  +  p(2.  -  o  ) 

zy  xx 


(17.21) 


9m  9m  9m 

*y  + _ xx.  + _ si 

9x  9y  9z 


i x2  +  P^y  -  °y> 


t  +  p(fc  -  o  ) 
yx  z  z 
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18.  CONSERVATION  OF  ENERGY 


An  Important  axiom  of  thermomechanics  is  the  principle  of  conservation 
of  energy  which  may  be  stated  as  follows: 

IV.  Principle  of  Conservation  of  Energy .  The  time  rate  of  change  of  the 
sum  of  the  kinetic  energy  and  the  internal  energy  is  equal  to  the  sum  of  the 
mechanical  energy .  heat  energy .  and  other  energies .  Here  we  exclude  chemical 
and  electrical  energies  so  that  we  may  express  this  law  mathematically  as 

•  • 

(18.1)  K  +  E  -  ID  +  & 


Here  K,  E,  ill  and  Q  are,  respectively,  the  kinetic  energy,  the  internal  energy, 
the  work  of  applied  loads  per  unit  time,  and  the  heat  energy.  For  a  micro- 
polar  continuum,  these  quantities  may  be  expressed  as: 

(18. 2>  K-i/°(Vk  +  JktVt>dv 

(18.3)  E  «  /pedv 

1/ 


(18. A) 


111 


f(t 

S 


lkvk  + 


mHkVk)d\  + 


/p(f 


kvk  +  Vk)dv 


(18.5) 


Q.  “  taida,  +  / phdv 

s  V 


The  physical  meaning  of  some  of  the  terms  occurring  in  these  equations  is 
known  to  us  in  the  classical  continuum.  For  example,  the  first  term  in  the 
integrand  of  (18.?.)  is  the  kinetic  energy  of  the  macromotion.  The  second 
term  is,  however,  new  and  it  is  the  kinetic  energy  of  the  microrotation  In 
(18.3),  e  is  the  internal  energy  density  per  unit  mass.  In  (18. A),  the  surface 
integral  is  the  work  of  surface  tractions  and  surface  couples  per  unit  time, 
while  the  volume  integral  is  the  work  of  the  body  force  and  body  couple,  per 


unit  time. 
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Finally,  in  (18.5)  the  surface  integral  gives  the  heat  input,  and  the  volume 
integral  the  heat  source. 

The  equation  of  local  energy  balance  is  obtained  by  postulating  that 
(18.1)  is  valid  for  any  arbitrary  volume  contained  in  the  body.  To  this 
end  we  first  carry  out  the  indicated  differentiation  with  respect  to  time 
for  K  and  E.  Thus 

(18.6)  K  -  /p(akvk  +  okvk)dv 

(18.7)  E  “  fpc dv 

1/ 


where,  in  anticipation  of  the  local  laws,  we  employed  the  equations  of 
conservations  of  local  mass  and  inertia  (17.5)2,  (17.7),  namely, 


(18.8) 


Dt 


0 


Next  we  convert  the  surface  integrals  of  (18.4)  and  (18.5)  into  volume 
integrals  by  use  of  the  Green-Gauss  theorem.  Hence 


W 


j,(ttkvk,t  + 


"tk\,t)du  + 


^k.t  +  P£k)vk  + 


<m«k.« +  °Vvkldv 


2  *  /(%  k  +  »h)dv 
l / 


Upon  carrying  (18.6),  (18.7),  and  the  above  equations  into  (18.1),  and  using 
the  equations  of  local  balance  of  momenta  (17.17)  and  (17.18),  we  obtain 


/(pc  -  t..  v  „  +  e  t  v  -  m..v  -  q  -  ph)dv  =  0 
y  £k  k,£  kmn  mn  k  ik  k,f  k,k 


This  is  assumed  to  be  valid  for  every  part  of  the  body.  Thus  we  must  have 
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(18.9)  pc  -  clkvk>t  -  ckmntm„vk  +  mlkvk_t  +  qk>k  +  Ph 


This  is  the  differential  equation  of  the  local  balance  of  energy  of  a  micro- 


polar  body  .  In  expanded  form  it  reads 


3v  3v  3v 

•  _ X  _ X  ,  .  _ X 

C  e  Cxx  3x  yx  3y  zx  3z 


+ 1  — — +  t  —r~  +  t  — u 

xy  3x  yy  3y  zy  3z 


3v  3v  3v 

_ z  _ z  z 

Cxz  3x  tyz  3y  Czz  3z 


(18.10) 


-(t  -t  )v  -  ( t  -t  )v 
yz  zy  x  zx  xz  y 


3v  3v  av 

_ X  ,  _ X  _ X 

+  m  ~~r~  +  m  t  t  m  “ 
xx  3x  yx  3v  zx  3z 


(t  -  t  )v 
xy  yx  z 


+  m  — +  m  — ~  +  m  -r~ 

xy  3x  yy  dy  zy  3z 


3v  3v  3v 

Z  _ z  _ z 

+  m  — r—  +  m  —r~  +  m  — r— 

xz  3x  yz  3y  zz  3z 


3q  3q  3q 

+-£  +  ^  +  Ph 

3x  dy  3z 


We  also  note  that 


(18.11) 


'  3c  ,  3e  ,  3c 
e  =  —  +  —  v  +  —  v  +t-v 
3t  3x  x  dy  y  3z  z 


1 


Again  the  energy  balance  equation  (18.9)  is  exact  in  this  form  and  is 
valid  for  the  nonlinear  theory- 
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19.  PRINCIPLE  OF  ENTROPY 

For  certain  classes  of  physical  phenomena  within  a  range  of  expected 
changes,  a  material  body  is  characterized  by  certain  constitutive  equations. 
These  equations  define  an  ideal  material  approximating  the  real  material 
under  consideration.  For  any  thermomechanical  change,  the  constitutive 
equations  assumed  must  not  violate  the  second  law  of  thermodynamics  In 
continuum  mechanics,  this  law  may  be  stated  as  follows: 

V .  Principle  of  Entropy  (Clauslus-Duhem  Inequality) .  The  time  rate  of 

change  of  the  total  entropy  H  _is_  never  less  than  the  entropy  influx  through 
the  surface  S  o_f  the  body  and  the  volume  entropy  supply  B  _in  the  body .  This 
is  postulated  to  be  true  for  all  parts  of  the  body  and  for  all  independent 
processes  (Eringen  [1966c]). 

Accordingly  we  write 

(19.1)  T  =  ~  -  B  -  •  da  >_  0 

s 

where  r  so  defined  is  the  total  entropy  production.  For  simple  thermomechani¬ 
cal  processes  we  have 

(19.2)  H  =  / pndv 

1/ 

(19.3)  B  =  /  ”  dv 

1/  6 

(19.4)  S  5  | 

where  n,  h,  q  and  6  are,  respectively,  the  entropy  density,  heat  source,  heat 
vector,  and  the  absolute  temperature.  Substituting  (19.2)  to  (19.4)  into  (19.1) 


gives 
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(19.5)  r  =  4r  /pndv  -  j  ^  dv  -  ^  •—  da,  >_  0 

dt  \J  l r  S  6  k 

Using  the  Green-Gauss  theorem  to  convert  the  surface  integral  to  a 
volume  integral  and  carrying  out  the  differentiation  with  respect  to  time, 
we  get 

(19.6)  r  ;/|p^-(^),k-|]dv>0 

Since  this  is  to  be  valid  for  all  parts  of  the  body,  we  must  have 

(19.7)  (T>>k- 


This  is  the  well-known  Clausius-Duhem  inequality  of  classical  continuum 
mechanics.  In  micropolar  bodies,  this  is  considered  to  be  unchanged. 

Upon  substituting  h  solved  from  (18.9),  we  can  rearrange  (19.7) 
into  the  following  form 


(19.8) 


py  i  o(n  -  f)  +i  tMv,ik  -  £  ik„ntmnYk 


+  f  mlkvk,l  +~2  V-k 

U 


Still  another  form,  convenient  for  some  cases,  is  found  by  introducing  the 
Helmholtz  free  energy 


(19.9) 


0  *  t  -  On 


Hence 


CY  2  -  £  U>  +  rO)  +  f  t.  „v 


t  • 


1 


6 


t  •), 


0  ki  £,k  6  kmn  mn  k 


+  imikuk,i  ‘I’iS-1 

U 


(19.10) 
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The  Clausius-Duhem  inequality  (19.8)  or  (19.10)  is  postulated  to  be  valid 
for  all  independent  thermomechanical  changes.  This  implies  that  we  must 
know  the  independent  variables  which  affect  n,  t^, ,  m^,,  and  q^  at  the 
outset.  This  in  turn  requires  writing  constitutive  equations  for  these 
variables.  An  example  of  this  is  to  be  found  in  the  following  article. 

We  also  note  that  the  wntropy  inequality  (19.5)  can  be  shown  to 
lead  to  an  inequality  restricting  the  normal  component  of  q/0  on  the  surface 
of  the  body  (Eringen  [1966c] )  ,i  .  e. , 

q 

(19.11)  (tl  ‘  n  >  0  on  S 

~  V  -  '  — 

where  a  boldface  bracket  indicates  the  difference  on  the  quantity  enclosed 
calculated  from  two  different  sides  of  S,  on  S. 
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20.  THEORY  OF  MICROPOLAR  ELASTICITY 

A  micropolar  elastic  solid  is  distinguished  from  an  elastic  solid 
by  the  fact  that  it  can  support  body  and  surface  couples.  These  solids 
can  undergo  local  deformations  and  microrotations.  Such  materials  may 
be  imagined  as  bodies  which  are  made  of  rigid  shorL  cylinders  or 
dumbbell  type  molecules. 

From  a  continuum  mechanical  point  of  view,  micropolar  elastic  solids 
may  be  characterized  by  a  set  of  constitutive  equations  which  define  the 
elastic  properties  of  such  materials.  A  linear  theory  as  a  special  case  of 
the  nonlinear  theory  of  microelastic  solids  was  first  constructed  by  Eringen 
and  Suhubi  [1964a,b].  Later,  Eringen  [1965],  [1966]  reorganized  and  extended 
this  theory.  Here  we  give  a  self-contained  account  of  this  theory. 

In  linear  micropolar  elasticity,  the  strain  measures  are  (cf.  equations 
(4.17)  and  (4.18)) 


(20.1)  e  =  e  +  e  (r  -  tp  )  *  u0  ,  +  e  $ 

k£  k£  k£m  m  m  £,k  1km  m 


(20.2)  Y.  =  e,  $ 

k£m  k£n  n,m 


Since  only  the  nine  components  .  of  y.  are  independent  and  non-vanishing 

K  |  X  KXITl 

(cf.  equation  (4.37)),  we  may  instead  of  Y^^m  use  the  axial  tensor  4>^  ^  for 
simplicity.  Upon  arbitrary  rotations  and  reflections  of  the  spatial  coordinates 
represented  by  Q^,  i.e., 


(20.3) 


*k  •  2k! N 


transform  according  to 
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Lkf  ^kmCmn‘'£n 

(20.5) 

<t  =  tQ  <♦>  0 

k,£  i<m 

where  in  the  last  equation  the  plus  sign  is  for  det  Q  *  +1  and  the  minus 
sign  for  det  =  -1.  This  is  because  4^  is  an  axial  vector.  Equations 
(20.5)  express  the  fact  that  both  and  4^  ^  are  objective  tensors  and 
are  appropriate  for  use  as  independent  constitutive  variables  To  this 
list  of  variables  we  also  include  the  temperature  6  so  that  the  material 
properties  of  these  materials  may  depend  on  the  temperature  as  well.  The 
constitutive  dependent  variables  are 

‘kt*  v  * and " 


We  now  propose  a  set  of  constitutive  equations  of  the  form 


Ekt  "  Fkt(Ers’  *r,s’  6) 


"k£  '  VErs'  6> 


(20.6) 


“k  '  Gk<Ers  >  *r,S-  6) 


4  *  ,4  ,  0) 

rs  r,s 


n  ’  S<Ers  '  *r,s-  6) 


The  above  equations  are  legitimate  for  linear  homogeneous  materials  whether 
isotropic  or  not.  For  nonlinear  isotropic  materia  s  they  are  acceptable  in 
form.  However,  since  we  are  employing  the  infinitesimal  strain  measures,  a 
nonlinear  constitutive  theory  in  terms  of  linear  strain  measures  would  be 
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inconsistent .  For  the  nonlinear  theory,  the  reader  is  referred  to  Eringen 
and  Suhubi  [1964a  &  b ] . 

The  constitutive  equation  (20.6)  must  be  consistent  with  the  second 
law  of  thermodynamics  as  expressed  by  (19.10).  Thus,  upon  substituting 
(20.6)  into  (19.10),  we  have 


£  ( 

6  v3e 


k£ 


Ck£  + 


3Y 


3 


k ,  H 


■  S  ui  •  •  i  • 

pk,Z  +  7e  9  +  r'6)  +  I  tk«.Ck£ 


(20.7) 


+  i\^n,k  +  i2  qk6’k-° 

y 


Consistent  with  the  linear  theory  we  write 


(20.8) 


— — -  (A  )  =  d) 

Dt  v\,r  \ ,z 
Ck£  Vk,£  Eki.mVm 


The  inequality  (20.7)  is  postulated  to  be  valid  for  all  independent  processes. 
•  •  • 

Here  0  and  0,^  can  be  varied  independently.  Since  this  inequality 
is  linear  in  all  these  variables,  we  must  set  the  coefficients  of  these  varia¬ 
bles  equal  to  zero.  Hence 


(20.9) 


dV 


~ki 


3e 


k£ 

3Y 


\i 


=  o 


3 


Z,k 


qk  = 


0 


n  = 


21 

30 
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We  therefore  see  that  for  jj  mlcropolar  elastic  solid  stress,  the  couple  stress 
and  entropy  density  are  derivable  from  _a  potential  and  the  heat  vector  vanishes . 
Since  we  did  not  consider  the  temperature  gradient,  we  have  no  heat  conduction. 
Nevertheless,  the  free  energy  Y  and  consequently  the  material  moduli  will 
depend  on  the  temperature  0.  Since  all  terms  in  (20.7)  vanish,  we  have  the 
entropy  production  density  y  also  vanishing. Thus . the  micropolar  elastic  solid 
is  in  thermal  equilibrium. 

Here  we  are  concerned  with  the  linear  theory.  We  therefore  consider  a 
polynomial  for  4>  which  is  second  degree  in  the  strain  measures  and  4^  i.e., 

k£  k£  2  x2mn  mn  k£  k,£  2  k£mnyk,£  m,n 

(20.10) 

■f  C  c.  0 

k£mn  ki  m,n 

where  Ao,  A,  ,  A,  ,  B  .  ,  ...  are  functions  of  0  only.  Since  4>  is  an  axial 
k  kiomn  k*.  fourth  and  the  k 

vector,  upon  a  reflection  of  the  spatial  axes  theAlast  termswill  change  sign 
while  Lhe  other  terms  do  not.  For  the  function  Y  to  be  invariant  B^  =  0, 

^k£mn  *  We  furCher  note  the  following  symmetry  conditions  which  are  clear 
from  various  summations  in  (20.10) 


(20.11) 


^kUmn 


mnkll 


B 


k£mn 


mnk£ 


which  shows  that  for  the  mos t  general  mir rppolar  anisotropic  elas tic  solid. 

the  number  of  A,  „  and  B  ,  „  is  45  each.  In  addition,  we  have  nine  A,  „ 
- k£mn  -  mnk£  —  -  K.S. 

which  give  rise  to  an  Initial  stress  in  the  undeformed  state  of  the  body. 
Upon  substituting  (20.10)  into  (20. 9)^  and  (20  9)^,  we  obtain 

(20.12)  t.  ■  A.  +  A  ,  s. 

ki  k£  Ukrnn  mn 


m,n 


(20.13) 


^£kmn* 
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These  are  the  linear  forms  of  the  stress  and  couple  stress  constitutive 
equations  for  the  anisotropic  micropolar  elastic  solids.  When  the  initial 
stress  is  zero,  then  we  must  also  have  *  0.  Thus,  for  the  micropolar 
solid  free  of  initial  stress  and  couple  stress,  we  have 


(20.14) 


'k£ 


mn£mn 


(20.15) 


"W  ^£kmn^m,n 


Various  material  symmetry  conditions  place  further  restrictions  on  the  con¬ 
stitutive  coefficients  A„,  and  B„.  .  These  restrictions  are  found  in 

£kmn  £kmn 

the  same  manner  as  in  classical  elasticity.  Here  we  obtain  the  case  of 
isotropic  solids.  If  the  body  is  isotropic  with  respect  to  both  the  stress 
and  couple  stress,  we  call  the  solid  microisotropic ■  In  this  case,  the 
constitutive  coefficients  must  be  isotropic  tensors.  For  the  second  and 
fourth-order  isotropic  tensors,  we  have  the  most  general  forms 

Tc£  ki.  \£mn  1  k£  mn  2  km  £n  3  kn  £m 

(20.16) 

B.  „  =  B.6,  .  6  +  B-6.  +  B_6,  6„ 

k£mn  1  ki  mn  2  km  in  3  kn  £m 


where  A,  A^  ,  A^ ,  A^,  B^,  B2,  and  B^  are  functions  of  9  only  In  this  case 
then,  (20.12)  and  (20.13)  take  the  special  forms 


<2° . 17)  tk(  -  ASkJ  *  Akr  rr6k|,  +  A^,  +  A^ 


(20.18) 


For  vanishing  initial  stress  A  =  0.  Introducing 
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(20.19) 


X  ,  A2  i  w  +  <  ,  A^  =  v 

B1  =  a  ,  B2  =  y  ,  B3  =  B 


the  above  equations  can  be  written  as 

(20 . 20)  tkl  •  »‘rr«kl  +  (u  +  E>Ekl  +  UE*k 


(20.21) 


\£  "  a*rtr\l  +  6*k,£  +  A,k 


For  the  free  energy  in  this  case  we  find 


(20.22) 


pl'-  2  <XEkkE«  +  (U+  <)Ek«Ek£+  “Vjk1 


T  (“*k,k*t,l  +  6*k,t*t,k  + 


An  alternative  form  to  (20.20)  to  (20.22)  is 


(20.23) 

Ck£  “ 

Ae  ' 

rr  k£ 

(20.24) 

%£  " 

“V.At 

(20.25) 

pY  = 

2  1 lekke 

-kit  kim  m  V 


"MM- 


Lk  Tk/V‘k  Yk' 


+  2  ia\k\l  +  B*k,A,k  + 


We  note  the  difference  between  isotropic  micropolar  elasticity  and  classical 
elasticity  by  the  presence  of  four  extra  elastic  moduli,  namely,  <,  a,  B  and 
Y.  When  these  are  set  equal  to  zero,  the  above  equations  (20.23)  to  (20.25) 
revert  to  Hooke's  law  of  the  linear  isotropic  elastic  solid. 
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21.  RESTRICTIONS  ON  MICROPOLAR  ELASTIC  MODULI 


The  stability  of  materials  requires  that  the  stored  elastic  energy 
be  nonnegative.  This  condition  is  also  essential  for  the  uniqueness  of 
the  solutions.  This  requirement  places  certain  restrictions  on  the  micro- 
polar  elastic  moduli.  The  following  theorem,  Eringen  [1966a ], provides  these 
conditions  for  <p  independent  of  6. 

Theorem.  The  necessary  and  sufficient  conditions  for  the  internal 
energy  to  be  nonnegative  are 

0^3A+2y+K,  0  <_  u  ,  0^_k 

(21.1) 

0  <_  3a  +  2y  ,  -y  ^  8  < ,  0  y 

The  sufficiency  of  (21.1)  is  proven  by  observing  that  when  these  inequalities 
hold, each  one  of  the  following  energies  constituting  the  internal  energy 
density  is  nonnegative 


(21.2) 


e  *  c  +  e  +  e 
E  R  M 


where 


PC-T 


1 


2  [lekk*«  +  <2“  +  ^‘kl'lk1 


(21.3) 


oer  5  <(rk  -  V(rk  -  V 


prM  ■  I  (c,*k,k‘i,t  +  6‘k,t*t,k  +  Y\.k‘e,k> 


The  fact  that  pe„  is  nonnegative  under  the  conditions  (21.1),  and  (21.1)- 
t  12 

is  well  known  for  the  classical  elasticity.  It  is  simple  to  observe  that 
peR  is  nonnegative  for  rR  4  whenever  <  _>_0.  To  see  the  same  for  pc^. 


we  write  this  expression  as 


-101- 


(21.4) 


where 


2(,£„  ■  t <3°  +  B  +  ,)tk.k*M +  (y  -  6)  vuV.m 


+  (Y  +  8)!*(k,l)  '  3  ',r,r\t',',(k,H)  "  3  is.s'kl1 


(21.5) 


*(k,l)  2  Uk,£  +  \,k} 


’[k.H]  1  (\,£  “  4>e,k) 


From  (21.4)  it  is  clear  that  when 


(21.6)  3a+B  +  Y2_0  *  y-6^_0  ,  v  +  B  2.  0 

we  have  pew  >  0  so  that  (21.1)  are  sufficient  for  pew  >  0 

M  ~  M  — 

Conditions  (21.1)  are  also  necessary  for  the  nonnegativeness  of  pe. 

To  prove  this,  we  recall  that  e^,  -  4^  and  4^  ^  can  be  varied  independ¬ 

ently  of  each  other.  Since  the  above  three  energies  are  uncoupled  with 
respect  to  these  variables,  each  one  of  these  energies  must  be  nonnegative 
independent  of  each  other.  The  fact  that  (21. 1)^  and  (21. 1)^  are  necessary 

for  pe  >  0  is  known  to  us  from  classical  elasticity.  Excluding  the  case  of 
t 

r^  “  4^  (Indeterminate  couple  s  tress  theory),  we  see  that  pe^  is  nonnegative 
if  and  only  if  <  >_  0.  For  the  case  of  r^  *  4^,  by  replacing  2„  +  <  by  a 
new  modulus  2 y  we  shall  have  <  disappear  from  the  constitutive  equations. 
Thus,  it  remains  to  prove  the  necessity  of  (21.1)  for  the  nonnegativeness 
of  pc^.  To  this  we  write  it  as  a  quadratic  form  in  a  nine-dimensional  space, 
i.e.  , 

(21. 7)  P£m  -  y±y  ,  a  -  a  ,  (i,j  =  1,2,  ,9) 


where 
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yl  '  *1  1  ’  ^  *2  2  ’  y3 

yA  :  2  *  ¥ 5  :  ^2  1  *  y6 

(21.8)  y 7  5  ^3,2  ’  y 8  *3,1  ’  y9 

all  "  a22  "  a33  "  0  +  6  +  Y*  a44  “  a55  " 

a12  “  a13  "  a23  "  °  ’  °45  "  a67 

all  other  ■  0 

The  characteristic  values  a^  of  a^  are  obtained 

(21.9)  det  (a  -  a6  )  «  0 

i  j  i  j 


=  4> 


3,3 


=  4> 


2,3 


=  <t> 


1,3 


a66  =  a77  “  a88  “  a99 


a89  "  6 


by  solving  the  equation 


Y 


The  nine  roots  a^  of 
al  -  a2  -  a3  -  Y 

a^  =  3a  +  8  +  y 

In  order  for  pe,.  >  0 
M  — 

cient  that)  (21.7)  be 


this  equation  are 

-  8  f  -  a5  =  a6  -  a?  =  ag  -  y  +  8 


to  be  satisfied  for  all  y  ,  it  is  necessary  (and 
an  ellipsoid  in  nine-dimensional  space,  i.e., 


suf fi- 


Y  -  8  >_  0  ,  y+B^O  ,  3a  +  8  +  >  >_  0 

This  set  of  conditions  is  the  same  as  the  last  three  conditions  of  (21.1). 
Hence  the  proof  of  the  theorem. 

The  nonnegative  character  of  the  internal  energy  density  has  important 
implications  in  regard  to  uniqueness  theorems  in  both  static  and  dynamic  micro 
polar  elasticity.  For  these  and  other  important  results  see  Eringen  [1966a, bJ 
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22.  FIELD  EQUATIONS,  BOUNDARY  AND  INITIAL  CONDITIONS 

The  field  equations  of  linear  micropolar  elasticity  are  obtained  by 
substituting  (20.23)  and  (20.24)  into  (17.17)  and  (17.18).  Hence 


(22.1)  (X  ♦  u)uI  tk  +  <u  +  <)uk>u  +  +  P(fk  -  uk)  -  0 


(22.2)  (a  +  6)«t>lk  +  Y*k>u  +  -  2,.,  +  p»k  -  J»k>  -  0 


where  we  have  taken  ■  j 6^  for  the  microisotropic  solid.  In  the  linear 
theory,  p  and  j  are  considered  constants  and  the  accelerations  and 

are  calculated  by  their  approximate  expressions 

.2 

C 

(22.3) 


at 


2 


324> 


♦k" 


3 1 


The  vectorial  forms  of  these  equations  are  found  to  be  convenient  for  the 
treatment  of  problems  in  curvilinear  coordinates.  These  are  readily  obtained 
from  the  above  equations  by  simply  multiplying  them  by  i^  and  observing  that 


(22.4) 


U£,£k*k 


VV 


EkJ-m^m, 


V  x  $ 


uk,U*k 


VV 


V  X  V  X  u 


where  V  is  the  gradient  operator  so  that 

V 4>  =  grad  <p  ,  V  •  u  =  div  u  ,  V  x  u  =  curl  u 


Hence 

(22.5)  (X  +  2p  +  k)VV  •  u  -  (u  +  <)V  x  y  x  u  +  <v  x  $  +  p(f  -  U)  *  0 

(22.6)  (a  +  8  +  y)VV  •  4>-yVxVx$  +  <Vxu-  2<<p  +  p(£  -  j$)  -  0 


-104- 


For  an  initial 

value  problem,  the 

initial  conditions  have  the  form 

(22.7) 

u(x,0)  -  u0(x) 

• 

u(x,0)  =  Vg(x) 

in  y 

(22.8) 

4>(x, 0)  -  $g (x) 

» 

<p  (  X ,0)  =  Vg  (x) 

where  ug,  vg,  0g,  and  \>q  are  prescribed  in  M  at  time  t  =  0. 

Many  different  types  of  boundary  conditions  are  suggested  in  applica¬ 
tions.  For  example,  we  may  prescribe 

(22.9)  u(x',t)  =  u’ 

x'  on  S 

(22.10)  <Kx\t)  = 

on  the  boundary  surface  S  of  the  body.  An  equally  permissible  set  of  boundary 
conditions  requires  the  prescription  of  the  tractions  and  couples,  i.e., 


(22.11) 

(22.12) 


£kn  £  =  C(n)k 

on  S 

m£k'V  =  m(n) k 


where  t^  and  m^  are  the  stress  and  the  couple  stress  tensors  given  by  (20.23) 

and  (20.24)  and  t,  and  m,  are  the  prescribed  tractions  and  couples  on 
(n)k  (n)k 

S  whose  exterior  normal  is  n. 

In  some  other  problems,  a  mixture  of  the  above  two  types  of  conditions 
occurs,  e.g.,  on  some  part  of  S  one  may  have  (22.9)  and  (22.10)  and  on  the 
remainder  S ^  ±  S  -  the  conditions  (22.11)  and  (22.12).  Still  other  types 
of  mixed  conditions  involving  some  components  of  one  set  and  the  remaining 
components  of  the  other  set  are  possible.  All  admissible  sets  of  boundary  con¬ 
ditions  allowing  unique  solutions  must  satisfy  (Eringen  (1966a]) 
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<22-13)  E(?)k\ +  >)k\  • 0  °nS  •  ‘i" 

where  u,  #,  t^,  and  aie,  respectively,  the  difference  of  u,  $,  t  ^  ^ 

and  from  their  respective  values  on  S. 

The  field  equations  (22.1)  and  (22.2)  are  valid  only  for  micropolar 
isotropic  solids.  Note  that  for  vanishing  k,  a,  6,  £,  and  j,  equation  (22.2) 
reduces  to  0  *  0,  and  (22.1)  gives  the  celebrated  equations  of  Navier  of 
classical  elasticity. 

For  the  anisotropic  micropolar  elastic  solid,  the  field  equations 
replacing  (22.1)  and  (22.2)  are  obtained  by  substituting  (20.14)  and  (20.15) 
into  (17.17)  and  (17.18) 


(22.14)  A  .  (u  „  +  c  *  „)  +  p  (f.  -  IL  )  -  0 

£kmn  n,m£  nmr  r,£  k  He 


(22.15)  B.  .  <|>  „  +  e.  A  (u  +  e  <t>  )  +  p(£.  -  ■).,<$„)  -  0 

k£mn  m,n£  kmn  mnpq  q,p  qpr  r  k  J£k  £ 


For  the  expressions  of  and  we  used  (20.1)  and  (15.18). 

Finally,  we  record  below  the  expanded  forms  of  the  field  equations 

(22.1)  and  (22.2)  in  rectangular  coordinates. 

~  9u  9u  9u  9^u  92u  9^u 

, ,  ,  \  9  x  ,  y  ,  z .  ,  ,  .  .  , _ x  ,  _ x  ,  _ xx 

(A  +  U)  97  (“+  "ly  +  ~ +  (U  +  <){~~2+~Z  +  ~l) 

9x  9y  9z 


9$  d<P  9  u 

+  +  0<£x  - 

,2  ,2  ,2 

9  u  9  u  9  u 


)  =  0 


~  9u  9u  9u 

(A  +  M)  T7  (“^  +  +  "lf)  +  (p  +  <)(~2  +  “~1>  +  mTz) 

J  9x  9y  9z 


(22.16) 


9  9$  92u 

+  <(— - rr-)  +  o(f - 

9z  9x  y  3  2 


■)  -  0 


2  2  2 

.  9u  9u  9u  9  u  9  u  9  u 

<i  +  “)17<lf  +  -#+-3f>  +  (“  +  '><— f  +  -i+— f> 

9x  9y  9z 


d<p  __  9 (f) 

+  <(■ 


92u 


S  -  ^7> +  0(£z  -  ~i>  ■  ° 
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2  2  2 

dtp  dtp  dtp  d  <P  d  p  d  4> 

<0  +  6)  u  (~+ 17+  ~,+Y(_77  +  “T +  — ) +  K(— 

3x  9y  3z 

-  2<<i  +  p(i.  -  j  - h  -  0 

x  x  dz 


9u  9u 


9z; 


2  2 

3  $  9'  p 


d2tp 


dtp  dp  dp 

<“  +  s)  TJ  <_af  +  +  — 2  +  — 2  +  ~ £>  *  K(lt 

dx  dy  3z 


3u  3u 


9x' 


(22.17) 


-  2**„  +  p<£  -  j 


d2p 


r)  -  o 


9t 


2 

9  4>  9  $  dp 

(»  +  6)  £  +  -r!)+Y( - T  +  - 1  +  —I)  +  <( 


3  4>  9$  dp 


3u  3u 


9z  9x  9y  9z 


3x 


9y 


3z 


9x  9y 


d2p 

-  2<p  +  p(i  -  j  - &  -  0 

2  2  9t 
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23.  INDETERMINATE  COUPLE  STRESS  THEORY 


In  Arc.  22  we  Indicated  that  classical  elasticity  is  a  special  case 
of  the  micropolar  theory  when  <*»a“S“Y*‘jM0  and  £  ■  0.  The  same 
is  also  true  for  motions  in  which  the  microrotation  vanishes.  There  exist 
other  classes  of  constrained  motions  which  have  attracted  the  attention  of 
research  workers.  The  most  popular  among  them  is  the  indetermine  couple 
stress  theory  which  is  contained  in  the  work  of  the  Cosserats  [1909].  Recently, 
Truesdell  and  Toupin  [1960],  Grioli  [1960],  Aero  and  Kuvshinskii  [1960],  Mindlin 
and  Tiersten  [1962],  Toupin  [1962],  and  Eringen  [1962]  independently  presented 
new  derivations  and  supplied  various  missing  parts  qf  the  theory.  This  theory 
can  be  obtained  as  a  special  case  of  the  micropolar  theory  if  the  constraints 


(23.1) 


A  «  r  «  —  r  II 

X  rk  2  klm  m,  l 


are  imposed.  In  this  case,  the  stress  constitutive  equations  (20.23)  reduce 
to 


(23,2)  C(k£)  “  Xerr6k£  +  (2u  +  <)ek£ 

where  a  parenthesis  enclosing  indices  as  usual  indicates  the  symmetric  part 
of  the  stress  tensor.  We  also  use  a  bracket  to  denote  the  antisymmetric  part 
of  tensors ,  e. g.  , 

a(k£)  =  2  (ak£  +  a£k^  *  a[k£]  =  7  (ak£  "  a£k) 

Thus,  when  (23.1)  is  valid,  the  antisymmetric  part  of  the  stress  disappears 
from  the  constitutive  equations.  We  can,  however,  remedy  this  situation  by 
another  artifice.  The  equations  of  moment  of  momentum  (17.18)  can  be  used 
to  solve  for  the  antisymmetric  part  of  the  stress  tensor.  Multiplying  (17.18) 
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by  Lkrs '  we  so^ve  ^or 


(23.3) 


1  1  /4  '  V 

-  r  p£  ,  ,  l  -  o  ) 

[ki]  2  rki  nr,n  2  rki  r  r 


where  we  used  the  identity 


(23.4) 


g  .  £  ■  6.  6.  -  6,  6. 

rki  rmn  km  in  kn  im 


Upon  carrying  (23.1)  into  (20.24)  we  have 


(23.5) 


"hci  2  £krsus,ri  +  2  eirsUs,rk 


Carrying  (23.5)  and  the  expression  of 


(23.6) 


j  -  j<*>  -  f  je  U 

r  J  r  2  J  rmn  n,m 


into  (23.3),  we  have 


(23.7) 


C[ki]  "  2  V  u[k, i]  "  2  p(£rki£r  +  Ju[k,i]) 


where  V  is  the  laplacian  operator  in  rectangular  coordinates 


(23.8) 


7  ut  5  ut.kk 


If  we  not  substitute  (23.2)  and  (23.7)  into 


(23.9) 


Ck£  "  C(k£)  +  C[ki] 


we  get  the  total  stress  tensor 


Ckt  '  iur,r4kl  +  <u  +  2)(uk,l  +  “l.k* 


+  2  v  u[k,l]  "  2  0<erkltr  +  JU[k,il) 


(23.10) 
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The  presence  of  the  body  couple  i  and  the  acceleration  in  this  equation 
is  certainly  disturbing  since  the  constitutive  equations,  in  general,  should 
not  contain  such  terms.  When  (23.9)  is  carried  into  the  equation  of  balance 
of  momentum  (17.17),  we  find 

(23.11)  (1  +  M  +  f  +  f  '2)\>lk  +  +  f  -  i  ’2)'‘l>kk 

*  2  £klrlr.t  +  o£l  ■  0(1  -  4  'X  +  ¥  “k,tk 

By  use  of  identities  (22.4),  we  may  also  obtain  the  vector  form  of 
these  equations: 

(23.12)  (A  +  2y  +  <)?V  •  u  -  (W  +  J  “  J  V2)?  *  7  x  y 

+  p(f  +  j  V  x  i)  ..  P(i  +  1  v  x  Vx)u  -  0 

where  for  the  lapalcian  operator  we  have  used 

(23.13)  V24  »  W  •  A  -  7  x  ?  x  a 

Equation  (23.12)  takes  the  form  (3.27)  obtained  by  Mindlin  and  Tiersten  [1962] 
in  an  entirely  different  way  if  we  write  u  for  y  +  </2  and  n  for  y/4  and  j  *  0. 
Thus,  these  author^  as  with  others,  have  neglected  the  micropolar  rotatory 
inertia.  Equations  (23.11)  or  (23.12)  are  the  field  equations  of  the  theory 
known  as  the  (indeterminate)  couple  stress  theory.  It  is  to  be  observed  that 
in  this  theory  the  skew- symmetric  part  of  stress  and  consequently  the  stress 
are  dependent  on  the  applied  loads  and  inertia,  and  they  are  not  determined 
solely  as  a  result  of  the  constitutive  character  of  the  medium.  This  violates 
the  axiom  objectivity  since  the  applied  loads  and  inertia  terms  involved  are 
not  objective  quantities.  A  second  relevant  point  is  that  while  according  to 
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(23.5)  one  can  determine  both  the  symmetric  and  the  skew-symmetric  parts 
of  the  couple  stress,  this  in  fact  is  not  the  case  of  the  couple  stress 
theory  in  which  constitutive  equations  are  derived  from  a  free  energy 
function,  as  is  done  in  Art.  20.  It  is  clear  from  (20.9)2>  for  example, 
that  all  components  of  ^  can  no  longer  be  used  as  independent  variables. 
In  fact,  if  one  uses  (23.1)  in  the  argument  of  free  energy  t,  one  finds 


symmetric  parts  of  the  stress  and  couple  stress  remain  indeterminate  . 

This  is  the  reason  for  the  use  of  the  terminology  "indeterminate". 

This  situation  has  certain  similarity  to  the  isochoric  motions  of 
compressible  bodies  as  compared  to  the  motions  of  incompressible  solids. 

As  is  well-known,  in  the  latter  case  the  pressure  is  not  determined  through 
the  constitutive  equations.  Finally,  in  the  indeterminate  couple  stress 
theory,  the  number  of  boundary  conditions  on  the  surface  tractions  and 
couples  must  be  reduced  from  six  to  five.  A  consistent  set  of  boundary 
conditions  must  not  violate  the  uniqueness  theorem.  Mindlin  and  Tiersten 
[1962]  have  obtained  a  uniqueness  theorem  for  the  following  set  of  boundary 
conditions . 

Let  (x^,  x2>  x^)  be  a  set  of  orthogonal  curvilinear  coordinates 
taken  in  such  a  way  thac  x^  *  :<3  locally  coincides  with  the  surface  S  of 
the  body.  The  boundary  conditions  consist  of  specifying  at  x^  *  x“  one 
factor  in  each  of  the  five  products 

(23. 1A)  Piui  ,  P2u2,  t^33ju3  ,  m31r1  ,  m32r2 


1 


In  this  regard,  see  the  discussion  given  in  Eringen  [1962,  Art.  40],  See 
also  Mindlin  and  Tiersten  [1962],  Toupin  [1962],  and  Eringen  [1964b]. 
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where 


(23.15) 


Pi  =  *(31)  "  2  (m12;i  “  m33;2* 
^2  =  *(32)  +  2  (mil;i  ‘  m33;l) 


Here  an  index  placed  after  a  semicolon  denotes  directional  differentiation 

along  the  corresponding  curvilinear  coordinate,  and  a  superposed  bar  the 

boundary  values  of  the  quantities  involved.  Curvilinear  components  of  the 

displacement  vector  on  S  are  denoted  by  u^,  the  couple  stress  by  m^,  and 

the  stress  tensor  by  t,  . . 

kH 

If  an  edge  is  an  intersection  of  two  orthogonal  surfaces  ■  x° 
and  ■  x^,  then  we  must  also  specify 


(23.16) 


[m  ]  -  [m  }  or  u 

33  -0  11  -o  2 

3  1 


The  reduction  of  the  number  of  boundary  conditions  from  six  to  the  above 
five,  (23.15),  is  similar  to  the  one  encountered  in  the  Bernoulli-Euler 
theory  of  thin  plates.  Conditions  (23.16)  are  the  analogs  to  the  corner 
conditions  . 

Under  suitable  regularity  assumptions,  the  above  five  boundary  condi¬ 
tions,  together  with  the  assignment  of  the  body  force  field  pf,  curl  (pH), 
and  the  initial  values  of  u  and  u  (with  j  =  0),  are  sufficient  for  the  unique 
determination  of  e^,  and  ^*e  displacement  field  u  and 

the  rotation  r  are  unique  too  within  an  arbitrary  rigid  body  displacement 
field. 

The  indeterminate  couple  stress  theory  described  above  has  many 
apparent  limitations.  Whether  or  not  bodies  can  undergo  such  constrained 
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motions  is  not  known.  The  appearance  of  nonobjective  quantities  in  the 
constitutive  equations,  the  limitations  on  the  spin  inertia  and  the  body 
couple  field,  and  physically  unnatural  boundary  conditions  leaves  much 
to  be  desired.  Experimental  comparison  for  certain  practical  applications 
exists  (cf.  Schijve  [1966]).  A  discussion  of  the  weakness  of  this  theory 
indicating  the  disagreement  of  the  theoretical  results  with  the  experiments 
has  been  given  by  Kaloni  and  Ariman  [1967).  Nevertheless,  recent  litera¬ 
ture  contains  a  large  number  of  solutions  in  the  field. 


2 
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24.  PROPAGATION  OF  WAVES  IN  AN  INFINITE 
MICROPOLAR  ELASTIC  SOLID1 

Here  and  in  the  following  several  articles  we  investigate  solutions 
of  certain  dynamic  and  static  problems  in  linear  isotropic  micropolar  elas¬ 
ticity.  Essential  to  these  problems  are  the  field  equations  (22.5)  and 
(22.6)  and  boundary  conditions  of  the  type  listed  in  (22.11)  to  (22.12). 

The  propagation  of  linear  isotropic  micropolar  elastic  waves  with 
vanishing  body  loads  is  governed  by  (22.5)  and  (22.6)  or 


(24.1)  (c^  +  c^)VV  •  u  -  (c*  +  c^)V  x  y  x  y  +  C^V  *  $  -  (i 

(24.2)  (c^  +  c^)VV  •  £  -  c^V  *  V  x  x  u  -  2iiK(J>  = 


4  b 

„  l 

0 

01 

i  nere 

2 

1+2  u 

2 

-  J± 

2 

(24.3) 

C1  = 

P 

’  C2 

P  ’ 

C3 

P 

2 

2  . 
C4  = 

JL 

Pj 

2 

*  C5 

III 

oft 

CM  O 

3 

_  C3 

=  j 

< 

Pj 

We  decompose  the  vectors 

u  and 

into 

scalar 

and 

vector 

potentials  as 

follows : 


(24.4) 


u  «  Vu  +  V  ‘  D  ,  V  •  U  ■*  0 

b  =  V$  +  V  *  0  ,  V  •  $  *  0 


Substituting  these  into  (24.1)  and  (24.2),  we  see  that  these  equations  are 
satisfied  if 


1 


The  present  section  is  based  on  the  work  of  Parfitt  and  Eringen  [1966] 
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(24.5)  (c2  +  c2)V2u  -  U 

(24.6)  (c2  +  c2)V2<f>  -  2to2<(i  -  $ 

(24.7)  (c2  +  c2)V2U  +  c2V  x  «  =  U 

(24.8)  c272$  -  2uj2<t  +  to2V  x  u  -  $ 

4  O'  O' 

It  may  be  observed  that  (24.5)  and  (24.6)  are  uncoupled  for  the  scalar 
potentials  u  and  while  equations  (24.7)  and  (24.8)  for  the  vector 
potentials  are  coupled. 

Plane  waves  advancing  in  the  positive  direction  of  the  unit  vector 
n  may  be  expressed  as 

(24.9)  (u,  <(>,  U,  <H  =  {a,  b,  A,  B}  exp  [ik  (n  *  r  -  vt)] 

where  (a,b)  are  complex  constants,  (A,  B)  are  complex  constant  vectors,  k 
is  the  wave  number,  and  r  is  the  position  vector.  Thus 

(24.1°)  k  =  ~~  r  =  xkiR 

in  which  £  is  the  wave  length  and  i^  are  the  unit  rectangular  base  vectors. 
Substituting  (24.9)  into  (24,5)  gives 

(24.11)  v2  =  c2  +  c2  ■  (X  +  2y  +  <)/p 

which  shows  that  a  plane  wave  with  the  displacement  vector 

(24.12)  u^  -  ik^an  exp  [ik^(n  •  r  -  v^t) ] 

may  exist  in  the  direction  of  propation  n,  Fig.  24.1.  These  waves  are  the 
counterpart  of  the  classical  irrotational  waves  and  reduce  to  them  when  <  =  0. 


-115- 


We  designate  these  waves,  longitudinal  displacement  waves. 

A  second  scalar  plane  wave  is  the  solution  of  (24.6)  in  the  form 
(24.9).  The  wave  speed  in  this  case  is 

2  2  2  2  -2 

(24.13)  v^  =»  cf  +  c^  +  2oj“k 

2  4  5  0 

If  we  introduce  the  angular  frequency  u 

(24.14)  w  -  2tt f ^  =  2iTV^/i.  ■  kv^ 

ori-B+y _ 

Pj(l  -  -^r) 

PJW2 

The  speed  of  these  waves  depends  on  the  frequency.  Hence  they  are  disper¬ 
sive.  Since 

(24.16)  a  +  8  +  y  ^0 

we  see  that  such  waves  can  exist  whenever 

(24.17)  u2  >  »/2u)0 

These  waves  will  be  called  longitudinal  microrotation  waves.  Fig.  24.1. 

The  microrotation  vector  is  given  by 

(24.18)  $  *  V<fr  =  ik^bn  exp  [ik2(n  •  r  -  v^t)] 

For  the  case  of  w„  =  *^2uj  ui  ,  we  have  v„  »  00  and  the  wave  does  not  exist. 
2  0  c  2 

When  <^2  <  ,  v£  becomes  purely  imaginary,  i.e., 


the  wave  speed  may  be  expressed  as 

2  _1 
2w 

(24.15)  V2  -  (c l  +  c^l - j)  = 
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(24.19)  v2  -  »l|v2  ,  i  /T 

It  may  be  seen  that  a  standing  wave  of  the  form 

(24.20)  *  b  exp  (-  -|  n  •  r)  exp  (-iu2t) 


is  possible. 

For  such  waves,  propagation  is  possible  if  w2  >  «/2u> •  Hence, 

=  is  a  cut-of f  frequency  for  these  waves. 

The  vector  wave  solutions  are  obtained  by  substituting  (24.9)  into 
(24.7)  and  (24.8).  This  results  in  two  simultaneous  vector  equations  for 
the  unknowns  A  and  B. 


(24.21) 


a  A  +  io  n  *  B  =  0 

A'  D ' 

18 .n  x  a  +  8„B  =  0 

A'  '  B' 


which  for  nonvanishing 

(24.22)  n 
resulting  from  (24 . 4 ) 2 

“a 

(24.23) 


c*a ,  aB»^»  and  Sg  are  subject  to 

•  A  =  0  ,  n  •  B  ■  0 

and  (24.4)  .  Here 
4 

.  2  ,  2  2  2,  .  .  2 

k  (v  -  c2  -  c3)  ,  r  kc3 

2  2  2  2  2  -2 
0  .  eB  £  Kv  -  C l  -  2«Jk  z) 


Equations  (24.22)  show  that  the  vectors  A  and  B  lie  in  a  common  plane  whose 
unit  normal  is  n.  Solving  from  (24.21)2  for  B  we  have 


_i  T~  ? 

6B 


x  A 


(24.24) 
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Hen  ce  ,  the  three  vectors  n,  A,  and  B  are  mutually  perpendicular.  Moreover, 
vanishing  A  implies  vanishing  B.  These  two  types  of  waves  are,  therefore, 
coupled  and  cannot  exist  without  each  other.  From  (24. 4)^  and  (24. 4)^  we 
see  that  u  and  corresponding  to  U  and  4>  are  normal  to  each  other  and  to 
the  direction  of  propagation  n.  Hence  these  waves  are  transverse  waves. 

We  call  the  waves  that  are  associated  with  U.  transverse  displacement  waves . 
and  those  that  are  associated  with  <t> ,  transverse  microrotation  waves.  Fig. 
24.2.  The  transverse  displacement  waves  have  their  classical  analogues  in 
equivoluminal  waves  and  in_the  limit  they  reduce  to  these  waves. 

The  velocities  of  propagation  of  these  waves  are  determined  by 
carrying  (24.24)  into  (24.21)2  and  using  (24.22)^.  This  gives 

(24.25)  av^  +  bv^  +  c  ■  0 


where 


a  = 


2  -2 

1  -  2i/w 
0 


(24.26) 


22  2-2  2  2-2 
b  =  —  [ c^  +  c2(l  -  2u)qu)  )  +  c3(l  -  u)Qu  )] 


-  2.2  2. 

c  =  c4(c2  +  c3) 


The  positive  real  roots  of  (24.25)  are 


V3  "  1 2a 


(-b  +  /(>2-4ac  ) 


1/2 


(24.27) 


v  **  r  — —  ( 
4  l2a  ' 


-L  (_b  _  tfl 


1/2 


4ac)  ] 


By  studying  the  discriminant 
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/T~.  /r  2  2  2  .  2  .  1  2.  2  -2/ 

n»  -4ac  -  Uc^,  -  c2  -  +  2(c2  +  ~  cyww  ] 

2  2  2  -2 
+  4c3c4w0a) 

under  the  conditions  k  >_0,  y  ^0  compatible  with  (21.1),  we  find  that 

is  real  when  w  >  oa  and  v.  is  real  for  all  values  of  w.  The  frequency 
c  4 

i  /Tw0  is  again  a  cut-off  frequency  for  waves  propagating  with  velocity 

V 

In  summary,  we  found  that  in  an  infinite  micropolar  elastic  solid 
six  different  types  of  plane  waves  propagating  with  four  distinct  speeds 
of  propagation  can  exist.  These  are: 

(a)  A  longitudinal  displacement  wave  propagating  with  speed  v^. 

(b)  A  longitudinal  microrotation  wave  with  speed  v^  propagating 

in  the  longitudinal  direction  whenever  the  frequency  of  these  waves  is 
above  the  cut-off  frequency  These  two  types  of  waves  are  uncoupled. 

(c)  Two  sets  of  coupled  transverse  displacement  waves  and  transverse 
microtation  waves  at  speeds  v^  and  v^ .  Of  these,  the  waves  having  velocity 

v^  can  exist  when  10  >  otherwise  they  degenerate  into  distance  decaying  sinu¬ 
soidal  vibrations. 

A  detailed  analysis  of  the  wave  speeds  V2,  v^,  and  v^  is  sketched  in 
Figs.  (24.3)  and  (24.4). 

By  use  of  (21.1),  Parfitt  and  tringen  [1966]  have  shown  that  a  con¬ 
sistent  solution  for  v^  and  v^  requires  the  existence  of 

2  2  ,  2 
r  >  c  +  c 

4-2  c3 

or  the  additional  inequality 

^  >  y  +  < 


(24.28) 
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In  addition,  the  wave  speeds  as  a  function  of  u  must  satisfy 


(24.29) 

v4  (0) 

<  V^(oj  ) 

—  4c 

2,  v 

1  v4 (") 

(24.30) 

2  .  . 
v2(“) 

2/  v 
>  v3(“) 

8 

cni  <r 
> 

A 

Also  they  found  that 

(24.31) 

2  2 
V1  *  V4 

The  study  of  the  relative  magnitudes  of  v„,  v  ,  and  v  requires 

sum  of  the^  J  4 

knowledge  of  the  relative  magnitude  of  the^constitutive  coefficients 
a  +  6  with  respect  to  j*.  Thus 


(24.32) 


2  2  2 
V2  >  V3  >  V4 


2  2  2 
V  >  v  >  V 

3  -  2  -  4 


2  2  2 
V2  -  V3  -  V4 


for  w  <  a) 
c 


•k 

for  u>  <  u>  <  u 


k 

for  w  <  u> 


a  +  B  >  \  jk 
a  +  6  <  j  jk 


where  u  is  a  solution  of  v2(w  )  “  v^Cw  )  in  the  range  <  w  <  00 . 
other  details  the  reader  is  referred  to  Parfitt  and  Eringen  [1966]. 


For 
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25.  REFLECTION  OF  A  LONGITUDINAL  DISPLACEMENT  WAVE 


In  this  section  we  study  the  reflection  of  a  plane  longitudinal 

displacement  wave  at  a  stress-free  plane  surface  of  a  micropolar  half 

space.  If  x  *  0  is  the  plane  of  the  incident  displacement  vector,  then 

the  reflected  waves  from  the  boundary  surface  z  =  0  can  be  shown  to 

remain  in  this  plane  so  that  the  problem  is  a  two-dimensional  one.  On 

the  boundary  plane  z  =  0  being  free  from  traction  we  must  have 

t.  \  ■  m,  \  “0.  Since  we  also  have  u  =  <p_  “  <p_  =  0  through  (22.11)  and 

~(n)  -(n)  -  l  5 

(22.12)  for  the  tractions  on  z  =  0,  we  will  have  t.  .  =  t  ,  t  .  .  »  t 

(z)z  zz  (z)y  zy 

and  m.  .  *  m  ,  or  using  the  constitutive  equations  and  (24.4)  we  have 

(Z  J  X  Z X 


t  =  A(u  +  u  )  +  (2y  +  <)(u  -  U  )  =  0 

zz  »yy  »zz  >zz  x,yz 


(25.1)  t  =  u(u  -  U  )  +  (u  +  k)(u  +  U  )  +  <(i>  -  <J>  )  =  0 

zy  ,yz  x,yy  ,yz  x,zz  z,y  y,z 


m  =  yi®  -  $  )  “  0 

zx  z,yz  y ,zz 


The  nonvanishing  components  of  the  solution  vectors  (24.9)  appropriate  to 
this  problem  are  given  by 


ua  "  aa  €Xp  ^i(klna  ’  -  ~  U)lt^ 


(25.2) 


U6  "  iA8x  6XP  [i(k8n8  *  r  "  V>] 


*6  "  (V  +  B8zk)  6XP  li(k8n8  *  r  "  “6°  1 


where  =  k^v^(a  =  1,2;8  =  3,4)  and  the  repeated  indices  are  not  summed 
The  coefficients  A  and  B  are  related  to  each  other  by 

iu,0A3x 

B3  “  "  ~  2  2-2  2  (n3z^  "  n3yk) 

k3(v3  "  2tuok3  ~  c4) 


(25.3) 
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with  a  similar  equation  for  B^. 

The  potentials  listed  in  (25.2)  satisfy  the  boundary  conditions  (25. 1)^ 
at  z  *  0  if 


CO,  ■  U)_  ■  CO,  *  10 

13  4 


(25.4) 


k.n.  »  k.n.  ■  k.n.  »  k.n, 

1  ly  1  2y  3  3y  4  4y 


k.n,  “  k.n„  *  k_n_  =  k.n. 

1  lx  1  2x  3  3x  4  4x 


and 


(25.5)  [  Ak^  +  (2p  +  ^Jk^n^Ja^  +  [  Ak^  +  (2p  +  k)  z  ]  a2 


-  (2u  +  <)K3"3y"32A3x  -  (2v  +  <)k4n4yn42A4x  -  0 


Since  the  incident  wave  is  in  the  x  -  0  plane,  n^  *  0  and  (25. 4) ^  yields 


n„  *  n,  ■  n.  -  0 
2x  3x  4x 


which  is  the  proof  of  our  statement  that  all  reflected  waves  are  on  the 
x  -  0  plane.  Waiting 


n,  =  cos  0,  ,  n„  =  cos  0»  ,  n.  cos  0_  ,  n.  =  cos  0. 

ly  1  2y  2  3y  3  4y  4 


(see  Fig.  25.1)  and  co^  =  k^v^  from  (25.4)  we  get 


V3  V4 

(25.6)  cos  0„  *  cos  0,  ,  cos  0„  ■  —  cos  0,  ,  cosO  .  «  —  cos  0, 

2  1  3  v  1  4  v^  1 


where  v^  are  the  speeds  of  the  various  waves  found  in  Art.  24.  From  (25.6) 
it  is  clear  that  0^  ■ 
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Subs  Ci  tuting  (25.2)  into  the  remaining  two  equations  of  (25.1),  we 
obtain  two  other  equations.  These  two  equations  and  (25.5)  are  adequate 
to  determine  the  amplitude  ratios  a2/a^,  A^/a^  and  A^/a^.  The  amplitude 
B  is  determined  from  equations  (25.3)  and  similar  ones  for  . 

If  we  set  n^  =  0  in  (25.5),  we  find  that  a^  =  -a^ .  This  shows 
that  for  normal  incidence  the  reflected  longitudinal  displacement  wave 
is  also  normal  to  the  boundary.  For  n^  f  0,  the  following  are  the  solu¬ 
tions  for  the  ratios  of  the  wave  amplitudes 


2  2  2  2 
(25.7)  —=  {[A  +  (A+2y+<)  tan  6-^  ]  [  ( ir*"<)  tan  0^  -  y  -  (y+<)Q^  tan  0^  tan  0^ 


-  q: 


«  „  tan  0 .  .  „ 

+  (uQj+Qj)  ta~—  ']  +  (2yf<)  tan  0^^  tan  0^  (Q1~l)  } 


2  2  2 

x{-[A+  (A+2y+<)  tan  0  J  [  ( y+  < )  tan  0  -  u  -  (y+<)Q  tan  0,  tan  0/  -  y|^ 

1  *  i  J 


9  9  t£in  0  .  9  9  1 

+  (pQ1+Q3)  tan  Q  ]  +  (2p+<)  (Q^-l)  tan  tan  0^} 


/  9 

(25.8)  — —  =  -2(2y+<)[A  +  (A+2y+<)  tan  0,]  tan  8. 

al  1  1 


*{-[A  +  (A+2u+<)  tan2  © ^ ] t ( U+  < )  tan2  0^  -  u  -  (u+<)Q^  tan  tan  0^ 


r\  n  n  tan  0  •  n  n  . 

-  +  (MQ3+Q3)  —  J  ]  +  (2u+<)  (Q^-l)  tan  tan 


(25.9) 


2  2 

2  2  ^  2  2  2<u  -1 

Qi  [v3(1  -  Hr  '  c4]lv'(1  ’  "  c4] 

ui  w 


2  2  1 

W.  v,  „  2u>  _1 

—  ~  I^1  2^  "  c4] 

CJ  V.  0) 

4y 

2  2  „  2 
o)  v_  _  2uj  „  _l 

-  — ^  [v2(l - 2-)  -  c2  1 

2  2  1  4U  2  C4 

0)  V.  U) 

3y 
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These  results  were  studies  in  detail  by  Parfitt  and  Eringen  [1966]  for 
various  special  cases.  A  summary  of  their  findings  is  given  below: 

A  longitudinal  displacement  wave  at  a  plane  stress  free  boundary 
in  general  produces  three  reflected  waves  (as  compared  with  the  two  waves 
of  the  classical  theory) ;  one  longitudinal  displacement  wave  reflected 
at  an  angle  which  is  the  same  as  the  angle  of  incidence,  and  two  sets 
of  coupled  transverse  waves,  one  travelling  with  speed  v^  and  the  other 
set  with  speed  v^  (Fig.  25.1).  Their  angles  of  reflection  are  calculated 
by  (25.6).  For  normal  incidence  (0^  =  90°)  the  waves  at  speeds  v^  and  v^ 
vanish  and  the  reflected  wave  is  a  longitudinal  displacement  wave  normal 
to  the  boundary.  The  amplitude  ratios  as  functions  of  the  angle  of  inci¬ 
dence  6^  are  given  by  (25.7)  to  (25.9)  and  similar  equations  for  B^. 

The  general  solution  prevails  as  0^  decreases  from  90°  to  a  critical 
* 

angle  0^  which  makes  0^  =  0  and  A^  =  0.  At  this  angle  of  incidence  we 

have  a  surface  motion  travelling  at  speed  v^  and  a  reflected  longitudinal 

*  * 

wave  at  speed  v^  and  angle  *  ®1  =  ®1*  As  ®l  decreases  from  6-^  Co  zero, 
the  angle  0^  becomes  complex.  The  interpretation  in  this  case  is  that  a 
longitudinal  wave  is  reflected  into  the  medium  at  an  angle  0^  and  a  coupled 
transverse  wave  (decaying  with  depth  and  travelling  at  a  speed  c,  in  the 
range  v^  <_  c  <_  v^,  is  reflected.  For  0^  0,  a  limit  solution  is  possible 

Parfitt  and  Eringen  have  also  studied  the  reflection  of  coupled 
transverse  shear  and  microrotational  waves  and  the  reflection  of  a  longi¬ 
tudinal  micro rotational  wave.  For  these  and  other  interesting  results,  the 
reader  is  referred  to  the  foregoing  reference. 
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26.  MICROPOLAR  SURFACE  WAVES1 


In  this  section  we  investigate  the  propagation  of  surface  waves 
in  a  micropolar  half  space.  In  Section  25  we  have  shown  that  the  inci¬ 
dent  and  reflected  waves  propagate  in  the  same  plane.  We  select  this 
plane  to  be  the  x  ■  0  plane.  Thus  we  take  u^  =  *2  =  ^  =  ^ 

and  U£  =  v(y,z,t)  u^  =w(y,z,t)  and  =  4>(y  ,z  ,t)  as  functions  of  y,  z, 
and  t  only.  In  this  case  (22.1)  and  (22.2)  reduce  to 


(X  + 


\  J_  ( 9v 

^  3y  3y 


3w. 

3z; 


+  (u  + 


k)  V^v 


+  K 


3*  3  v 


(26.1) 


(X  + 


.  3  ,  3v 

u)  sT  (s7 


3w. 

3z; 


+  (p  + 


k)  V^w 


3  <p  3  w 


3y 


3t 


0 


+  <(|^  .  |X}  .  -  pj 


0 


where  V  is  the  two-dimensional  laplacian  operator,  e.g., 


We  consider  waves  which  are  propagating  in  the  plane  x  =  0  with  an  amplitude 
decay  in  the  z  direction. 

v  *  A  exp  (-4z)  exp  [iq(y  -  ct)  ] 

(26.2)  w  =  B  exp  (~4z)  exp  [iq(y  -  ct) ] 


*  ■  C  exp  (-;z)  exp  [ iq (y  -  ct)  ) 


1 


Suhubi  and  Eringen  [1964b] 
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Substituting  (26.2)  Into  (26.1)  gives  three  equations 

[-(A  +  2 u  +  <)q2  +  (y  +  <) ;2  +  pq2c2]A  -  (A  +  y)iq;B  -  kt;C  -  0 
(26.3)  -  (A  +  u)iq^V  +  [(A  +  2y  +  <)  42  -  (y  +  <)q2  +  pq2c2]B  -  iq<C  »  0 
<CA  +  iqKB  +  [y(t;2  -  q2)  -  2<  +  pjq2e2]C  -  0 


A  nonvanishing  solution  for  A,  B,  and  C  may  exist  if  the  determinant  of 
the  coefficients  is  zero.  This  gives 


c  c  «•  c  2 

(26. A)  [(e  +  -|);2  -  (e  +  - )q2]  {[j6;2  -  J(0  -  “)q2  -  2c]  x 

C2  C2  C2 

[e  +  IK2  -  (e  1  -  “jOq2 ]  +  e2(c2  -  q2)}  -  0 

C2 

where 


(26.5)  e  =  c2/c2  ,  0  =  c2/c2 

and  c^  to  c^  are  given  by  (24.3).  A  set  of  approximate  roots  of  (26.4)  is 

2 

obtained  by  neglecting  the  terms  containing  c  .  Hence 


4 


2 

1 


(1 


(26.6) 


[1  -  (1  -  c)  £j]q2 
C2 


4 


2 

3 


In  order  for  the  waves  to  be  surface  waves,  we  must  consider  only  the  posi¬ 
tive  values  of  the  roots  »  anc*  4-j*  The  displacement  and  microrotation 

field  can  now  be  expressed  as 
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v  -  £  A,  exp  (-:z)  exp  [iq(y  -  ct)] 

k-1 


(26.7) 


w  -  [  \ ^  exp  ("?kz)  exp  [iq(y  -  ct)  ] 

k-1 


where 

(26.8) 


-  w3A3  exp  (-;3z)  exp  (iq(y  -  ct)) 


q 

ii 

C3  ■ 


2  5, 


1  r/1  ,  \/2e  c2  2.  c2  2. 

7T  [(1+  e  (^e  2  q  )+“^] 

J  c4  c2 


On  stress-free  boundary  surface  we  must  have 


(26.9) 


t2z-il7+<x  +  2,,+  '0£ 


3w  .  ,  .  .  3v  . 

Czy  '  "17  +  (u  +  37  +  «♦ 


at  z  -  0 


m 


zx 


■  yii.O 

3z 


Substituting  (26.7)  into  (26.9),  we  obtain  a  set  of  three  homogeneous  equa¬ 
tions  for  A^,  and  Ay  The  determinant  of  the  coefficients  must  vanish. 
Hence 

(26.10)  u3  -  8o)2  +  8(3  -  k)u  -  16(1  -  k)  -  16e(l  -  ku>)  -  0 

(26.11)  c3u3  -  0 
where 

(26.12)  u  =  c2/c2  ,  k  =  c2/c2 

For  e  -  0,  (26.10)  reduces  to  the  classical  expression  of  the  Rayleigh  surface 

waves.  Denoting  the  values  of  u  for  this  case,  with  u  to  a  first-order 

R 
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appi animation  in  e,  we  obtain 


16(l-ku)  ) 

(2c  13)  to  ■  u_  +  — 5 -  c 

3ior  -  16toR  +  8(3-2k) 


For  k  *  1/3,  which  correspond  to  Poissons  ratio  1/A,  and  for  the  incom¬ 
pressible  solids  (k  ■  0)  we  have,  respectively, 


(26. 1A) 


c  -  0.919  (1  +  0.932e)c2  ,  (k  -  1/3) 
c  -  0.955  (1  +  0.783e)c2  ,  (k  -  0) 


The  tertrs  ccntaining  e  are  the  first-order  corrections  to  the  Rayleigh 
wave  velocity  in  each  case. 

Equation  (26.11)  gives  the  speed  of  propagation  of  a  new  type 
surface  wave  not  encountered  in  classical  elasticity.  This  is  given  by 


(26.15)  ~  =*  [2e/j  (1  +  e  -  6)]1/2q_1  =  (2e/j) 1/2q_1 

2 

This  new  wave  speed  depends  on  Jt  and  it  is  dispersive. 


-128- 


27.  STRESS  CONCENTRATION  AROUND  A  CIRCULAR  HOLE 

For  the  determination  of  the  stresses  in  a  plate  with  a  circular 

hole,  we  need  to  express  the  basic  equations  of  micropolar  elasticity 

in  plane  polar  coordinates  (r,  8),  Fig.  27.1.  For  the  equations  of 

balance  of  momenta,  parallel  to  classical  elasticity,  we  have 

2 

3trr  ,  X  3t8r  ,  V  '  t6e  ,  .  3  “r 

— +7  — + - 7 +pfr  ■°~ 


...  ..  3tre  l  nee  cre  +  ter  3  ue 

(27'1)  ~  +  7  — +  ;  °  8  ” 


3m  ,  3m  „  m 

rz  1  _ 0z  rz  .  .  „ 

3r  r  30  r  r8  8r  z  J 


d2<p 


3t 


For  the  nonzero  components  of  the  strain  tensor  given  by  (A. 17)  through 
the  methods  presented  in  Eringen  [1962,  Appendix}  we  obtain 

_  j^r 

err  3r 


(27.2) 


c00  r  ^  30  +  Ur 


m^e  _  . 

Cr 0  3r  z 


..  3u 

V  •  7  <1e  -  V  + 


The  constitutive  equations  for  the  stress  read 


t  *  (A  +  2u  +  k)c  +  X 
rr  rr  00 


'ee  ■  urr  +  (X  +  2u  +  *>cee 


cre  ■  (“  +  E)Ere  +  UEer 


t„  ■  pe  +  (p  +  k)c. 
0r  r0  0r 


(27.3) 
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Xn  these  equations  (trj.,  tr0,...),  (e  ,  and  uf  and  uQ  are, 

respectively,  the  physical  components  of  the  stresses,  strains,  and 
displacements . 

The  equations  of  compatibility  (8.9)  in  plane  polar  coordinates 
take  the  form 


3c0r  £6r  *  £r0  _  1  3err  _ 

3r  r  r  30  3r 


-  0 


(27.4) 


3e69  C06  ~  err  _  1  3er0  _  1  3*z 

3r  r  r  30  r  36 


3m 


6z 


m 


0z 


3r 


3m 

rz 

30 


For  static  problems  and  vanishing  body  loads  upon  introducing  the 
stress  functions  F(r,  0)  and  G(r,  0)  by 

1  3F  ;  1  32F  1  32G  1  3G 

rr  "  r  3r  2  ..2  "  r  3r30  2  30 

r  30  r 


32F  1  32G  _1_  ^G 

06  =  ,  2  r  3r30  2  30 

3r  r 


1  32F  1  3F  1_3G  1  32G 

Cr0  "  "  r  3r30  2  30  '  r  3r  '  2  .n2 

r  r  30 


1  32F  1  3F  32G 

C0r  "  "  r  3r36  2  30  .  2 

r  3r 


3G  1  3G 

rarz  3r  *  m0z  "  r  30 


we  see  that  equations  (27.1)  are  satisfied  identically.  From  (27.3),  we 
solve  for  the  strains  in  terms  of  the  stresses  and  then  substitute  (27.5) 
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for  the  stress  components.  If  we  now  use  the  expressions  (27.4),  we  get 


(27.6) 


where 


“•  (G  -  cVg)  -  -2(1  -  v)b2  “  (72F) 


7  3?  (G  "  cVG)  “  2(1  "  V)b2  £  ^ 


.2  -  Y(y  ±  <1 

k(2w  +  <) 


2  . 


2(2u  +  <) 


(27.7) 


V  5 - A -  ,2  5  _£  +  i  J.  +  i_  _£ 

2X  +  2tJ  +  k  lr!  t  >r  r2  se2 


Solutions  of  the  following  equations  are  also  solutions  of  (27.6) 


(27.8) 


2  2  2  2 
v  F  ■  0  ,  V*(G  -  cVg)  -  0 


For  the  problem  of  a  circular  hole  in  a  plate  (or  cylindrical  cavity  in  an 
infinite  solid)  subject  to  a  field  of  simple  tension  at  infinity,  an  appro¬ 
priate  solution  of  (27.8)  is 


(27.9) 


A 

F  ■  r2(l  -  cos  26)  +  log  r  +  (— j  +  A^)  cos  20 

r 


^4 

G  -  [— j  +  A^K^r/c)]  cos  29 
r 


where  K2  is  the  modified  Bessel  function  of  the  second  kind  and  second  order. 
This  set  satisfies  (27.6)  if 


(27.10) 


A.  -  8(1  -  v)c2A_ 
4  3 


The  remaining  four  constants  A^,  A2 ,  A^,  and  A,.  are  to  be  determined  from  the 
boundary  conditions 
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(27.11) 


t  ■  t  «  ■  m  "0  for  r  -  a 
rr  r9  rz 


t  ■  ^  (1  +  cos  2  0) 
rr  2 

for  r  ■  00 

tr6  -  -  f  sin  20 


m  ■  0 
rz 


Here  T  is  the  constant  tension  field  at  a  plane  x  -  const  at  infinity,  Fig. 
27.1.  By  use  of  (27.5)  and  (27.9),  we  calculate  the  components  of  the  stress 
and  couple  stress  censors 


rr 


T  A.  6A_  4A.  6A, 

—  (1  +  cos  20)  +  —  -  ( — jr  +  — — - — )  cos  2  0 

r  r  r  r 


2A  2 

+  K0(r/c)  +  (1  +  te-K^r/c)  ]  cos  20 

r 


■00 


T  0A»  6A . 

J  (1  -  cos  20) - j  -  — -fp  cos  20 

r  r  r 


(27.12)  t 

r0 


t 


0r 


m 


rz 


m 


0z 


-  “  K0(r/c)  +  (1  +  ■^y)K1  (r/ c)  ]  cos  20 

r 

6A_  2A_  6A, 

“(2+  2  '  4)  8in  26 

r  r  r 

A  2 

+  [“ ■  K0(r/c)  +  (1  +  i—)K1(r/c)]  sin  20 

r 

_  6A„  2A_  6A, 

"(2  +  4  +  2  ~  ~V  Sln  26 

r  r  r 

A-  ,2  -  . _  3 

+  -f  [(1  +  ~)K0(r/c)  +  (— ■  +  —-jK^r/c)]  sin  20 
c  r  r 

2A  A  2 

- j  sin  20  -  -f[^-  K0(r/c)  +  (1  +  ^j)K  (r/c)]  sin  20 

r  r 


2A,  2A 

{ — j  +  (K0(r/c)  +  —  (r/c) )}  cos  20 
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Using  the  boundary  conditions  (27.11),  we  find  that 

T  2  Ta  (1  ~  Fl) 

1  “  "  2  a  *  A2  “  4(1  +  F1) 


(27.13) 


2  2  2 
Ta  4(1  -  v)a  b  T 

A3  "  2(1  +  Fx)  A4  1  +  Fx 
TacF^ 

A5  “  ‘  (1  +  F1)K1(a/c) 


where 


(27.14) 


.2  2 

8(1  -  v)  \  [4  +  + 

c  c 


2a  K0(a/c) 
c  K1(a/c)^ 


Substituting  (27.13)  into  (27.12),  we  obtain  the  stress  and  couple  stress 
fields.  The  value  of  toa  at  the  perifery  of  the  circular  hole  is  of  great 
interest.  For  this  we  obtain 

(27.15)  tM  -  1(1  + 


The  maximum  value  of  this  t„„  occurs  at  0  ■  iv/2. 

66max 


(27.16) 


tee  max  .  3  *  F1 

T  =  c  “  1  +  F1 


The  quantity  Sc  so  defined  is  the  stress  concentration  factor.  From  (27.14) 

it  is  clear  that  Sc  depends  on  v,  a,  b,  and  c. 

The  above  result,  (27.16),  was  given  by  Kaloni  and  Ariman  [1967]  who 

adopted  the  solution  of  the  same  problem  for  the  indeterminate  couple  stress 

2  2 

theory  given  by  Mindlin  and  Tiersten  [1962],  If  we  set  b  /c  *  1,  we  obtain 
for  F^  the  one  given  by  Mindlin  [1963],  namely, 


2  2 

where  H  is  a  physical  length  for  a  given  material.  The  case  of  b  /c  -  1 
gives  <  ■  -2\i  which  is  not  acceptable  on  grounds  of  uniqueness  and  physical 
reasonableness  since  <  cannot  be  as  great  as  twice  the  shear  modulus.  Other 
discussion  of  these  results  and  comparison  of  the  stress  concentration  factors 
of  the  indeterminate  couple  stress  theory  and  micropolar  theory  are  to  be 
found  in  Kaloni  and  Ariman  [1967],  Below  we  reproduce  several  of  their  curves. 
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28.  GALERKIN  AND  PAPKOVITCH  REPRESENTATIONS1 

A  useful  mathematical  representation  for  the  displacement  and  micro¬ 
rotation  fields  can  be  constructed  as  follows:  Let 

Xi  ^  ■  T  ^  -k  •  «  5  XJ  +  X2  +  X3  '  Q1  E  (i  +  2“  +  “>«  -  ^ 

(28.1) 

Q2  =  (W  +  <)Q  -  PT2,Q3  =  (W-&fY)Q-pjT2-2<,  Q4=  yQ  -  pjT2  -  2< 

and  let  the  matrix  L«  ^ijl  '  ”  1*2,.  ..,6)  be  given  by 

(28.2) 


q2  +  (*+W)X2 

(A+P  )X1X2 

(x-*-p)xxx3 

0 

-<X3 

<x2 

(X+p)X1X2 

Q2  +  (A+w)X2 

(x+p)x2x3 

<x3 

0 

-<xi 

(A+y)X  x3 

(x+p)x2x3 

Q2  +  (X+P)X2 

_ICX2 

<X1 

0 

0 

-kX3 

<X2 

Q4  +  (o+B)X2 

(o+6)X1X2 

(o+6)X1X 

<X3 

0 

-kXl 

(o+BJX^ 

Q4  +  (cH-6)X2 

(o+6) X2X 

-<X2 

0 

(o+B)  X^ 

(o+6)X2X3 

Q4+(ct+6)X: 

The  equations  (22.1)  and  (22.2)  can  be  expressed  in  the  matrix  form 


1 


Sandru  [1966] 
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(28.3) 


■  U1 

1 

h-» 

_ 1 

U2 

f2 

U3 

.  -p 

f3 

♦l 

*1 

♦2 

*2 

.  '3  . 

.  ‘3  . 

For  the  Inverse  of  the  matrix  we  have,  formally, 


N, 


(28.4) 


L-l.  _  i.1  .. 

ij  QiQ3(Q2Q4  +  <^Q) 


where 


2, .,2. 


Nii  ’  VQ1Q4  ‘  l <^+W>Q4  -  <  ]XJ}  ,  (i  -  i» 2, 3) 


Nii  "  Q1{Q2Q3  *  l<a+6>Q2  ‘  K2lXi-3}'  (i  "  4’5’6) 


Nij  "  Nji  "  Q3t“CX+y)Q4  +  k  ]XAX  ,  i  i  j,  (i.j  -  1,2, 


3) 


(28.5) 


Nij  "  Nji  *  Q1l-(oH-6)Q2  +  <  ]Xi_3Xj_3,  i  4  j,  (i.j  -  4,5,6) 


W  m  N  m  N  •  N  »  M  +  N  >0 

14  25  36  W41  52  w63 


N15  -  ‘Si  '  N42  “  -N24  “  “Ws 


N16  -  -N61  “  N43  ■  -N34  "  ~*W2 


N26  "  -N62  "  N53  "  -N35  "  *Wl 


Consider  now 
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(28.6) 


'  U1  ■ 

F‘  1 

u2 

CM 

u3 

-  N 

F3 

* 

♦l 

F1 

*2 

F* 

2 

.  *3  . 

F* 

L  3  J 

If  ^  *  Q^F  and  #2  “  Q^F  through  (28.5)  and  (28.6),  we  find  that 


(28.7) 


u(x,t)  -  Q1  □4$1  -  [(X  +  u)04  -  <2]?V  •  — k  D3V  x  $2 
$(x,t)  -  □2n3?2  -  [(a  +  6)D2  -  <2)VV  •  -kDx7  x  ^ 


where 


D.  =  (X  +  2p  +  <)7  -  p  — r  ,  □  ,  =  (y  +  <)V  -  p  — - 

X  3t  3t 

(28.8) 

2  2 

□  =  (a  +  6  +  Y)V2  -  pj  -i-  -  2< ,  □  i  yV2  -  pj  -2< 

3t  H  3t 


From  (28.3),  (28.4),  and  (28.7),  it  follows  that  <I>3  and  {>2  satisfy  the  follow¬ 
ing  uncoupled  equations 


(28.9) 


□  x( D2  +  <2v2)$x  «  -pf 


□  3(D2D4  +  <V)?2  -  -pf 


If  we  take  <  -  0  in  (28.9)^,  we  obtain  a  representation  known  in  the  classical 
theory  of  elasticity. 

In  the  static  case,  we  set  T  ■  0  in  (28.7)  and  (28.8)  and  obtain  the 
Galerkin  representation  for  the  micropolar  elasticity,  namely, 


(28.10) 


u(x)  *  (X  +  2 vi  +  k)72(y72  -  2k)$^  -  [ y ( X  +  y)72  -  «(2A  +  2y  +  <)] 

2 

77  • ^  -  <[(o  +  6  +  Y>v  "  2k]7  x  *2 

$(x)  *  (y  +  k)72[(o  +  8  +  y)V2  -  2«)f2  -  [(y  +  k)  (a  +  B)72  -  k2]77  •  4>2 
-  tc  ( X  +  2y  +  k)  72  (7  x  j  ) 

where  and  $2  satisfy  the  equations 

(A  +  2y  +  <)7^[(y  +  k)y72  -  <(2y  +  tc)]$^  *  -pf 

(28.11) 

[(a  +  8  +  y)^2  “  2< ] [ (y  +  k)y7^  -  <(2y  +  <)]?2  *  “P£ 

For  k  ■  0.  (28. 10) ^  gives  Galerkin's  representation  of  classical  elasticity. 

We  decompose  the  body  force  and  body  couple  into  irrotational  and 
solenoidal  parts  as  follows: 


(28.12) 


pf  ■  7IIq  +  7  x  n 

*  * 

pi.  •  vn o  +  7  x  n 


and  note  that 

+  k272)$1  -  7A0  ,  EU^  -  7  x  A 

(28.13) 

(□2CV  <2V2)?2  -  7A*  ,  D^2  -  V  x  A* 

From  (28.7)  and  (28.9)  it  now  follows  that 

(28.14)  O^Ao  “  - H 0  (  □3aJ  -  -n* 


and  we  obtain 
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(28.15) 


u  -  7A0  +  7  x  (  □  A)  -  <V  x  (7  x  A*) 
$  -  -t(V  x  (V  x  A)  +  7  x  (  □  2A*) 


Thus,  if  we  determine  Ao,  A  and  A  for  a  given  problem,  the  displacement 
and  microrotation  fields  can  be  calculated  from  (28.15). 

Another  useful  decomposition  can  be  deduced  from  (28.11)  with 
f  ■  4  «  Q  by  selecting  4^  and  $2  as  the  sum  of  three  special  vector 
functions.  This  has  the  form 


u 


Ax  +  A2 


21+  2u+< 
A(l+2y+<) 


V(x  • 


+  Aq) 


fell 

k(2w+ic) 


+  a2) 


(28.16) 


+  2V  x  B.  -  r-T~  ^  x  B_ 
-1  2p+<  -  -3 


$  -  vv 


'2  + 


B  -  -  luteh.  77 

-3  <(2u+<) 


-3  + 


VV 


?i  +  ?? 


<  A.  t  Jj+£ 

'1  K 


V  x  A2 


where  A0,  A^,  A2 ,  ,  B2,  and  B^  satisfy 

v2a0-o.  72ax.o,  (1  -  v2)A2  -  0 


(28.17) 


'B  -  0  (1  -  a  *  ^  +  K  V^)B  ■  0  (1  -  (V- ■  V^)B  *» 

'1  '•  u  2<  '2  '»  u  <(2y  +  <)  '3 


These  results  when  ic  -*■  <*>  reduce  to  the  corresponding  decomposition  for  the 
indeterminate  couple  stress  theory  obtained  by  Mindlin  and  Tiersten  [1962]. 
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29.  MICROPOLAR  INFINITE  SOLID  SUBJECT  TO 
CONCENTRATED  FORCE  AND  CONCENTRATED  COUPLE 


The  problem  of  an  infinite  solid  subject  to  a  concentrated  force  F 
and  a  concentrated  couple  C  at  the  origin  of  coordinates  x  ■  0  it;  of  fund¬ 
amental  interest.  In  classical  theory,  this  problem  is  known  as  Kelvin's 
problem.  This  problem  for  the  static  case  was  treated  by  Sandru  [1966]. 

(i)  Concentrated  Force.  Let  F  be  a  concentrated  force  acting  at 
the  origin  of  the  coordinates  x  ■  0.  We  write 


(29.1) 


pf  -  F5(x) 


where  6(x)  is  the  Dirac  delta  function. 

The  solution  of  (28.12)^  is  given  by 


(29.2) 


n0(x)  -  -  —  /pf(pY(~)dv(p 
n(x)  -  -  ~  ^ Pf(p  *  ?(7)dv(p 


where 


(29.3)  r  =  [(Xl  -  q)2  +  (x2  -  52)2  +  (x3  -  C3)2] 


1/2 


Using  (29.1)  we  find 


(29.4) 


?  ■  -  -t  ?  ‘  ?<?> 


where 


2  2  2 
R  £  (Xj.  +  x2  +  x3) 


1/2 


(29.5) 
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From  (28.13).  _  for  the  static  case  we  get 
1  ,  / 


(29.6) 


(A  +  2w  +  k)V2A0  -  ^  F  •  7(i) 


[y(p  +  < -  <(2y  +  <)  ] A  -  —  F  x  V(^) 


These  equations  possess  the  solutions 

1 


A0  “ 


F  .  x 
8tt(A+2v+<)  R 


(29.7) 


-  "  8ttk(2ij+>:)  -  *  *?R*  +  4tt k ( 2 u-K}  -  “  (1 


e-R/l)l 


where 


(29.8) 


,2  =  Y(w  ±  <) 

*  ‘  (2w  +  k)< 


With  A  -  0,  (28.15)  gives  for  the  static  case 


5A+6y+3* 


F 
r  + 


2A+2u+< 


F  •  x 


8tt (2u+k)  (A+2u+x)  R  8tt(2vi+k)  (A+2p+<)  ^3 


(29.9) 


4tt(2u+<) 


1 - 7  V  x  {7  x  [I  (e“RA-l)  ]  } 


*  - 


47r(2tJ+K)  ' 


V  x 


tf  (1  -  e'R/S] 


For  k  =  y  ®  0,  (29.9).,  reduces  to  the  well-known  solution  of  Kelvin's  problem. 
Love,  [1944,  p.  185 J.  The  solution  of  this  problem  for  the  indeterminate 
couple  stress  theory  is  obtained  by  letting  k  -*  «.  This  gives  the  result 
obtained  by  Mindlin  and  Tiersten  [1962], 

(ii)  Concentrated  Couple.  Let  M  be  a  concentrated  couple  acting  at 
x  ■  0.  We  write 


pi  -  M6(x) 
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In  this  case,  we  express  the  solution  of  (27.12)  in  the  same  way  as  in  (28.2). 
Using  (28.13)2  3  ^or  static  ca9e  we  have 


(29.10) 


where 


((a  +  8  +  y)V2  -  2k ] A*  -  -nj 


[ Y ( W  +  <)7^  -  <(2p  +  k)V2]A  ■  -II 


(29.11) 


-  -  47  M  •  7(i)  ,  n*  -  -  i  M  X  7(i, 


Equations  (29.10)  have  the  solutions 


(29.12) 


where 


a*  1  xi  nr1  (  “R/h  ,vi 

A°  '  8^  ‘  (e  -  1)1 


arsko.  ?  -  <^> +  -(-2^-  ?  -  £ (1  -  «'**» 


(29.13) 


.  2  _  a  +  8  +  Y 
h  =  2^ - 


Substituting  these  results  into  (28.15)  for  the  static  case,  and  with  Aq  ■  0 
and  A  »  0,  we  obtain 


4tt(2u+<)  ~ 


M  _d  /  0 

l  *  It  (1  -  e  K/*)] 


(29.14) 


*  -  8^  <e"R/h-l)])  8-  a.fc*  ).  ;  >  <?  «tf  (1  -  e'R/t)l) 


These  results  may  be  used  in  (20.23)  and  (20.24)  to  obtain  the  stress  and 
couple  stress  fields. 

It  is  worth  noting  here  that  the  concentrated  couple  is  a  fundamental 
problem  which  is  not  deduced  as  the  limiting  solution  of  two  equal  parallel 
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forces  directed  in  opposite  sense.  We  note  that  in  micropolar  elasticity, 
the  concept  of  body  couple  is  totally  independent  of  the  force  and  it  can 
exist  even  when  the  body  force  is  absent.  For  micropolar  elasticity,  there¬ 
fore,  force  and  moment  singularities  will  have  totally  different  natures. 
This  is  then  expected  to  affect  uniqueness  theorems  for  infinite  and  for 
the  finite  regions  considerably. 
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SUMMARY 

The  theory  of  micropolar  elasticity  presented  here  should  provide 
an  adequate  background  for  pursuing  analytical  work  and  starting  a  badly 
needed  experimental  program.  We  believe  the  theory  is  well-posed  with  field 
equations,  boundary  and  initial  conditions.  Certain  wide  classes  of  unique¬ 
ness  theorems  have  been  proved  though  not  presented  here  (cf.  Eringen  [1966a]). 
The  important  implications  of  the  theory  are  brought  to  the  surface  especially 
in  connection  with  the  problems  in  the  field  of  wave  propagation.  Existence 
of  additional  waves  over  those  existing  in  classical  elasticity  should  be 
attractive  to  workers  in  the  field  of  experimental  wave  propagation.  The 
linear  theory  is  simple  enough  to  lend  itself  to  the  solution  of  non-trivial 
boundary  and  initial  value  problems.  With  the  additional  Internal  degrees 
of  freedom  provided  by  the  microrotation  and  spin  inertia,  it  incorporates 
the  problems  involving  concentrated  body  and  surface  couples  into  the  funda¬ 
mental  singularities  of  the  field. 

The  field  is  rather  new  and  not  even  partly  explored.  Experimental 
works  are  badly  needed.  Nevertheless,  it  seems  to  us  that  the  logical  founda¬ 
tion  of  the  theory  is  solid  and  promising  for  the  understanding  of  the  mechanics 
of  micropolar  solids. 
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POSSIBLE  DIRECTIONS  FOR  FUTURE  RESEARCH 

The  theory  of  micropolar  elasticity  presented  here  is  one  of  the 
simplest  extensions  of  the  classical  theory  of  elasticity  for  the  treat¬ 
ment  of  materials  with  microstructure.  Extensions  of  this  theory  to 
micropolar  fluids  and  viscoelasticity  exist  (cf.  Eringen  [1964],  [1966b], 
[1967]).  A  nonlinear  theory  is  also  contained  in  the  more  general  theory 
of  microelasticify  given  by  Eringen  and  Suhubi  [1964a, b].  The  theory  of 
micromorphic  materials,  of  which  the  theory  of  microelasticity  is  a  repre¬ 
sentative  field  for  the  oriented  solids,  possesses  promise  for  entering 
the  granular  and  molecular  world  of  materials  from  the  continuum  side.  It 
should  not  be  surprising  if  such  a  theory  would  be  in  wide  use  for  the  full 
description  of  material  properties  of  composites,  granular  and  fiberous 
solids.  The  connection  of  these  theories  to  continuum  dislocation  theory 
has  already  been  exhibited.  The  intimate  ties  between  the  plasticity  theory 
and  continuum  dislocations  has  been  recognized  for  some  time  by  seme  research 
workers  in  this  field  (cf.,  Kondo  [1962],  [1963],  Kroner  [1963],  Bilby  [1960], 
Bilby,  Gardner  and  Stroh  [1967]).  Existing  theories  of  micTomechanics ,  multi¬ 
polar  theories,  and  the  continuum  theory  of  dislocations  have  not  been 
amalgamated  into  a  unified  structure  as  yet,  although  some  attempts  exist 
in  this  direction.  Presently,  serious  efforts  are  being  made  to  bring  some 
order  to  the  world  of  the  microcontinuum.  The  field  of  micromorphic  materials 
needs  and  deserves  attention  from  both  theoretical  and  experimental  workers. 
Especially,  the  need  for  rational  experiments  is  felt  badly.  The  theory  of 
micropolar  elasticity  is  certainly  ready  for  such  ^  test  The  theory  of 
micromorphic  materials  with  its  logical  structure  ^nd  wide  possible  applica¬ 
tions  opens  new  and  promised  rich  lands  for  futu-e  workers. 
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